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Abstract. Quasinormal modes are eigenmodes of dissipative systems. Pertur- 
bations of classical gravitational backgrounds involving black holes or branes nat- 
urally lead to quasinormal modes. The analysis and classification of the quasinor- 
mal spectra requires solving non-Hermitian eigenvalue problems for the associated 
linear differential equations. Within the recently developed gauge-gravity dual- 
ity, these modes serve as an important tool for determining the near-equilibrium 
properties of strongly coupled quantum field theories, in particular their trans- 
port coefficients, such as viscosity, conductivity and diffusion constants. In as- 
trophysics, the detection of quasinormal modes in gravitational wave experiments 
would allow precise measurements of the mass and spin of black holes as well as 
new tests of general relativity. This review is meant as an introduction to the 
subject, with a focus on the recent developments in the field. 
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1. Introduction 



"The mathematical perfeetness of the black holes of Nature is [...] revealed at every 
level by some strangeness in the proportion in conformity of the parts to one another 
and to the whole. " S. Chandrasekhar, "The Mathematical Theory of Black Holes" 

Characteristic modes of vibration are persistent in everything around us. They 
make up the familiar sound of various musical instruments but they are also 
an important research topic in such diverse areas as seismology, asteroseismology, 
molecular structure and spectroscopy, atmospheric science and civil engineering. All 
of these disciplines are concerned with the structure and composition of the vibrating 
object, and with how this information is encoded in its characteristic vibration modes: 
to use a famous phrase, the goal of studying characteristic modes is to "hear the 
shape of a drum" [1]. This is a review on the characteristic oscillations of black 
holes (BHs) and black branes (BHs with plane-symmetric horizon), called quasinormal 
modes (QNMs). We will survey the theory behind them, the information they carry 
about the properties of these fascinating objects, and their connections with other 
branches of physics. 

Unlike most idealized macroscopic physical systems, perturbed BH spacetimes 
are intrinsically dissipative due to the presence of an event horizon. This precludes 
a standard normal-mode analysis because the system is not time-symmetric and the 
associated boundary value problem is non-Hermitian. In general, QNMs have complex 
frequencies, the imaginary part being associated with the decay timescale of the 
perturbation. The corresponding eigenfunctions are usually not normalizable, and, 
in general, they do not form a complete set (see Refs. [21 E] for more extensive 
discussions). Almost any real- world physical system is dissipative, so one might 
reasonably expect QNMs to be ubiquitous in physics. QNMs are indeed useful in 
the treatment of many dissipative systems, e.g. in the context of atmospheric science 
and leaky resonant cavities. 

Two excellent reviews on BH QNMs [4j[5] were written in 1999. However, much 
has happened in the last decade that is not covered by these reviews. The recent 
developments have brought BH oscillations under the spotlight again. We refer, in 
particular, to the role of QNMs in gravitational wave astronomy and their applications 
in the gauge-gravity duality. This work will focus on a critical review of the new 
developments, providing our own perspective on the most important and active lines 
of research in the field. 

After a general introduction to QNMs in the framework of BH perturbation 
theory, we will describe methods to obtain QNMs numerically, as well as some 
important analytic solutions for special spacetimes. Then we will review the QNM 
spectrum of BHs in asymptotically flat spacetimes, asymptotically (anti-)de Sitter 
(henceforth AdS or dS) spacetimes and other spacetimes of interest. After this general 
overview we will discuss what we regard as the most active areas in QNM research. 
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Schematically, we will group recent developments in QNM research into three main 
branches: 

(i) AdS/CFT and holography. In 1997-98, a powerful new technique known as the 
AdS/CFT correspondence or, more generally, the gauge-string duality was discovered 
and rapidly developed [6]. The new method (often referred to as holographic 
correspondence) provides an effective description of a non-perturbative, strongly 
coupled regime of certain gauge theories in terms of higher-dimensional classical 
gravity. In particular, equilibrium and non-equilibrium properties of strongly coupled 
thermal gauge theories are related to the physics of higher-dimensional BHs and 
black branes and their fluctuations. Quasinormal spectra of the dual gravitational 
backgrounds give the location (in momentum space) of the poles of the retarded 
correlators in the gauge theory, supplying important information about the theory's 
quasiparticle spectra and transport (kinetic) coefhcients. Studies of QNMs in the 
holographic context became a standard tool in considering the near-equilibrium 
behavior of gauge theory plasmas with a dual gravity description. Among other 
things, they revealed the existence of a universality of the particular gravitational 
frequency of generic black branes (related on the gauge theory side to the universality 
of the viscosity-entropy ratio in the regime of infinitely strong coupling), as well as 
intriguing connections between the dynamics of BH horizons and hydrodynamics . 
The duality also offers a new perspective on notoriously difficult problems, such as 
the BH information loss paradox, the nature of BH singularities and quantum gravity. 
Holographic approaches to these problems often involve QNMs. This active area of 
research is reviewed in Section \8\ 

(ii) QNMs of astrophysical black holes and gravitational wave astronomy. 

The beginning of LIGO's first science run (SI) in 2002 and the achievement of design 
sensitivity in 2005 marked the beginning of an era in science where BHs and other 
compact objects should play a prominent observational role. While electromagnetic 
observations are already providing us with strong evidence of the astrophysical reality 
of BHs [8] , gravitational wave observations will incontrovertibly show if these compact 
objects are indeed rotating (Kerr) BHs, as predicted by Einstein's theory of gravity. 
BH QNMs can be used to infer their mass and angular momentum [S] and to test 
the no-hair theorem of general relativity [lOl [11]. Dedicated ringdown searches in 
interferometric gravitational wave detector data are ongoing [121 E] • The progress on 
the experimental side was accompanied by a breakthrough in the numerical simulation 
of gravitational wave sources. Long-term stable numerical evolutions of BH binaries 
have been achieved after 4 decades of efforts [TH [TSl [H] , confirming that ringdown 
plays an important role in the dynamics of the merged system. These developments 
are reviewed in Section [9l 

(iii) Other developments. In 1998, Hod suggested that highly-damped QNMs 
could bridge the gap between classical and quantum gravity [17] . The following years 
witnessed a rush to compute and understand this family of highly damped modes. 
The interest in this subject has by now faded substantially but, at the very least. 
Hod's proposal has contributed to a deeper analytical and numerical understanding 
of QNM frequencies in many different spacetimes, and it has highlighted certain 
general properties characterizing some classes of BH solutions. These ideas and other 
recent developments (including a proposed connection between QNMs and BH phase 
transitions, the QNMs of analogue BHs, the stability of naked singularities and its 
relation with the so-called algebraically special modes) are reviewed in Section [TOl 

The present work is mostly intended to make the reader familiar with the new 
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developments by summarizing the vast (and sometimes confusing) bibliography on the 
subject. We tried to keep the review as self-contained as possible, while avoiding to 
duplicate (as far as possible to preserve logical consistency) material that is treated 
more extensively in other reviews on the topic, such as Refs. [H IHl [21 [TS] . A detailed 
understanding of BH QNMs and their applications requires some specialized technical 
background. QNM research has recently expanded to encompass a very wide range of 
topics: a partial list includes analogue gravity, alternative theories of gravity, higher- 
dimensional spacetimes, applications to numerical relativity simulations, explorations 
of the gauge-gravity duality, the stability analysis of naked singularities and ringdown 
searches in LIGO. Because of space limitations we cannot discuss all of this material in 
detail, and we refer the reader to other reviews. Topics that are treated in more detail 
elsewhere include: (1) a general overview of gravitational radiation [T9', '20' and its 
multipolar decomposition [2l]; (2) BH perturbation theory [22, 23, 24, 25, 26, 27,; (3) 
the issue of quantifying QNM excitation in different physical scenarios (see e.g. for 
an introduction pre-dating the numerical relativity breakthroughs of 2005, and [5^ for 
a more updated overview of the field); (4) tests of general relativity and of the no- hair 
theorem that either do not make use of ringdown [30] , or do not resort to gravitational 
wave observations at all [311 [32l [33] ; (5) BH solutions in higher dimensions [34] ; (6) 
many aspects of the gauge-gravity duality [T] [35] [36l [37l [38]. The reviews listed 
above provide more in-depth looks at different aspects of QNM research, but we tried 
to provide concise introductions to all of these topics while (hopefully) keeping the 
presentation clear and accessible. 

Chandrasckhar's fascination with the mathematics of BHs was due to their 
simplicity. BHs in four-dimensional, asymptotically flat spacetime must belong to the 
Kerr-Newman family, which is fully specified by only three parameters: mass, charge 
and angular momentum (see e.g. Ref. [39], or Carter's contribution to Ref. pQj). 
One expresses this by saying that BHs have no hair (or more precisely, that they 
have three hairs). A consequence of the no- hair theorem is that all perturbations 
in the vicinities of a BH must decay to one and the same final state, i.e. that all 
hairs (except three) must be lost. Perturbative and numerical calculations show that 
the hair loss proceeds, dynamically, via quasinormal ringing. The gravitational wave 
signal from a perturbed BH can in general be divided in three parts: (i) A prompt 
response at early times, that depends strongly on the initial conditions and is the 
counterpart to light-cone propagation; (ii) An exponentially decaying "ringdown" 
phase at intermediate times, where QNMs dominate the signal, which depends entirely 
on the final BH's parameters; (iii) A late-time tail, usually a power-law falloff of the 
field [41]. Mathematically, each of these stages arises from different contributions to 
the relevant Green's function (see Section [3T2|l . QNM frequencies depend only on the 
BH's parameters, while their amplitudes depend on the source exciting the oscillations. 

Numerical and analytical analysis of processes involving BHs confirm these 
expectations. QNMs were observed for the first time in numerical simulations of the 
scattering of Gaussian wavepackets by Schwarzschild BHs in 1970, soon after the BH 
concept itself was introduced and popularized by John Wheeler. Vishveshwara [42] 
noticed that the waveform at late times consists of a damped sinusoid, with ringing 
frequency almost independent of the Gaussian's parameters. Ringdown was observed 
again in the linearized approximation to the problem of a test mass falling from infinity 
into a Schwarzschild BH [33]. By now, decades of experience have shown that any 
event involving BH dynamics is likely to end in this same characteristic way: the 
gravitational wave amplitude will die off as a superposition of damped sinusoids. 
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Figure 1. Four different physical processes leading to substantial quasinormal 
ringing (see text for details). With the exception of the infalling-particle case 
(where M is the BH mass, the particle's mass and 1/12 the Zerilli wavefunction) , 
^22 is the I = m = 2 multipolar component of the Weyl scalar >I'4, M denotes 
the total mass of the system and r the extraction radius (see e.g. Ref. [44]). 



Figure [T] shows four different processes involving BH dynamics. In all of them, 
quasinormal ringing is clearly visible. The upper- left panel (adapted from Ref. 44J ) 
is the signal from two equal-mass BHs initially on quasi-circular orbits, inspiralling 
towards each other due to the energy loss induced by gravitational wave emission, 
merging and forming a single final BH [14L The upper-right panel of Fig. [T] shows 
gravitational waveforms from numerical simulations of two equal-mass BHs, colliding 
head-on with v/c = 0.94 in the center-of-mass frame: as the center-of-mass energy 
grows (i.e., as the speed of the colliding BHs tends to the speed of light) the waveform 
is more and more strongly ringdown-dominated [45j . The bottom- left panel shows 
the gravitational waveform (or more precisely, the dominant, I = 2 multipole of 
the Zerilli function) produced by a test particle of mass falling from rest into a 
Schwarzschild BH [33j: the shape of the initial precursor depends on the details of 
the infall, but the subsequent burst of radiation and the final ringdown are universal 
features. The bottom-right panel (reproduced from Ref. 46J) shows the waves emitted 
by two massive neutron stars (NSs) with a polytropic equation of state, inspiralling 
and eventually collapsing to form a single BH. 

QNM frequencies for gravitational perturbations of Schwarzschild and Kerr BHs 
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have been computed by many authors. Rather than listing numerical tables of well- 
known results, we have set up a web page providing tabulated values of the frequencies 
and fitting coefficients for the QNMs that are most relevant in gravitational wave 
astronomy |47| . On this web page, we also provide Mathematica notebooks to compute 
QNMs of Kerr and asymptotically AdS BHs [47] . 

1.1. Milestones 

QNM research has a fifty-year-long history. We find it helpful to provide the reader 
with a "roadmap" in the form of a chronological list of papers that, in our opinion, 
have been instrumental to shape the evolution of the field. Our summary is necessarily 
biased and incomplete, and we apologize in advance for the inevitable omissions. A 
more complete set of references can be found in the rest of this review. 

• 1957 - Regge and Wheeler [48| analyze a special class of gravitational perturbations 
of the Schwarzschild geometry. This effectively marks the birth of BH perturbation 
theory a decade before the birth of the BH concept itself. The "one-way membrane" 
nature of the horizon is not yet fully understood, and the boundary conditions of the 
problem are not under control. 

• 1961 - Newman and Penrose [49] develop a formalism to study gravitational radiation 
using spin coefficients. 

• 1963 - Kerr [SD] discovers the mathematical solution of Einstein's field equations 
describing rotating BHs. In the same year, Schmidt identifies the first quasar ("quasi- 
stellar radio source"). Quasars (compact objects with luminosity ~ 10^^ that of our 
sun, located at cosmological distance [51]) are now believed to be supermassive BHs 
(SMBHs), described by the Kerr solution. 

• 1964 - The UHURU orbiting X-ray observatory makes the first surveys of the X-ray 
sky discovering over 300 X-ray "stars" , most of which turn out to be due to matter 
accreting onto compact objects. One of these X-ray sources, Cygnus X-1, is soon 
accepted as the first plausible stellar- mass BH candidate (see e.g. [52 ). 

• 1967 - Wheeler [131 154] coins the term "black hole" (see the April 2009 issue of 
Physics Today, and Ref. [55] for a fascinating, first-person historical account). 

• 1970 - Zerilli [5B^, "37] extends the Regge- Wheeler analysis to general perturbations 
of a Schwarzschild BH. He shows that the perturbation equations can be reduced to a 
pair of Schrodinger-like equations, and applies the formalism to study the gravitational 
radiation emitted by infalling test particles. 

• 1970 - Vishveshwara |42j studies numerically the scattering of gravitational waves 
by a Schwarzschild BH: at late times the waveform consists of damped sinusoids (now 
called "ringdown waves" ) . 

• 1971 - Press [58] identifies ringdown waves as the free oscillation modes of the BH. 

• 1971 - Davis et al. [43] carry out the first quantitative calculation of gravitational 
radiation emission within BH perturbation theory, considering a particle falling 
radially into a Schwarzschild BH. Quasinormal ringing is excited when the particle 
crosses the maximum of the potential barrier of the Zerilli equation, which is located 
at r ~ 3M (i.e., close to the unstable circular orbit corresponding to the "light ring"). 

• 1972 - Goebel [59] points out that the characteristic modes of BHs are essentially 
gravitational waves in spiral orbits close to the light ring. 

• 1973 - Teukolsky [601 decouples and separates the equations for perturbations in 
the Kerr geometry using the Newman-Penrose formalism [49 . 

• 1974 - Moncrief [61J introduces a gauge-invariant perturbation formalism. 
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• 1975 - Chandrasekhar and Detweiler |62j compute numerically some weakly damped 
characteristic frequencies. They prove that the Regge- Wheeler and Zerilli potentials 
have the same spectra. 

• 1978 - Cunningham, Price and Moncrief [S^llMllSSj study radiation from relativistic 
stars collapsing to BHs using perturbative methods. QNM ringing is excited. 

• 1979 - Gerlach and Sengupta give a comprehensive and elegant mathematical 
foundation for gauge- invariant perturbation theory ^66. ,67] . 

• 1983 - Chandrasekhar's monograph [22] summarizes the state of the art in BH 
perturbation theory, elucidating connections between different formalisms. 

• 1983 - York [68] attempts to relate the QNM spectrum to Hawking radiation. To 
our knowledge, this is the first attempt to connect the (purely classical) QNMs with 
quantum gravity. 

• 1983 - Mashhoon [69] suggests to use WKB techniques to compute QNMs. Ferrari 
and Mashhoon [7^ analytically compute QNMs using their connection with bound 
states of the inverted BH effective potentials. 

• 1985 - Stark and Piran |71| extract gravitational waves from a simulation of rotating 
collapse to a BH in numerical relativity. QNM excitation is observed, as confirmed by 
more recent work [72; . 

• 1985 - Confirming the validity of Goebel's arguments j59j, Mashhoon [73] regards 
QNMs as waves orbiting around the unstable photon orbit and slowly leaking out, and 
estimates analytically some QNM frequencies in Kerr-Newman backgrounds. 

• 1985 - Schutz and Will develop a WKB approach to compute BH QNMs. 

• 1985 - Leaver [75l [76l [77] provides the most accurate method to date to compute 
BH QNMs using continued fraction representations of the relevant wavefunctions, and 
discusses their excitation using Green's function techniques. 

• 1986 - McClintock and Remillard [78] show that the X-ray nova A0620-00 contains 
a compact object of mass almost certainly larger than 3Mq, paving the way for the 
identification of many more stellar-mass BH candidates. 

• 1989 - Echeverria [9] estimates the accuracy with which one can estimate the 
mass and angular momentum of a BH from QNM observations. The formalism is 
later improved by Finn [79j and substantially refined in Ref. [lOj . where ringdown- 
based tests of the no-hair theorem of general relativity are shown to be possible. An 
Appendix of Ref. [TO] provides QNM tables to be used in data analysis and in the 
interpretation of numerical simulations; these data are now available online [47j . 

• 1992 - Nollert and Schmidt [8^ use Laplace transforms to compute QNMs. Froman 
et al. [ST] first introduce phase-integral techniques in the context of BH physics. 

• 1993 - Anninos et al. [82] first succeed in simulating the head-on collision of two 
BHs, and observe QNM ringing of the final BH. 

• 1993 - Bachelot and Motet-Bachelot [83] show that a potential with compact support 
does not cause power-law tails in the evolution of Cauchy data. Subsequently Ching 
et al. [531 m] generalize this result to potentials falling off faster than exponentially. 

• 1996 - Gleiser et al. [85] extend the perturbative formalism to second order and 
use it to estimate radiation from colliding BHs employing the so-called "close limit" 
approximation, quantifying the limits of validity of linear perturbation theory 86J . 

• 1997 - Maldacena i6] formulates the AdS/CFT (Anti-de Sitter/Conformal Field 
Theory) duality conjecture. Shortly afterward, the papers by Gubser, Klebanov, 
Polyakov ^87] and Witten [SB] establish a concrete quantitative recipe for the duality. 
The AdS/CFT era begins. 
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• 1998 - The AdS/CFT correspondence is generalized to non-conformal theories in a 
variety of approaches (see [35] for a review). The terms "gauge-string duality", "gauge- 
gravity duality" and "holography" appear, referring to these generalized settings. 

• 1998 - Flanagan and Hughes [53] show that, under reasonable assumptions and 
depending on the mass range, the signal-to-noise ratio for ringdown waves is potentially 
larger than the signal-to-noise ratio for inspiral waves in both Earth-based detectors 
(such as LIGO) and planned space-based detectors (such as LISA). 

• 1998 - Hod [17] uses earlier numerical results by NoUert [90J to conjecture that 
the real part of highly-damped QNMs is equal to Tin 3 (T being the Hawking 
temperature), a conjecture later proven by Motl [91] using the continued fraction 
method. Hod also proposes a connection between QNMs and Bekenstein's ideas on 
BH area quantization. 

• 1999 - Creighton [12] describes a search technique for ringdown waveforms in LIGO. 

• 1999 - Two reviews on QNMs appear: Quasinormal modes of stars and black holes, 
by Kokkotas and Schmidt [3] and Nollert's Quasinormal modes: the characteristic 
"sound" of black holes and neutron stars [5]. 

• 1999 - Horowitz and Hubeny [92] compute QNMs of BHs in AdS backgrounds of 
various dimensions and relate them to relaxation times in the dual CFTs. 

• 2000 - Shibata and Uryu [93] perform the first general relativistic simulation of the 
merger of two neutron stars. More recent simulations confirm that ringdown is excited 
when the merger leads to BH formation [46]. 




• 2001 - Birmingham, Sachs and Solodukhin [S^ point out that QNM frequencies of 
the (2 -I- l)-dimensional Bafiados-Teitelboim-Zanelli (BTZ) BH [95] coincide with the 
poles of the retarded correlation function in the dual (1 -I- l)-dimensional CFT. 

• 2002 - Baker, Campanelli and Lousto [SS] complete the "Lazarus" program to 
"resurrect" early, unstable numerical simulations of BH binaries and extend them 
beyond merger using BH perturbation theory. 

• 2002 - Dreyer [97j proposes to resolve an ambiguity in Loop Quantum Gravity using 
the highly damped QNMs studied by Hod [17] . 

• 2002 - Son and Starinets [92 formulate a recipe for computing real-time correlation 
functions in the gauge-gravity duality. They use the recipe to prove that, in the 
gauge-gravity duality, QNM spectra correspond to poles of the retarded correlation 
functions. 

• 2002 - QNMs of black branes are computed [99]. The lowest QNM frequencies of 
black branes in the appropriate conserved charges channels are naturally interpreted 
as hydrodynamic modes of the dual theory [100] . 

• 2003 - Motl and Neitzke [101] use a monodromy technique (similar to the phase 
integral approaches of Ref. [81]) to compute analytically highly damped BH QNMs. 

• 2003 - In a series of papers |1021 11031 1104j , Kodama and Ishibashi extend the Regge- 
Wheeler-Zerilli formalism to higher dimensions. 

• 2003 - In one of the rare works on probing quantum aspects of gravity with 
gauge theory in the context of the gauge-gravity duality (usually, the correspondence 
is used the other way around), Fidkowski et al. |105| study singularities of BHs 
by investigating the spacelike geodesies that join the boundaries of the Penrose 
diagram. The complexified geodesies' properties yield the large-mass QNM frequencies 
previously found for these BHs. This work is further advanced in Ref. |106| and 
subsequent publications. 

• 2004 - Following Motl and Neitzke [101] , Natario and Schiappa analytically compute 
and classify asymptotic QNM frequencies for d-dimensional BHs |107] . 
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• 2005 - The LIGO detector reaches design sensitivity |108j . 

• 2005 - Pretorius [l^ achieves the first long-term stable numerical evolution of a BH 
binary. Soon afterwards, other groups independently succeed in evolving merging BH 
binaries using different techniques [151 116j . The waveforms indicate that ringdown 
contributes a substantial amount to the radiated energy. 

• 2005 - Kovtun and Starinets [109j extend the QNM technique in the gauge-gravity 
duality to vector and gravitational perturbations using gauge-invariant variables for 
black brane fluctuations. A classification of the fluctuations corresponding to poles 
of the stress-energy tensor and current correlators in a dual theory in arbitrary 
dimension is given. These methods and their subsequent development and application 
in |110[ lllli I112[ I113j become a standard approach in computing transport properties 
of strongly coupled theories from dual gravity. 

• 2006-2008 - An analytic computation of the lowest QNM frequency in the shear 
mode gravitational channel of a generic black brane [114] reveals universality, related to 
the universality of the shear viscosity to entropy density ratio in dual gauge theories. 
This and further developments |115[ I116j also point to a significance of the QNM 
spectrum in the context of the BH membrane paradigm (for a recent review of the 
membrane paradigm approach, see |117j ). 

• 2008-2009 - QNM spectra are computed in applications of the gauge-gravity duality 
to condensed matter theory [371 [3H] • 

1.2. Notation and conventions 

Unless otherwise and explicitly stated, we use geometrized units where G = c = 1, so 
that energy and time have units of length. We also adopt the ( — h ++) convention 
for the metric. For reference, the following is a list of symbols that are used often 
throughout the text. 



d Total number of spacetime dimensions (we always consider one timelike 

and d — 1 spatial dimensions) . 
L Curvature radius of (A)dS spacetime, related to the negative 

cosmological constant A in the Einstein equations {G^i, + Ag^jy — 0) 

through £2 ^ ^[d - 2){d - I) / {2K). The - sign is for AdS, + for dS. 
M Mass of the BH spacetime. 

a Kerr rotation parameter: a = J /M £ [0, M]. 

Radius of the BH's event horizon in the chosen coordinates. 
u! Fourier transform variable. The time dependence of any field is ~ e~*"*. 

For stable spacetimes, Im(a;) < 0. Also useful is u) = lo/2ttT. 
LOji, LUj Real and imaginary part of the QNM frequencies. 
s Spin of the field. 

I Integer angular number, related to the eigenvalue Aim — l{l + d — 3) 

of scalar spherical harmonics in d dimensions, 
n Overtone number, an integer labeling the QNMs by increasing |Im(cj)|. 

We conventionally start counting from a "fundamental mode" with n = 0. 
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2. A black hole perturbation theory primer 

Within general relativity (and various extensions thereof involving higher-derivative 
gravity), QNMs naturally appear in the analysis of linear perturbations of fixed 
gravitational backgrounds. The perturbations obey linear second-order differential 
equations, whose symmetry properties are dictated by the symmetries of the 
background. In most cases, these symmetries allow one to separate variables with 
an appropriate choice of coordinates reducing the system to a set of linear ordinary 
differential equations (ODEs) or a single ODE. The ODEs are supplemented by 
boundary conditions, usually imposed at the BH's horizon and at spatial infinity. 
QNMs are the eigenmodes of this system of equations. The precise choice of 
the boundary conditions is physically motivated, but it is clear that the presence 
of the horizon, acting for classical fields as a one-sided membrane, is of crucial 
importance: it makes the boundary value problem non-hermitian and the associated 
eigenvalues complex. The methods used to reduce the problem to a single ODE 
depend on the metric under consideration; some of them are discussed and compared 
in Chandrasekhar's book [22]. Given the progress in the field in recent years and 
the vast literature on the subject, we will not attempt to describe these techniques in 
detail. As a simple example illustrating the main extensions of the formalism described 
in [22] we discuss field perturbations in d-dimensional, non- rotating geometries. For 
the interested reader. Sections O [H and [7] provide references on other background 
geometries. 



2.1. Perturbations of the Schwarzschild-anti-de Sitter geometry 

Consider the Einstein-Hilbert gravitational action for a dimensional spacetime with 
cosmological constant A: 

S = [ d'^Xy/^{R-2A)+ [ d'^Xy^Cm, (1) 



16nG „ 

where Cm is the Lagrangian representing a generic contribution of the "matter fields" 
(scalar, Maxwell, p— form, Dirac and so on) coupled to gravity. The specific form of 
Cm depends on the particular theory. The Einstein equations read 

Gf,. + Ag^, ^ %nGT^,, , (2) 

where T^^ is the stress-energy tensor associated with Cm- Eq. ^ should be 
supplemented by the equations of motion for the matter fields. Together with Eq. ([2]), 
they form a complicated system of non-linear partial differential equations describing 
the evolution of all fields including the metric. A particular solution of this system 
forms a set of background fields gj^^,^^''^, where $ is a cumulative notation for all 
matter fields present. By writing g^j. — g^^ + h^^, $ — + 0, and linearizing the 
full system of equations with respect to the perturbations /i^^ and </>, we obtain a set 
of linear differential equations satisfied by the perturbations. 

Maximally symmetric vacuum {T^J^ = 0) solutions to the field equations are 
Minkowski, de Sitter (dS) and anti-de Sitter (AdS) spacetimes, depending on the 
value of the cosmological constant (zero, positive or negative, respectively). Generic 
solutions of Eq. ^ are asymptotically flat, dS or AdS. We will be mostly interested 
in asymptotically flat or AdS spacetimes. AdS spacetimes of various dimension arise 
as a natural groundstate of supergravity theories and as the near-horizon geometry 
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of extremal BHs and p-branes in string theory, and therefore they play an important 
role in the AdS/CFT correspondence [35l ITTSl ITTOl ITM [T2T] . 

BHs in asymptotically AdS spacetimes form a class of solutions interesting 
from a theoretical point of view and central for the gauge-gravity duality at finite 
temperature. Their relation to dual field theories is discussed in Section[8l In addition 
to the simplest Schwarzschild-AdS (SAdS) BH, one finds BHs with toroidal, cylindrical 
or planar topology [1221 lESl HSll [1251 lESl [l27] as well as the Kerr-Newman- 
AdS family |128j . The standard BH perturbation theory [22] is easily extended to 
asymptotically AdS spacetimes |1021 11291 11301 1131j . For illustration we consider the 
non-rotating, uncharged d-dimensional SAdS (or SAdS^) BH with line element 

ds^ = -fdf + f-^dr^ +r^dnl_2, (3) 

where /(r) = 1 -f jh"^ — rp^'^/r'^^'^, d^^_^ is the metric of the (rf — 2)-sphere, 
and the AdS curvature radius squared I? is related to the cosmological constant by 
I? = —{d — 2){d — 1)/2A. The parameter Tq is proportional to the mass M of the 
spacetime: M ^ (d - 2)^d-2ro"Vl67r, where Ad-2 = 277^''-'^^^ /T[{d - l)/2]. The 
well-known Schwarzschild geometry corresponds to L oo. 



Scalar field perturbations Let us focus, for a start, on scalar perturbations in vacuum. 
The action for a complex scalar field with a conformal coupling is given by Sm = 
J d'^Xy/—g Cm, where 

/:m = -(a^$)^9^$-^^^7^*^*-"i'*^*- (4) 

For 7 = 1, TO = the action is invariant under the conformal transformations 
g^j^i, —> fl^gi^i, , $ 17^"'*/^$, and for 7 = 0, to = one recovers the usual minimally 
coupled massless scalar. The equations of motion satisfied by the fields g^i, and 
(massless) $ are 

= ^^^7i?$ , G,,, + A g^, = SttG T^, , (5) 

where T^^, is quadratic in $. Considering perturbations of the fields, g^^^, = g^^ + h^^j 
and $ = + (j) with = 0, we observe that the linearized equations of motion 
for hfj^n and (j) decouple, and thus the metric fluctuations h^n can be consistently set 
to zero. The background metric satisfies G^^ + Ag^^ = 0. We choose g^^ to be the 
SAdSd metric ([3]). The scalar fluctuation satisfies the equation 

The time-independence and the spherical symmetry of the metric imply the 
decomposition 

<^(t,r,0) = ^e-*H^y,„(0), (7) 

where Yim{d) denotes the d— dimensional scalar spherical harmonics, satisfying 
^Q.d-2^im = —lil + d — 3)yim, with Ao^_2 the Laplace-Beltrami operator, and the 
"s = 0" label indicates the spin of the field. Here and in the rest of this paper, 
for notational simplicity, we usually omit the integral over frequency in the Fourier 
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transform. Substituting the decomposition into Eq, 
for ^'s=o(r-): 



we get a radial wave equation 



+ {uj^ - V,=o) *.=o = . 



(8) 



dr^ '' '' dr 

We will see shortly that perturbations with other spins satisfy similar equations. In 
the particular case of s = 0, the radial potential Vs is given by 

'l{l + d~3) ^ d-2f{d-A)f 2/' 



^.=0 = / 



(9) 



Finally, if we define a "tortoise" coordinate r* by the relation dr^/dr = 1//, Eq. ^ 
can be written in the form of a Schrodinger equation with the potential Vs 

^^ + {.^-Vs)^s^O. (10) 

Notice that the tortoise coordinate —^ —oo at the horizon (i.e. as r ^ r_|_), but its 
behavior at infinity is strongly dependent on the cosmological constant: — > +c» for 
asymptotically-flat spacetimes, and r» —* constant for the SAdS^ geometry. 



Electromagnetic, gravitational and half-integer spin perturbations Equations for 
linearized Maxwell field perturbations in curved spacetimes can be obtained along the 
lines of the scalar field example above. To separate the angular dependence we now 
need vector spherical harmonics [^ 11321 [133] . In d = 4, electromagnetic perturbations 
can be completely characterized by the wave equation PH)) with the potential 

'1(1 + 1 



1 



/ 



(11) 



A comprehensive treatment of the four-dimensional case can be found in Ref. |132| 
for the Schwarzschild spacetime, and in Ref. [130J for the SAdS geometry. Higher- 
dimensional perturbations are discussed in Ref. jl34j . 

The classification of gravitational perturbations hf^,j{x) on a fixed background 
g^^{x) is more complicated. We focus on the SAdS4 geometry. After a decomposition 
in tensorial spherical harmonics, the perturbations fall into two distinct classes: odd 
(Regge- Wheeler or vector- type) and even (Zerilli or scalar-type), with parities equal 
to (-1)'+^ and (-1)', respectively [HI HIl EZl [T3n|. In the Regge- Wheeler gauge 
[H [5l [27l HHl 1136] . the perturbations are written as h^^^ — e~*'^*/i^i,, where for odd 
parity 



whereas for even parity 







ho{r) 

Ho{r)f 
Hi{r) 











hi{r) 








de 



YloiO); 



(12) 



Hi{r) 
H2{r)/f 





''K{r) 
r'^K{r) sin^ 



Yio{9). (13) 



The angular dependence of the perturbations is dictated by the structure of tensorial 
spherical harmonics [27l [21] [57l I135j . Inserting this decomposition into Einstein's 
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equations one gets ten coupled second-order differential equations that fully describe 
the perturbations: three equations for the odd radial variables, and seven for the 
even variables. The odd perturbations can be combined in a single Regge- Wheeler or 
vector- type gravitational variable ^'7=2' the even perturbations can likewise be 
combined in a single Zerilli or scalar- type gravitational wavefunction ■ The Regge- 
Wheeler and Zerilli functions (^'7=2 ^^'^ ^^=2' respectively) satisfy the Schrodinger- 
like equation (fTO|) with the potentials 



Vs=2 - fir 



1(1 + 1) 



6AI 



(14) 



and 



2j(^) 9M3 + SX^Mr^ + {1 + X) + 9M^ ( 



Xr + 



r-" [Ml + Xry 

The parameters /iq and hi of the vector-type perturbation are related to ^'J^2 by 



(15) 



*7-. = 



I{r) 



hi{r), 



i d 
ha = — -J— 

uj dr.. 



(16) 



For the scalar-type gravitational perturbation, 
expressed through ^'^^2 

6AP 



the functions Hi and K can be 



where X = (l 



K 
Hi 



+ A (1 + A) r2 + 3Af (Xr - f^) 



r2 (3Af + Xr) 
'iXMr - Xr^ + 3Af 



d^Z 



dr^ 



r(3Af + Ar)/(r) ^ /(r) dr, ' 

2)/2, and Ha is then obtained from the algebraic relation 



(17) 



(18) 



2) + 

r 



Ha 



.l{l + \) 



{l-l){l + 2)+rf - 



Auj^r^ + r^f 



Lo r^ 
2 fa 



2/ 



(M + r^ /L^)-2iujr Hi 



K = 0. (19) 



A complete discussion of Regge- Wheeler or vector-type gravitational perturbations 
of the four-dimensional Schwarzschild geometry can be found in the original papers 
by Regge and Wheeler [48] as well as in Ref. |137j . where some typos in the 
original work are corrected. For Zerilli or scalar-type gravitational perturbations, the 
fundamental reference is Zerilli's work |561 157] : typos are corrected in Appendix A of 
Ref. [138| . An elegant, gauge- invariant decomposition of gravitational perturbations of 
the Schwarzschild geometry is described by Moncrief [HI] (see also [551 1571 [T^ Il40j ). 
These papers are reviewed by Nollert |5j and Nagar and RezzoUa [23 • For an 
alternative treatment, see Chandrasekhar's book [55]. Chandrasekhar's book and 
papers [1411 I142j use a different notation, exploring mathematical aspects of the 
relations between different gravitational perturbations (see [Appendix A[ ). Extensions 
to the SAdS4 geometry can be found in Ref. [130] . while the general d— dimensional 
case has been explored in a series of papers by Kodama and Ishibashi }102[ 11031 1104| . 

The case of Dirac fields seems to have been discussed first by Brill and Wheeler 
[143] ■ with important extensions of the formalism by Page [144] . Unruh [145| and 
Chandrasekhar [52]. For the treatment of Rarita-Schwinger fields, see [146] . 
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To summarize this Section: in four-dimensional Schwarzschild or SAdS 
backgrounds, scalar (to = 0, 7 = 0, s = 0), electromagnetic (s — ±1) and Regge- 
Wheeler or vector-type gravitational (s = 2) perturbations, can be described by the 
master equation (fTU|) with the potential 



The potentials for the scalar-type gravitational perturbations and half-integer spin 
perturbations have forms different from ((20| . see for instance Refs. |1021 1147) . 
However, the vector-type (Regge- Wheeler) and scalar-type (Zerilli) potentials have the 
remarkable property of being isospectral, i.e. they possess the same QNM spectrum. 
The origin of this isospectrality, first discovered by Chandrasekhar \22^ , is reviewed in 




2.2. Higher- dimensional gravitational perturbations 

The literature on gravitational perturbations can be quite confusing. Naming 
conventions were already unclear in 1970, so much so that Zerilli decided to list 
equivalent terminologies referring to odd and even tensor spherical harmonics (cf. 
Table II of Ref. [S7j). The situation got even worse since then. Chandrasekhar's 
book, which is the most complete reference in the field, established a different 
terminology: "odd" (Regge- Wheeler) perturbations were called "axial" and described 
by a master variable while "even" (Zerilli) perturbations were renamed "polar" 
and described by a master variable VP^. In recent years, Kodama and Ishibashi 
|102[ 11031 I104j extended the gauge-invariant perturbation framework to higher- 
dimensional, non-rotating BHs. In higher dimensions three master variables are 
necessary to completely describe the perturbations [102j . Two of them (the vector- 
type gravitational perturbations and the scalar-type gravitational perturbations) reduce 
to the Regge- Wheeler and Zerilli master variables in d = 4. Kodama and Ishibashi 
refer to the third type of perturbations, which have no four-dimcnsional analogue, 
as tensor-type gravitational perturbations. In this review we will usually adopt the 
Kodama-Ishibashi terminology. 

2. 3. Weak fields in the Kerr background: the Teukolsky equation 

In four-dimensional asymptotically flat spacetimes, the most general vacuum BH 
solution of Einstein's equations is the Kerr metric. In the standard Boyer-Linquist 
coordinates, the metric depends on two parameters: the mass M and spin J — aAI . 
The spacetime has a Cauchy horizon at r = r_ — M—\J — and an event horizon 
at r = r-i- = M -|- \J ivP- — . The separation of variables for a minimally coupled 
scalar field in the Kerr background was first reported by Brill et al. |148j . 

Teukolsky [6OI I149j showed that if one works directly in terms of curvature 
invariants, the perturbation equations decouple and separate for all Petrov type-D 
spacetimes. He derived a master perturbation equation governing fields of general 
spin, including the most interesting gravitational perturbations (see [22l 1150] for 
reviews). Teukolsky's approach is based on the Newman-Penrose [IHI formalism. In 
this formalism one introduces a tetrad of null vectors 1 , n , m , m* at each point in 
spacetime, onto which all tensors are projected. The Newman-Penrose equations are 
relations linking the tetrad vectors, the spin coefficients, the Weyl tensor, the Ricci 
tensor and the scalar curvature [49j . The most relevant perturbation variables, which 




(20) 




[Appendix A[ 
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both vanish in the background spacetime, are the Weyl scalars 'i'o and 4^4, obtained 
by contracting the Weyl tensor C^^Xa |151j on the tetrad legs (roughly speaking, these 
quantities describe ingoing and outgoing gravitational radiation): 

*o = - Ci3i3 = -C^^xJ^m'^l^m'' , (21) 

*4 = - ^2424 = -C^.Aan''m*''n^m*'^ . (22) 
Two analogous quantities $o f^nd $2 describe electromagnetic perturbations: 

*o = F^J'^m'^ , *2 = F^,m* ^n" . (23) 

By Fourier-transforming a spin-s field ijj{t ,r ,9 ,(j)) and expanding it in spin- weighted 
spheroidal harmonics as follows: 



1 

2^ 



00 / 

-iujt \^ 

l—\s\ m—'-l 



i){t,r ,9 ,(, 

Teukolsky finds separated ODEs for gSim and Rim [60lll49 

sSlm 



'""^ sSim{9)Rim{r)du . (24) 



9 2^ 9 



2 2 2 

a uj u 



2aLUSU + S + sAln 



(m + suY 

1-m2 



AdlRim + {s + l)(2r - 2M)drRim + VRw, 



sSlm 

= 0. 



0, 



(25) 



Here u = cos 9, A = (r — r_)(r — r+) and 



V = 2isujr — a UJ 



2, ,2 



{r^ + a^Yuj^ - AManiLur + a m 



+ is {am{2r - 2M) - 2Muj{r'^ - a^)) 



(26) 



The solutions to the angular equation (j25p are known in the literature as spin- 
weighted spheroidal harmonics: sSim = sSimio.^ Foi' o,'-^ — the spin- 
weighted spheroidal harmonics reduce to spin- weighted spherical harmonics sYim{(^ j '/') 
|152j . In this case the angular separation constants s-^im are known analytically: 
sAimici = 0) = 1(1 + 1) — s{s + l). The determination of the angular separation constant 
in more general cases is a non-trivial problem (see |153j and references therein). 



Table 1. Teukolsky wavcfunction ip, as in l|24|l . for each value s of the spin. 
The spin-coefficient p = — l/(r — iacosB). The quantities xo and xi refer to the 
components of the neutrino wavcfunction along dyad legs. 







(+1/2,- 
{xa,P 



-1/2) 



(+1, 
(*o,P 



-1) 
-2*2) 



(+2, 
(*o,P" 



-2) 

4*4 



The field's spin weight s is equal to 0,±1/2,±1,±2 for scalar, Dirac, 
electromagnetic and gravitational perturbations, respectively. The Teukolsky master 
variable is related to the perturbation fields by the relations listed in Table [T] (see 
also Appendix B of Ref. |149j ). Relations between the Regge- Wheeler- Zerilli and 
the Teukolsky variables are explored in Ref. [22] . Reconstructing the metric from the 
Teukolsky functions is a highly non-trivial problem which is still not completely solved 
(see e.g. [TMllTCSllTCSlfWllTSglfTCg] ). 
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3. Defining quasinormal modes 

3.1. Quasinormal modes as an eigenvalue problem 

In a spherically symmetric background, the study of BH perturbations due to 
linearized fields of spin s can be reduced to the study of the differential equation 
(jlOp. Henceforth, to simplify the notation, we will usually drop the s-subscript in 
all quantities. To determine the free modes of oscillation of a BH, which correspond 
to "natural" solutions of this unforced ODE, we must impose physically appropriate 
boundary conditions at the horizon (r» — > — oo) and at spatial infinity {r^ oo). 
These boundary conditions are discussed below. 

Boundary conditions at the horizon For most spacetimes of interest the potential 
V ^ as — + — cxD, and in this limit solutions to the wave equation (|10|) behave 
as ^' ~ g"4i^(t±r,)^ Classically nothing should leave the horizon: only ingoing modes 
(corresponding to a plus sign) should be present, and therefore 

^ ^-i-^it+r.) ^ r^^-oo {r^ r+) . (27) 

This boundary condition at the horizon can also be seen to follow from regularity 
requirements. For non-extremal spacetimes, the tortoise coordinate tends to 

r.,= J f-Ur^if {r+)]-' log ir-r+), r ^ r+ , (28) 

with f'{r+) > 0. Near the horizon, outgoing modes behave as 

g-i(^(t-r,) _ ^-iujv^2iujr, ^ e~^'^'" {r — r+)'^^'^^^'^^'+^ , (29) 

where v = t + r^. Now Eq. ([29|) shows that unless 2iu/ f'{r+) is a positive integer the 
outgoing modes cannot be smooth, i.e. of class C°°, and they must be discarded. An 
elegant discussion of the correct boundary conditions at the horizon of rotating BHs 
can be found in Appendix B of Ref. [160j . 

Boundary conditions at spatial infinity: asymptotically flat spacetimes For asymp- 
totically flat spacetimes, the metric at spatial infinity tends to the Minkowski metric. 
From Eq. (|20p with L — s- cx) we see that the potential is zero at infinity. By requiring 

e-"^^'-'^'\ r^oo, (30) 

we discard unphysical waves "entering the spacetime from infinity" . 

The main difference between QNM problems and other prototypical physical 
problems involving small oscillations, such as the vibrating string, is that the system 
is now dissipative: waves can escape either to infinity or into the BH. For this reason 
an expansion in normal modes is not possible [H O [771 [HO] ■ There is a discrete infinity 
of QNMs, defined as eigenfunctions satisfying the above boundary conditions. The 
corresponding eigenfrequencies wqnm have both a real and an imaginary part, the 
latter giving the (inverse) damping time of the mode. One usually sorts the QNM 
frequencies by the magnitude of their imaginary part, and labels them by an integer n 
called the overtone number. The fundamental mode n = is the least damped mode, 
and being very long-lived it usually dominates the ringdown waveform. 

A seemingly pathological behavior occurs when one imposes the boundary 
conditions (P7)) and ((50)) . When the mode amplitude decays in time, the characteristic 
frequency wqnm must have a negative imaginary component. Then, the amplitude 
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near infinity (r* — > +00) must blow up. So it is in general impossible to represent 
regular initial data on the spacetime as a sum of QNMs. QNMs should be thought 
of as quasistationary states which cannot have existed for all times: they decay 
exponentially with time, and are excited only at a particular instant in time (see 
|161| for an alternative viewpoint on "dynamic" QNM excitation). In more formal 
terms, QNMs do not form a complete set of wavefunctions [5]. 



Boundary conditions at spatial infinity: asymptotically anti-de Sitter spacetimes 
When the cosmological constant does not vanish, by inspection of Eq. (fTO|) we see 
that 

'is=o-Ar-^ + Br, ^=1.2 - A/r + B , r ^ 00 . (31) 

Regular scalar field perturbations should have B — 0, corresponding to Dirichlet 
boundary conditions at infinity. The case for electromagnetic and gravitational 
perturbations is less clear: there is no a priori compelling reason for a specific boundary 
condition. A popular choice implements Dirichlet boundary conditions for the Regge- 
Wheeler and Zerilli variables [130j . but other boundary conditions were investigated, 
e.g., in Ref. |162j . A discussion of preferred boundary conditions in the context of the 
AdS/CFT correspondence can be found in Refs. [MIIIMIISS] (see also SectionH^]). 



3.2. Quasinormal modes as poles in the Green's function 

The QNM contribution to the BH response to a generic perturbation can be identified 
formally by considering the Green's function solution to an inhomogeneous wave 
equation [771 HOl I166|. I161|. I167j . Consider the Laplace transform of the field, 
/:*(t,r) = *(w,r) = ^{t,r)e''^'^dt, which is well defined if w/ > c (the usual 
Laplace variable is s —iuj; we use to for notational consistency with previous 
works). The problem of computing the gravitational waveform produced when a BH 
is perturbed by some material source (such as a particle of mass m <C M falling into 
the BH) can be reduced to a wave equation of the form (fTO|) with a source term: 

^ + {^^-V)^ = Iiu;,r). (32) 

We can solve this equation by the standard Green's function technique [168j (see 
Refs. [771 [501 [TB51 11611 1167j for applications in this context), focusing for definiteness 
on asymptotically flat spacetimes. Take two independent solutions of the homogeneous 
equation: one has the correct behavior at the horizon, 

lim , ~ e-^'^''- , (33) 
hm - A„(^)e-''"'-- + Aut(c^)e'""* , (34) 

r — >C30 

and the second independent solution '5 00+ ^ e*"*"* for large r. The Wronskian of these 
two wavefunctions is = 2iLuAin, and we can express the general solution as [77] 

The time- domain response is obtained by inversion of the Laplace transform: 

5f(i, r) = — *(w , r)e-*"*dw . (36) 

■^"^ J—oo+ic 
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Figure 2. Integration contour for Eq. II36I I. The liatclied area is the branch cut 
and crosses mark zeros of the Wronskian W (the QNM frequencies). 



The frequency integral can be performed by the integration contour shown in Fig. [21 
There are in general three different contributions to the integral. The integral along 
the large quarter-circles is the flat-space analogue of the prompt response, i.e. waves 
propagating directly from the source to the observer at the speed of light. Depending 
on the asymptotic structure of the potential, there is usually a branch point at w = 
|41| . To prevent it from lying inside the integration contour, we place a branch cut 
along the negative imaginary-w axis and split the half circle at \uj\ — > oo into two 
quarter circles [1691 11701 [51] . The branch-cut contribution gives rise to late-time 
tails |1711 [77l [84] : physically, these tails are due to backscattering off the background 
curvature [171] , and therefore they depend on the asymptotics of the spacetime. Tails 
are absent for certain backgrounds, such as the SAdS [9? or the Nariai spacetime [172] . 
The third contribution comes from a sum-over-residues at the poles in the complex 
frequency plane, which are the zeros of Ai^. These poles correspond to perturbations 
satisfying both in-going wave conditions at the horizon and out-going wave conditions 
at infinity, so (by the very definition of QNMs) they represent the QNM contribution 
to the response. Although we used the asymptotic behavior of the solutions of the 
wave equation, this discussion can be trivially generalized to any spacetime. 

Far from the source, the QNM contribution in asymptotically flat spacetimes can 
be written as (TMl [HZ! 



^i(t,r) = -Rc 



-iuJn {t—r^ ) 



(37) 



where the sum is over all poles in the complex plane. The C„'s are called quasinormal 
excitation coefficients and they quantify the QNM content of the waveform. They 
are related to initial-data independent quantities, called the quasinormal excitation 
factors (QNEFs) and denoted by _B„, as follows: 

/(c^,r)*,^ ^ / X -1 



Br, 



A-ont ( dA\^ 



2lu \ duj 



(38) 



In general the QNM frequencies ix>„, the -B„'s and the C„'s depend on ^, m and on the 
spin of the perturbing field s, but to simplify the notation we will omit this dependence 
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whenever there is no risk of confusion. 

The QNEFs play an important role in BH perturbation theory: they depend only 
on the background geometry, and (when supplemented by specific initial data) they 
allow the determination of the QNM content of a signal, i.e. of the C„'s. This has been 
known for over two decades, but relatively little effort has gone into understanding 
how these modes are excited by physically relevant perturbations. The QNEFs 
have long been known for scalar, electromagnetic and gravitational perturbations of 
Schwarzschild BHs [77l I173|. I166| I161j . and they have recently been computed for 
general perturbations of Kerr BHs [167j . In [167j it was also shown that for large 
overtone numbers {n —> oo) J5„ oc for all perturbing fields in a large class of 
non-rotating spacetimes. The excitation factors C„ have mainly been computed for 
the simple case where the initial data are Gaussian pulses of radiation |166|, I167j . The 
only work we are aware of studying the excitation factors C„ for a point particle falling 
into a Schwarzschild BH is Leaver's classic paper [77] . 

Besides the theoretical interest of quantifying QNM excitation by generic initial 
data in the framework of perturbation theory, excitation factors and excitation 
coefficients have useful applications in gravitational wave data analysis. First of all, a 
formal QNM expansion of the BH response can simplify the calculation of the self- force 
acting on small bodies orbiting around BHs. It was shown recently, using as a model 
problem the Nariai spacetime, that a QNM expansion of the Green's function can be 
used for a matched expansion of the "quasi-local" and "distant-past" contributions 
to the self force |172j . The problem of quantifying QNM excitation is of paramount 
importance to search for inspiralling compact binaries in gravitational wave detector 
data. All attempts to match an effective-one-body description of inspiralling binaries 
to numerical relativity simulations found that the inclusion of several overtones in the 
ringdown waveform is a crucial ingredient to improve agreement with the numerics 
|174[ 1175] . So far the matching of the inspiral and ringdown waveforms has been 
performed by ad hoc procedures. For example, the amplitudes and phases of the C„'s 
have been fixed by requiring continuity of the waveform on a grid of points, or "comb" 
|176] . These matching procedures have their own phenomenological interest, but a 
self-consistent estimation of the excitation coefficients within perturbation theory is 
needed to improve our physical understanding of the inspiral-ringdown transition. 

Finally, we point out that recent investigations have addressed the issue of mode 
excitation in the gauge/gravity duality |177| I178| . There it was found, for example, 
that the residue of the diffusion and shear mode decays at small wavelength, so these 
modes effectively cease to exist. 

4. Computing quasinormal modes 

To determine the QNMs and compute their frequencies we must solve the eigenvalue 
problem represented by the wave equation (jlOp , with boundary conditions specified by 
Eq. (|27p at the horizon and Eq. pop at infinity. There is no universal prescription to 
compute QNMs. In this section we discuss various methods to obtain such a solution, 
pointing out that different methods are better suited to different spacetimes. 

For a start we consider the exceptional cases where an exact, analytical solution 
to the wave equation can be found. In some spacetimes the potential appearing in the 
wave equation can be shown to reduce to the Poschl- Teller potential [179j , for which an 
exact QNM calculation is possible [70] . These spacetimes and their QNM spectra are 
reviewed in Section [4. II In the general case, QNM calculations require approximations 
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or numerical methods |j. Some of these (includmg WKB approximations, monodromy 
methods, series solutions in asymptotically AdS backgrounds and Leaver's continued 
fraction method) are reviewed in Sections 14.2114.61 



4-.1. Exact solutions 

In general exact solutions to the wave equation are hard to find, and they must be 
computed numerically. There are a few noteworthy exceptions, some of which we 
summarize here. 

We begin our review of exact solutions by sketching the analytical derivation of 
the QNMs of the Poschl- Teller potential. Many of the difficulties in computing QNMs 
in BH spacetimes arise from the slow decay of the potential as r oo, which is 
due (mathematically) to the presence of a branch cut and gives rise (physically) to 
backscattering of gravitational waves off the gravitational potential and to late-time 
tails. Ferrari and Mashhoon realized that these difhculties can be removed and exact 
solutions can be found if one considers instead a potential that decays exponentially 
as r — > oo, while recovering the other essential features of the Schwarzschild potential. 
Such a potential is the Poschl- Teller potential [7^. After reviewing QNM solutions 
for the Poschl- Teller potential we briefly review the modes of pure AdS and dS 
spacetimes. Then wc show that perturbations of the near-extreme Schwarzschild de- 
Sitter (SdS) and of the Nariai spacetime reduce to a wave equation with a Poschl- Teller 
potential, so they can be solved analytically. We also discuss two asymptotically 
AdS BH spacetimes which have improved our understanding of the role of QNMs 
in the AdS/CFT correspondence: the (2 + l)-dimensional BTZ BH [184] and the 
d-dimensional topological (massless) BH [122 [ml [ml [ESI [l26l [H^ . The BTZ BH 
is the first BH spacetime for which an exact, analytic expression for the QNMs has 
been derived ^95j, and it offers interesting insights into the validity of the AdS/CFT 
correspondence [Ml [113 [US] • 



The Pdschl-Teller potential In this section we analytically compute the QNMs of 
the Poschl- Teller potential [179| [70] , which will serve as a prototype for several BH 
spacetimes to be discussed below. Consider the equation 



dr'i 



Vn 



* = 0. 



(39) 



cosh a(r, — i 

The quantity r* is the point at which the potential attains a maximum, i.e. 
dV/dr^.{f^,) = 0, and Vq is the value of the Poschl- Teller potential at that point: 
Vq = ^(^*)- The quantity a is related to the second derivative of the potential at 
r* = f,, = ~{2VQ)''^(fV/drl{f^). The solutions of Eq. ^ that satisfy both 



boundary conditions (|27p and (|30p are the QNMs [70]. To find these solutions we 



define a new independent variable = [l 



+ e-2a(r.-r.)j rewrite Eq. Q as 

4a2 



C(l-0(2e-l)^ + 



9C 



k(i-e) 



= 0.(40) 



X Fiziev |180l I181| actually showed that the Regge-Wheeler and Teukolsky equations can be 
solved analytically in terms of confluent Heun's functions. These solutions allow a high-accuracy 
determination of the Schwarzschild quasinormal frequencies. The calculation of Heun's functions is 
quite contrived, but it already proved useful to address interesting open problems in perturbation 
theory [T82l[T83l . 
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Near spatial infinity 1 — ^ ^ e ^aC^, r-.)^ near the horizon ^ ^ g^^^r, r,) ^ j£ 
we define a = [a + \/q;2 - Wq ~ 2iuo\ /(2a) ,b = [a - ^Jd^ ~ 4Vb - 2iuj\ /(2a) , c = 
1 — I w/a and set 5* = (^(1 — ^)) '"^/(^"^ ■y^j^g get a standard hypergeometric equation 
for y [187] : 

+ [c - (a + 6 + l)e]9ey - afoy = , (41) 

and therefore 

* = Af" /(2a)^^ _ ^)-»c./(2a)^(^ - C + 1, 6 - C + 1, 2 - C, 

+ B (C(l-e))-*"/''")na,6,c,0. (42) 
Recalling that F(ai, a2, aa, 0) = 1 and that ^"^/(2") ^ e*'^'"' near the horizon, we 
see that the first term represents, according to our conventions, an outgoing wave 
at the horizon, while the second term represents an ingoing wave. QNM boundary 
conditions require A = 0. To investigate the behavior at infinity, one uses the z ^ 1 — z 
transformation law for the hypergeometric function I187j : 

F{aAc,z) = (l-z)--''£ffl^±^F(c-a,c-6,c-a-6+l,l-z) 

r(a)r(6) 

T(c)T(c — a — b) , , , , , 

+ ^ j^Fia,b-c+a+b+l,l-z). 43 

1 (c — a)i (c — b) 

The boundary condition at infinity implies either l/r(a) = or l/r(6) = 0, which are 
satisfied whenever 



Lu ^ ±^/Vo - a2/4 - ia{2n + l)/2 , n = 0, 1, 2, .... (44) 

where n is the overtone index. 

A popular approximation scheme to compute BH QNMs consists in replacing the 
true potential in a given spacetime by the Poschl- Teller potential. This approximation 
works well for the low-lying modes of the Schwarzschild geometry, ft predicts 
Muj = 0.1148 - 0.1148?, 0.3785 - 0.0905i for the fundamental ^ = s = 0, 2 
perturbations, respectively |70j . This can be compared to the numerical result 
[731 [ini EZ] Mw = 0.1105 - 0.1049ii, 0.3737 - 0.0890i. In the eikonal limit {I oo) 
the Poschl- Teller approximation yields the correct solution: it predicts the behavior 
SySMw = ±(/ + 1/2) - i{n + 1/2), in agreement with WKB-based calculations 
[SHUUHKHS] (see also [T90] 'l. 

The Poschl- Teller approximation provides a solution which is more and more 
accurate for near-extremal SdS BHs, since the event horizon and the cosmological 
horizon coalesce in the extremal limit |162[ I191[ 1192] . 



Normal modes of the anti-de Sitter spacetime A physically interesting analytical 
solution concerns the QNMs of pure AdS spacetime, which can be obtained by setting 
ro = in the metric ([3]). In this case, QNMs are really normal modes of the spacetime, 
and have been computed for scalar field perturbations by Burgess and Liitken [193j . 
They satisfy 

Luj = 2n + d^l-l, 71 = 0,1,2,... s = 0, (45) 

where I is related to the eigenvalue of spherical harmonics in d dimensions by 
-A-im = l{l + d—3) \153\. Normal modes of electromagnetic perturbations in rf = 4 were 
shown to be the same as normal modes for gravitational perturbations (see Appendix 
in Ref. [194j ): for general d, they can be computed from the potentials for wave 
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propagation derived by Kodama and Ishibashi [104 . Ref. [107j studies the normal 
modes of gravitational perturbations. Tensor gravitational perturbations obey the 
same equation as s = fields, and therefore their normal modes are 

Luj = 2n + d + l-l, n = 0, 1, 2,... s = 2 tensor - type . (46) 

Vector-type gravitational perturbations have normal modes with the following 
characteristic frequencies 

Luj = 2n + d + l-2, n = 0, 1, 2,... s = 2 vector - type . (47) 

Finally, scalar-type gravitational perturbations have a somewhat surprising behavior. 
For d = 4 they are given by Eq. (I47p . For d ^ 5 they have a continuous spectrum, 
and for d > 5 one finds [lOJI 

Luj = 2n + d + l-3, n = 0, 1, 2,... s = 2 scalar - type . (48) 



Quasinormal modes of the de Sitter spacetime The dS spacetime is an extensively 
studied solution of the Einstein field equations, most of the early investigations being 
motivated by cosmological considerations. It satisfies Eq. ^ with ro — and 
/(r) = 1 — r"^ /L"^ . Natario and Schiappa found that no QNM solutions are allowed 
in even-dimensional dS space [107]. For scalar fields and tensor- type gravitational 
perturbations in odd-dimensional dS backgrounds the QNM frequencies are purely 
imaginary, and given by 

Luj ^ ~i{2n + d + l) , 7i 0, 1, 2, ...tensor - type . (49) 

For the other types of gravitational perturbations Ref. |107] finds 

Llo ^ ~i{2n + d + l + l) , n = 0, 1, 2, ...vector - type , (50) 

Lu)=-i{2n + d + l + 2), ti = 0, 1, 2, ...scalar - type . (51) 

Notice the striking similarity with the pure AdS results when one replaces L — > iL. 
Fields of other spins were considered in Rcfs. |1951 1192] . 



Nearly- extreme SdS and the Nariai spacetime The metric for d-dimensional 
Schwarzschild de-Sitter (SdSd) BHs can be obtained from Eq. fS]) by the replacement 
L iL, i.e., /(r) = 1 — r"^ jL? — r^^^ /r'^~'^. The corresponding spacetime has two 
horizons: an event horizon at r = r+ and a cosmological horizon at r = Tc. It was 
observed in Ref. |196j that for rc/r^ — 1^1 perturbations of this spacetime satisfy a 
wave equation with a Poschl- Teller potential. In particular, setting fcb = (rc — r+)/2r^, 
for near-extreme SdSd one finds M/r^'^ - l/(c^-l), rl/L"^ - (d-3)/(d-l) [l9k[T97] . 
QNM frequencies for scalar field and tensor-like gravitational perturbations are then 
given by [Ml [197] 



^ = Y;(Z + d- 3) - i - i , tensor -type, (52) 

where k = {d — 3)(rc — r+)/{2r\) is the surface gravity of the BH. For gravitational 
perturbations the result is 



^ = J(/-l)(/ + d-2)-i - Jn+i) , vector -type, (53) 



— J {I — 1){1 + d — 2) — d + — i (n + -] , scalar — type . 
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Through an appropriate limiting procedure |196| I198j , the nearly-extreme SdSd 
geometry can yield a spacetime with a different topology, the Nariai spacetime 
[TMl nm EOOl EOl], of the form 

ds^ = - {-r^/L^ + 1) dt^ + {~r^/L^ + l)'^ dr^ + r^dnl_2 . (54) 

This manifold has topology dS2 x Sd-2 and two horizons (with the same surface 
gravity). The QNM frequencies are the same as for nearly-extreme SdS^, if k is 
replaced by the surface gravity of each horizon [198] . 

The BTZ black hole Ichinose and Satoh |202j were the first to realize that the wave 
equation in the (2 -I- l)-dimensional, asymptotically AdS BTZ BH [184) can be solved 
in terms of hypergeometric functions. An analytical solution for the QNMs of this BH 
was first found in Ref. 95J. The non-rotating BTZ BH metric is given by ^184 

ds^ = (r^/L^ - M) dt^ ~ (r^/L^ - M) dr^ - r^d(f^ , (55) 

where M is the BH mass and r_|_ = M^^^L is the horizon radius. In 2 -I- 1 dimensions 
gravity is "trivial" : the full curvature tensor is completely determined by the local 
matter distribution and the cosmological constant. In particular, in vacuum the 
curvature tensor R^^Xp — A {g^xg^p — g^ipgvx) and R — 6A. Curvature effects 
produced by matter do not propagate through the spacetime. There are no dynamical 
degrees of freedom, and no gravitational waves [2031 1204] , Therefore, we will focus 
on scalar fields and assume an angular dependence of the form e™'^. Scalar QNM 
frequencies are given by 

ujL = ±m-2iln+l)r+/L , (56) 

with TO the azimuthal number, and n the overtone number |95| . This result has been 
generalized by Birmingham, Sachs and Solodukhin [94] to the rotating BTZ BH. For 
general massive scalar perturbations with mass parameter /i they find 

ujL = ±TO-i[2n+ (1 + y/^FTl)] (r+ -r^)/L. (57) 

The BTZ background has provided a first quantitative test of the AdS/CFT 
correspondence: the QNM frequencies (|57)) match the poles of the retarded correlation 
function of the corresponding perturbations in the dual CFT [94] . Recently the QNMs 
of BTZ BHs were shown to be Breit-Wigner-type resonances generated by surface 
waves supported by the boundary at infinity, which acts as a photon sphere [205] . 
This interpretation is highly reminiscent of work in asymptotically flat spacetimes, 
interpreting QNMs as null particles slowly leaking out of circular null geodesies (see 
Refs. [5i]73l[70l]206llMl[2Qi[209] and Section below). 

An alternative to Einstein's gravity in three dimensions is the so-called 
"topologically massive gravity" , obtained by adding a Chern-Simons term to the action 
[21011211] . Topologically massive gravity allows for dynamics, i.e. gravitational waves. 
The QNMs of BTZ BHs in this theory have recently been computed, providing yet 
another confirmation of the AdS/CFT correspondence [185[ 1186] . 

Massless topological black holes It is also possible to obtain exact solutions in a 
restricted set of higher-dimensional BH spacetimes. These asymptotically AdS 
solutions are known as topological BHs. The horizon is an Einstein space of positive, 
zero, or negative curvature [1221 11231 11241 11251 11261 1127] . In the negative-curvature 
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case there is a massless BH playing a role quite similar to the BTZ BH in three 
dimensions. Consider the exterior region of the massless topological BH |212j 

ds^ = - {r^/L^ - 1) dt^ + {r^/L^ - l)"^ dr^ + r^da^ . 

This is a manifold of negative constant curvature with an event horizon at r = _L. 
Here da^ stands for the line element of a {d— 2)-dimensional surface T,d-2 of negative 
constant curvature. 

The wave equation for a massive scalar field with non-minimal coupling can be 
solved by the ansatz $ — \E'(r)e~''^*F, where F is a harmonic function of finite 
norm with eigenvalue — Q = — (d — 3) /4 — xhe parameter ^ is generically 
restricted, assuming only discrete values if S(j-2 is a closed surface. If the effective 
mass = — jd{d — 2)/ ALP' (with /i the mass of the field and 7 the conformal 
coupling factor) satisfies the bound m^^L^ > — [(d— 1)/2]2, one set of QNM frequencies 
is given by |212j 

ujL = ±i-i{2n+l + ^/{d- 1)2/4 + m2„X2^ , = 0, 1, 2... . (58) 

If the mass and the coupling constant 7 satisfy the relations \/{d— 1)2/4 + nil^^L^ = 

^t'" ~ 2 d-i-7(^d-2) ^ (^t'") ^ ^ mlffL"^ < 1 — (^^) ^, there is another set of modes 
for which the QNM frequencies are given by [212] 

uoL = ±S,-i{2n+l - {d - If / A + ml.L'^^ . (59) 

This computation represents the first exact analytic determination of QNMs in four 
and higher dimensions. The generalization to other fields (and in particular to 
gravitational perturbations) can be found in Ref . [213] . 



^.2. The WKB approximation 

Normal modes of vibration of an object usually have a simple interpretation in terms 
of waves traveling across or around the object. For example, the Earth's free modes of 
oscillation were highly excited and measured for the first time in the 1960 earthquake 
in Chile f214] . These (roughly) one-hour long periodic oscillations correspond to waves 
traveling around the globe, and carry information about the Earth's interior. Just like 
the Earth's free modes of oscillation, BH QNMs can be thought of as waves traveling 
around the BH [551 [SHI [7D1 [73 [Ml [Ml [Ml [Mj- More precisely, QNMs can be 
interpreted as waves trapped at the unstable circular null geodesic (also known as 
the light-ring) and slowly leaking out. The instability timescale of the geodesic is the 
decay timescale of the QNM, and the oscillation frequency lo ~ c/rLR, with c the 
speed of light and tlr the light-ring radius [ 208^ . 

This intuitive picture, first proposed by Goebel [53], is related to a more rigorous 
WKB approximation developed by Mashhoon [SS] and by Schutz and Will [74] (see 
also [111 [701 m [Ml [Ml [Ml IM])- Their derivation and results closely parallel the 
Bohr-Sommerfeld quantization rule from quantum mechanics. The procedure involves 
relating two WKB solutions across a "matching region" whose limits are the classical 
turning points, where oj'^ = V{r). The technique works best when the classical turning 
points are close, i.e. when lo^ ^ Vmax, where Vmax is the peak of the potential. Under 
these assumptions we can expand in a Taylor series around the extremum of the 
potential f*: 

Q = w'-V^~go + Qo('-*-^*)V2, Ql = d^Qldrl. (60) 
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In this region, the wave equation + Q"^ 



can be approximated by 
f = 0. 



(61) 



This equation has an exact solution in terms of parabolic cylinder functions [187[ I215j : 

^ = AD,{z) + BD_,_^{iz), z= (2Q(;)3e^f (r, -n), (62) 

with u = — «(5o/-\/2(5o — 1/2. Using the asymptotic behavior of cylinder functions 
|187[ I215j and demanding only outgoing waves at spatial infinity we get, near the 
horizon, 

,2 



* - Ae-'^'^z'^e-- - i^lTiA [r(-iy)]"' e^^/^z^^^^e^ . (63) 

QNM boundary conditions imply that the outgoing term, proportional to e^^/^, should 
be absent, so \ jV{—v) = 0, or z/ = n(= 0,1,2, ...). As we anticipated, the leading- 
order WKB approximation yields a "Bohr-Sommerfeld quantization rule" defining the 
QNM frequencies: 

Qo/V2Q[{ = i(" + l/2), n = 0,l,2,... (64) 




n (overtone number) n (overtone number) 

Figure 3. Percentage errors for the real (left) and imaginary part (right) of the 
QNM frequencies as predicted from WKB calculations. Thick lines: third-order 
WKB approximation; thin lines: sixth-order WKB approximation. 



Higher-order corrections to Eq. (pO)) have been computed [5H]. Iyer and Will 
|188U189] carried out a third-order WKB expansion, and more recently Konoplya [216] 
pushed the expansion up to sixth order. There is no rigorous proof of convergence, 
but the results do improve for higher WKB orders. This is shown in Fig. [Sj where we 
compare numerical results for the QNMs of Schwarzschild BHs from Leaver's continued 
fraction method (to be discussed in Section 15.11 below) against third- and sixth-order 
WKB predictions. The WKB approximation works best for low overtones, i.e. modes 
with a small imaginary part, and in the eikonal limit of large I (which corresponds to 
large quality factors, or large lor/oji). The method assumes that the potential has a 
single extremum, which is the case for most (but not all) BH potentials: see Ref. |103| 
for interesting counterexamples. 

Ref. [217] introduces a new, WKB-inspired asymptotic expansion of QNM 
frequencies and eigenfunctions in powers of the angular momentum parameter / -I- 1/2. 
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Their asymptotic expansion technique is easily iterated to high orders, and it seems 
to provide very accurate results in spherically-symmetric spacetimes. The asymptotic 
expansion also provides physical insight into the nature of QNMs, nicely connecting 
the geometrical understanding of QNMs as perturbations of unstable null geodesies 
with the singularity structure of the Green function. 

4-3. Monodromy technique for highly- damped modes 

A powerful variant of the WKB approximation, which is particularly useful in the 
highly-damped limit ui — s- oo, is the so-called monodromy technique [101] . The 
basic idea is related to Stokes' phenomenon in the theory of asymptotic expansions 
(see e.g. |218| for an excellent introduction to the topic). As shown by Andersson 
and Howls [219j . the monodromy technique is a simple variant of the phase-integral 
approach, whose application to BH physics dates back to the work by Froman et al. 



Let us consider the wave equation pop for the Schwarzschild geometry, but 
allowing r and to be complex variables. In the complex-r plane, solutions to 
Eq. (|10|) may be multivalued around the singular points (r = and r = 2M). To 
deal with the singular points we introduce branch cuts emanating from r = and 
r = 2M. The relation (r) is also multi- valued: in the Schwarzschild geometry 
r*(r) = r-|-2Af log {r/2M — 1), and we choose the branch such that log(— 1) = iir. We 
can now define a variable z = r^/2RI — itt which tends to zero as r — > 0. The Stokes 
lines are defined as the lines for which Re(r,) = '218], and they are shown in Fig. |4] 
(near the singular point r = 0, the Stokes lines form an angle 7r/4). The idea now is 

f \ / \ 



Figure 4. Left panel: contour for calculation of the QNM frequencies in the 
complex-r plane. The different regions are separated by the associated Stokes 
lines. Not shown in the plot are branch cuts from r = — oo to the origin and from 
r = 1 to point A. Right panel: integration contour in the complex-^ plane, with 
z = r, /2M — ITT. For more details see Ref. [101 . 



to equate the monodromy as computed in two different ways. The first computation 
takes the contour in Fig. |4l starting and ending at point A and following the Stokes 
lines, joined at infinity by the large semi-circles shown in Fig. [H If we start from A 
with a plane wave e*"^, we can extrapolate the behavior to near the singularity z — 0. 
This is because for large imaginary w the term is much larger than the other terms 




V 
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in the potential. Near the origin (r = 0) the solutions to the wave equation can be 
written as 



= A+\/47r Mwz J| {2Mloz) + A_ \/47r MujzJ^ | (2Mwz) , (65) 

where A± are constants and Js/2 is a Bessel function. From the asymptotics of Bessel 
functions for large loz one has 

*(z) - e'^'^''^'[A+e-'"++A^e-'°'-]+e''^^'''^'-[A+e'"++A^e'°'-] , (66) 

with a± — 7r(l ± s)/4. The second term must vanish because of the boundary 
conditions. Continuing along the contour, we make a 3/27r turn around r — 0. In 
terms of the z coordinate, this means a Stt rotation. Using the representation for the 
Bessel functions 

71—0 ^ ' ^ 

one can see that J±s/2(e^*''2Mw2) = e^'"'^'"''!'^ Jj^^i^^'^Mujz). Therefore 

^ _ g»2Afc.j:j^^g7»a+ + ^_ g^^- ] + g-^^^'^"^ [^+6^*"+ + A_ 6^*"- ] .(68) 

The e"*^*^'^^ term is exponentially small, since we are in the region where Re(2:) ^ 
and LOi <C 0. Using Eq. ([66|) we get for the coefRcient of e*'^^ the monodromy 

■ ^ - - . (69) 

Requiring that the second term in Eq. vanishes, Eq. ([55)1 yields — (1 + 2cos(7rs)) 
for the monodromy. On the other hand, the only singularity inside the contour is at 
r = 1, giving a factor e^***'^^^") and leading to the condition 

87rAfw = ±log3-i7r(2n+l) (70) 

for the highly damped QNM frequencies of Schwarzschild BHs with s = , 2 (s = 1 is 
predicted to have a vanishing real part). This formula is in agreement with numerical 
studies by NoUert [220j for the Schwarzschild spacetime. The procedure has been 
generalized to several other backgrounds, including Reissner-Nordstrom (RN) [lOlj . 
Kerr [HIllllll, SAdS [153] and other BH spacetimes [W], always finding excellent 
agreement with numerical calculations [2241 12251 1194] . 

The application of the monodromy method requires extreme care in identifying 
the appropriate Stokes lines and integration contours. Natario and Schiappa |107| 
point out some instances of inappropriate applications of the technique. A more 
rigorous mathematical treatment of these methods is highly desirable. 



Asymptotically anti-de Sitter black holes: a series solution 

With the exception of the monodromy technique, the methods discussed so far 
work for asymptotically flat or (with some minor modifications) dS spacetimes. 
For asymptotically AdS spacetimes, the perturbation equations exhibit regular 
singularities both at the horizon and at spatial infinity. Local solutions in the vicinity 
of the regular singular points are represented by convergent Frobenius series [226j . 
In many cases, the radius of convergence of the series is equal to or larger than the 
interval of interest. In this simple numerical procedure (implemented e.g. by 

Horowitz and Hubcny in 92]) is possible. Consider, for example, a scalar field 5" 
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in the SAdS4 spacetime, satisfying the wave equation (|TO|l . Defining $ for a generic 
wavefunction as $ = e"^''*^', we find 

/«^ + [/'-2^-]^-^<i>^0, (71) 

where /' = df /dr. To find the frequencies u that satisfy the boundary conditions we 
first note that Eq. ((7T|l has only regular singularities in the range of interest. To deal 
with the point at infinity, we first change the independent variable r to x = l/r and 
get 

(a;-x+)s(x)— +tx — + ^^$ = 0, 72 

dx^ dx X — x+ 

where a;+ = l/r+ = h, s = —fx'^/{x — x^), t = — [2/x''^ — f'x^ — 2iujx'^'^, u = 
V{x — x^)/ f. Eq. ([7^ admits a local (near x — h) Frobenius solution of the form 

oc 

$(a;) = (.T-/i)"^a„(x-/i)", (73) 

n=0 

where a„(ci;) is a function of the frequency. The radius of convergence of the series is 
limited by the distance to the next nearest singular point of the equation. The index 
a is determined by imposing the boundary condition at the horizon. Writing 

oc oo oo 

s{x) = ^ Si{x-hY , t{x) = ^ U{x-hy , u{x) = ^ Ui{x-hy , (74) 

1=0 j=0 i=0 

and substituting Eqs. ([75)1 and ([7^ in Eq. ([7T|) we get an indicial equation for a: 
a{a — l)so + Oito + mq — 0. We also have sq = 2h^K , = 2/i^(/c — iu) and uq — 0, 
where the surface gravity is k = /'(r+)/2. Therefore, either a = or a = iu/K. One 
can check that the second option corresponds to outgoing waves at the horizon, so we 
choose a = 0. Inserting the decomposition (|73p with a = in Eq. (j7ip and comparing 
powers of {x — h) we get a recurrence relation for the a„'s: 

^ n — l 

a„ = -, — ■ — — [iii - l)s„-i + itn~t + Un-i] at . (75) 

n{n - l)so + nto 

Since the differential equation is linear, oq is just an arbitrary normalization constant. 
Using the boundary condition $ = at infinity {x = 0) we finally obtain a simple 
relation determining the QNM frequencies: 

oc 

^a„(-/i)" = 0. (76) 

Tl = 

In practice one computes the A^-th partial sum of the series ([75)1 and finds the roots lu 
of the resulting polynomial expression, checking for convergence by comparison with 
(say) the roots obtained from the N + 1-th partial sum. The method can be easily 
implemented: see e.g. the publicly available Mathematica notebook [47j that was used 
in [130j to compute low-lying modes of the SAdS4 geometry. This notebook can be 
easily modified to deal with other geometries. 

Some mathematical aspects of the series solution method as applied to QNMs in 
asymptotically AdS backgrounds are discussed in Ref. [59] • 
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4-5. Asymptotically anti-de Sitter black holes: the resonance method 

The series solution described in the previous subsection converges very slowly for small 
SAdS BHs. Fortunately, in this regime there is a simple alternative: the resonance 
method [227]. This method requires a numerical integration of the wave equation, 
but unlike the original numerical procedure by Chandrasekhar and Detweiler ^62] , one 
only needs to search for roots on the real-w line. 

It is well known that quasi-bound states manifest themselves as poles in the 
scattering matrix, and as Breit-Wigner resonances in the scattering amplitude. 
Chandrasekhar and Ferrari made use of the form of the scattering cross section near 
these resonances in their study of gravitational wave scattering by ultra-compact stars 
|2281 1229j . For quasi-bound, trapped modes of ultra-compact stars, the asymptotic 
wave amplitude at spatial infinity ^ a cos Lor + (3 sin cor has a Breit-Wigner-type 
behavior close to the resonance: 

-f « constant [{uj - lorY + Loj] , (77) 

where loJ^ is the lifetime of the quasi-bound state and its characteristic "energy". 
The example of ultra-compact stars shows that the search for weakly damped QNMs 
corresponding to quasi-bound states {oj = lor — i^i with oj/ <C (^r) is extremely 
simplified. One locates the resonances by looking for minima of -\- 0^ on the 
u> = ujR line, and the corresponding damping time wi can then be obtained by a fit to 
a parabola around the minimum [2281 1229] . 

It was shown in [227j that these ideas can be used very successfully in BH 
spacetimes, if one integrates in from spatial infinity to the BH horizon (i.e., the 
appropriately re-defined quantities -\- are taken close to the horizon). The 
resonance method is extremely simple and accurate, especially for small SAdS BHs 
|227| (see also Section I6.ip . It is presently unclear whether it can be applied 
successfully to study QNMs in other BH spacetimes. 

Jf..6. The continued fraction method 

Applications of the continued fractions to eigenvalue problems have an interesting 
history dating back to Jaffe's 1933 paper on the spectrum of the hydrogen molecule 
ion, or perhaps to even earlier times (some details and relevant references can be found 
in [99]). 

The method was applied to gravitational problems by Leaver, leading to what 
is possibly the most successful algorithm to date to compute QNM frequencies 
[751 [75] [77] . Leaver's approach is based on the observation that the Teukolsky equation 
is a special case of a general class of spheroidal wave equations that appear in the 
calculation of the electronic spectra of the hydrogen molecule ion [76] . These equations 
can be solved through a three-term recurrence relation, and the boundary conditions 
lead to a continued fraction equation characterizing the QNMs. Originally devised 
for the Schwarzschild and Kerr geometries |2201 12301 12311 1225 1 , the method has also 
been applied to RN BHs [232], Kerr-Newman BHs [233], higher-dimensional BHs 
[2341 [2M1 (2361 [237] . SdS BHs [2381 ESS] and acoustic BHs [240] (see Section MM 
below), among others. We illustrate the main ideas by considering the Kerr case. Start 
with the following series solution for the angular eigenfunctions defined in Eq. (|25p : 

oc 

sSUu) = e'^"" (1 + uf- (1 - uf^ + ' (^8) 

p=0 
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where k± = |to±s|/2. The expansion coefficients Up are obtained from the three-term 
recurrence relation 

aoai + /3oao = , apUp+i + PpUp + ■ypUp-i = , p=l,2... (79) 

where the constants a^, Pp, jp are given in the original work [75 . Given a complex 
argument uj, the separation constant s^im can be obtained solving numerically the 
continued fraction f3o — ^^^^^z---- = or any of its inversions [75] , 

A solution Rr^ of the radial equation (|25|) should satisfy the appropriate boundary 
conditions: 

lim Rr. - (r+ - r^)-'^-'+"^+"'+e"^'^+{r - r+)-'-"'+ , (80) 

lim Rr^ - ^L(w)r-l-2s+zc^gzc^r (g^) 

In these relations (j+ — (a;r+ — am) /6, and b = \/l — 4a^. A convenient series solution 
close to the horizon can be found by methods due to Jaffe (see [75]): 

oo y- \ n 

R^^ = e"^''{r-r_)-^-'+'"^+"'+ {r-r+)-'-"'^ ("—^j .(82) 

n=0 ' 

The coefficients aj'j can be obtained from a recurrence relation similar to Eq. (j79p . 
The continued fraction method is very powerful at computing overtones because the 
n-th overtone turns out to be the most stable numerical root of the n-th inversion of 
the radial continued fraction (which in principle should be completely equivalent to 
any other inversion) [75j . This observation makes it (relatively) easy to numerically 
compute Kerr overtones up to n 50, and Schwarzschild overtones up to n '--^ 100; 
refinements of the technique to compute even higher overtones will be discussed 
in Section 15.31 below. Given the QNM frequencies and the corresponding angular 
eigenvalues sMm it is a simple matter to compute i?^^ for moderate values of r 
(the convergence of the expansion gets worse for large values of r). We provide an 
implementation of the method to compute the QNM frequencies, as well as the radial 
and angular wavefunctions, in a publicly available Mathematica notebook |47| . On 
the same web page we also provide numerical data for the QNM frequencies of the 
first 8 Kerr gravitational modes (n — 0, . . . , 7) with 2 < Z < 7, as computed by an 
independent Fortran code which will be made available upon request [241j . 

5. The spectrum of asymptotically flat black holes 

In this section we briefly review the structure of the QNM spectrum of BHs belonging 
to the Kerr-Newman family. More details can be found in Ref. |241j . We focus our 
discussion on gravitational perturbations, and refer to relevant works on other kinds 
of perturbations in the appropriate sub-sections. 

5. 1 . Schwarzs child 

For Schwarzschild BHs, scalar-type and vector-type gravitational perturbations can 
be shown to give rise to the same QNM spectrum. This remarkable result, due to 
Chandrasekhar and collaborators [S^ [141| . is reviewed in [Appendix A| below. 

The vector-type potential is simpler than the scalar-type potential, so it is 
customary to compute QNMs for vector-type perturbations. We computed QNM 
frequencies using the continued fraction technique (see Section 14. 6p as improved by 
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Figure 5. Top: QNM frequencies for gravitational perturbations with I = 2 
(black circles) and I = 3 (red diamonds). In both cases we mark by an arrow 
the algebraically special mode, given analytically by Eq. 1 183 I I; a more extensive 
discussion of this mode is given in |Appendix A| Notice that as the imaginary part 
of the frequency tends to infinity the real part tends to a finite, I -independent limit. 
Bottom: comparison of the I = \s\ QNM frequencies for scalar, electromagnetic 
and gravitational perturbations. 
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Nollert [220]. Schwarzschild QNM frequencies with I = 2 and I ~ 3 are shown in 
Fig.m The data are also available online [47] , 

It is apparent from Fig. [5] that the gravitational QNMs of Schwarzschild BHs 
are naturally divided in two categories. A mode whose frequency is (almost) purely 
imaginary separates the lower QNM branch from the upper branch. This mode is very 
close to (and may actually coincide with) the so-called "algebraically special mode" 
[142] . discussed in [Appendix A| It is located at 



and it can be taken as roughly marking the onset of the asymptotic high-damping 
regime. The algebraically special mode corresponds to an overtone index n = 9 when 
Z = 2, to n = 41 when I = 3 and to even larger values of rt as Z increases. This 
means that the asymptotic high-damping regime is very hard to probe numerically as 
we approach the eikonal limit {I ^ 1). The existence of an almost purely imaginary 
QNM frequency is unique to gravitational perturbations: for other known fields, the 
approach to the asymptotic regime is monotonic (see the bottom panel of Fig. [5] and 
Ref. [234j ). For gravitational perturbations we will, somewhat arbitrarily, define the 
weakly damped (highly damped) regime as corresponding to imaginary parts smaller 
(larger) than the algebraically special frequency. For spins of other fields there is no 
such clear marker (see Fig. [5]), but we will usually call "weakly damped" modes those 
with n < 10, and "highly damped" modes those with n ^ 10. 

The weakly damped modes The weakly damped gravitational modes were computed 
numerically by Chandrasekhar and Detweiler [G^, Leaver [75] and many others, 
and they are available online |47j. Iyer [189] computed the first few modes for 
scalar, electromagnetic and gravitational perturbations using the WKB technique, and 
compared results against the more accurate continued fraction method. In geometrical 
units, the fundamental I — \s\ mode is AIui — 0.1105 — 0.1049* for s = and 
Mu! = 0.2483 — 0.0925i for s — 1. For ringdown detection from astrophysical BHs the 
most relevant QNM is, in most situations, the fundamental gravitational mode with 
I — 2, with Mu! — 0.3737 — 0.0890i. This mode has oscillation frequency and damping 
time given by 



where M ~ lOM© is a typical value for a stellar-mass BH (see Section l973l below) . 

Weakly damped QNMs for Schwarzschild perturbations with half-integer spins 
have been computed in the WKB approximation by Cho and collaborators |1471 1242j . 

The highly damped modes Using a variant of Leaver's method, Nollert carried out the 
first reliable numerical calculation of highly damped QNM frequencies for gravitational 
perturbations |220| . The real part of the QNM frequencies is well fitted, for large n, 
by a relation of the form ujn — ujoo + \s,i/\/n. The leading-order fitting coefficient 
2Mti;oo — 0.08742 is independent of I, and it has the same value for s = 2 and 
5 = |224j . Nollert 's numerical results have been confirmed by various analytical 



Mloi k ±i{l - l)l{l + l){l + 2)/12 , 
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calculations. Motl pjl] used Leaver's continued fraction conditions to show that in the 
limit 71 — > CX3 the following asymptotic expansion holds: 

tj-rin3-(2n+l)7riT + 0(71^1/2), (86) 

where T = (87rAf)~^ is the Hawking temperature. This conclusion was confirmed by 
complex-integration [lOlj and phase-integral techniques pi9j . 

Corrections of order ^ n~^^^ to Eq. ([55)) were first obtained by Neitzke |243| 
and Maassen van den Brink [244j . Ref. [245] developed a systematic perturbative 
approach to determine subdominant corrections as n —^ oo, using a solution of the 
Regge- Wheeler equation in terms of Bessel functions. Ignoring contributions of order 
0{l/n), the result for s = and s = 2 can be written as 

uj , , , 1 + i 3Z(/ + 1) - (s^ - 1) , 
- = - 2n+l 7ri + ln3 + ^^^= ^ ' V — ^ \ - \ - 87 

The subdominant coefficients are in agreement with fits of numerical results [2201 12241 
I234j . For electromagnetic perturbations, analytic and numerical results suggest that 
the real part lor^ — > in the asymptotic limit [101[ I234j , which is also apparent from 
Fig-O The above calculation has been generalized to massless fermionic perturbations. 
By solving the Teukolsky equation in terms of confluent hypergeometric functions, 
Refs. [2461 1247| confirmed the known results for integer spins, and found in addition 
that the result for s = 1/2 and s = 5/2 can be written as: 

For s — 3/2 there is no correction at order 0{n^^^'^). Numerical studies on s = 1/2 
fields confirm that lur ^ and the spacing Suji 2'kT = 1/4M as tt, ^ oo [248) . 

In summary, numerical and analytical results for Schwarzschild BHs are in perfect 
agreement. As \luj\ — + oo, cijj^ — )■ Tin 3 for scalar and gravitational oscillation 
frequencies, and ujji —^ for perturbations of other spins. The spacing of the imaginary 
parts for large n is always given by 2ttT. By considering the scattering amplitude in 
the Born approximation, Padmanabhan showed that this universality in the spacing 
arises from the exponential redshift of the wave modes close to the horizon [249] . 

The eikonal limit By their own nature, WKB methods become increasingly accurate 
for large I, and they can be used to compute the QNM frequencies analytically when 
Z > 1. Up to order 0{l-^) the result is [58l [T88l [T89l [T90] 
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The convergence of the series gets worse with increasing n (cf. also Fig. [3|). In the 
eikonal limit the asymptotic behavior of the potential is not important, so it should 
not be surprising that a Poschl- Teller approximation of the Schwarzschild potential, 
discussed in Section yields the correct result (to leading order) for large I. The 
eikonal regime of Schwarzschild BHs is related to the properties of unstable circular 
null geodesies |58| [59l [TOl HSl [Ml [Ml Iffl [209] : the leading order of Eq. [89] can be 
written as uj = ^Ic I — i {n + 1/2) A, where fie is the orbital frequency and 1/A is the 
instability timescale of the unstable circular null geodesic. Such a connection can be 
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generalized to any asymptotically fiat, static spacetime in four and higher dimensions 
[208] , 

5. 2. Reissner-Nordstrom 

With a few exceptions, BH charge is usually considered astrophysically negligible. 
Despite this, the RN metric is of more than academic interest: for example, charged 
naked singularities have been proposed as classical models for elementary particles (see 
[233 and references therein). Handling scalar fields in the background of a charged 
BH requires only a straightforward generalization of the uncharged case, resulting in 
a wave equation of the form (jlOp . Electromagnetic and gravitational perturbations 
are more technically challenging, since they are coupled through the Einstein-Maxwell 
equations. It is still possible to reduce the problem to the study of two wave equations 
of the general form (llOp for two wavefunctions , f22]: 

_ [t' - 2Mr + Q^) 



u/)zt^^V<-^zt\ (90) 



Here gi +53 = 6M, -qiq2 = 'iQ^{l-l){l + 2) and i,j = 1,2 (i ^ j). In the limit when 
the charge Q/M goes to zero, Z^ , Z2 describe pure electromagnetic and pure vector- 
type gravitational perturbations of a Schwarzschild BH, respectively. In the limit of 
maximally charged BHs (Q/M = 1) the wave equations have a different singularity 
structure, and deserve a special treatment [250j . 



0.20 




0.17 



0.50 



0.60 



2Mco, 



0.70 



0.80 



R 



Figure 6. Trajectory the complex-frequency plane described by the fundamental 
QNM as the charge Q/M is increased. The solid line corresponds to i = 2 and 
Z2 (which reduces to vector- type perturbations as Q/M — > 0); the dashed line 
corresponds to / = 1 and (purely electromagnetic perturbations in the limit 
Q/M — > 0). The two modes coalesce in the extremal limit {Q/M — > 1). 



The weakly damped modes The behavior of the weakly damped modes of 
gravitational/electromagnetic perturbations is exemplified in Fig. [B] The solid line is 



Quasinormal modes of black holes and black branes 



36 



the trajectory described in the complex-frequency plane by the fundamental QNM 
with I — 2 corresponding to Zj" (perturbations that reduce to the vector-type 
gravitational Schwarzschild case as Q/M — 0). The dashed line is the trajectory 
of the fundamental QNM with I — 1 corresponding to (which reduces to purely 
electromagnetic perturbations as Q/M — > 0). The Schwarzschild limit corresponds 
to the bottom left of each curve, and the trajectories are described clockwise as 
Q increases. The real part of the frequency grows monotonically with Q, and the 
imaginary part shows an extremum. A striking feature is that modes of Z2 with 
angular index / coalesce with modes of with index (l — 1) in the extremal limit. 
We will elaborate on the significance of this result later in this section. In general 
there are two algebraically special frequencies for RN BHs [142j . located at 

(/-1)Z(? + !)(/ + 2) 

^^1 2 = ±7: , (91) 

2 3M T V^AP + 16Q2(/ -!)(/ + 2) ^ ' 

Using the numerical procedure described by Chandrasekhar and Detweiler |62| , Gunter 
computed the lowest- lying QNMs of electromagnetic/gravitational perturbations [251j . 
Extensive comparisons of Gunter's results against WKB predictions were done by 
Kokkotas and Schutz |252| . The continued- fraction method can be generalized to 
charged BHs with relatively minor modifications |232| . and it was used to compute 
numerically the weakly and highly damped QNMs of gravitational/electromagnetic 
perturbations ISMl [MS ^M- Weakly damped modes of massive 

charged fields were computed by WKB techniques in Ref. [255J. Dirac QNMs in 
the Kerr-Newman metric (which includes RN as a special case) have been computed 
by continued fractions in Ref. [256j . For generalizations to higher-dimensional charged 
BHs, see j257j and references therein. 

The highly damped modes Let us focus our attention on the high-damping regime for 
Z2] results for Z^ are similar. For a more detailed discussion of highly-damped modes, 
see [2241 1241] . The modes' behavior is better understood by looking separately at the 
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Figure 7. Numerical and analytical predictions for the real and imaginary part 
of the RN QNM frequencies as a function of charge for n = 5 X 10'^. 



real and imaginary parts of their frequencies as functions of charge. Numerical results 
are shown in Fig. [71 Frequencies and damping times oscillate as a function of charge 
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whenever Q is larger than some Qcrit, and Qmt/M decreases with overtone number. 
For moderate-to-large values of Q/M , numerical data are in excellent agreement with 
the analytical prediction by Motl and Neitzke [1011 I219|, I107j , valid for scalar and 
electromagnetic- gravitational perturbations of a RN BH in the large-damping limit: 

e^/T + 2 + 3e«3' ^) = , T = ^ - QV [2^r^] , (92) 

where r+ = M + yM^ — is the radius and T is the Hawking temperature of the 
outer horizon. The complex solutions of Eq. ((5^ exactly overlap with the oscillations 
observed in Fig. [7] for large enough values of Q/M. This agreement gives support to 
the asymptotic formula (j92p . while confirming that the numerics are still accurate for 
large values of n and Q/M |243j . However, Eq. ([5^ presents us with some striking 
puzzles. The predicted asymptotic RN QNM frequencies do not reduce to the expected 
Schwarzschild limit of Eq. dB!]): one finds instead lur ^ Tlnb as Q/M ^ PUTIIMH] . 
Eq. (|92p should also be taken with care in the extremal limit {Q/M — > 1): in this 
case the inner and outer horizons coalesce, the topology of the singularities in the 
differential equation changes, and the problem requires a separate analysis 107 . Such 
an analysis shows that the asymptotic oscillation frequency for extremal RN BHs is 
not given by the limit of Eq. ([^^ as Q/Ad 1. Instead, the modes are purely damped 
{uiR 0). An interesting classification of the solutions of Eq. ((5^ can be found in 
Fig. 3 and Section 4.4 of [219] . Besides "spiralling" solutions, the equation also admits 
periodic solutions when y'l — Q^ /M"^ is a rational number, and even pure imaginary 
solutions that may not be QNMs at all. A survey of highly-damped QNM spectra of 
charged BHs can be found in Ref. |107j . For a complete account of the asymptotic 
QNM spectrum of several perturbing fields in the RN metric, see Ref. |258) . 

The eikonal limit The eikonal regime (/ 3> 1) is well described by a WKB analysis. 
The lowest-order WKB approximation yields 

ujR^ {l + l/2)n,, Loi ^ -1/2 17, y^3A//ro ~ ^Q'^/rl , (93) 

where and f2c are the radius and frequency of the unstable circular null geodesic, 
Tc = (3M + ^9M2 - 8Q2^ /2 and VL^ = ^jM/rl ~ Q'^/r^ [2521 1254]. QNMs in this 
regime can also be interpreted, from a geometrical-optics point of view, as waves 
trapped at the unstable circular null geodesies [73l I208j . 

Extremal Reissner- Nordstrom A direct application of Leaver's method fails in the 
extremal limit. In this limit the two horizons coalesce, and the radial wave equation 
has irregular singular points at the horizon and at infinity. Onozawa et al. |250j 
managed to reduce the problem to a five-term recurrence relation. As anticipated from 
Fig. |6l the QNM spectrum for extremal RN BHs is characterized by an isospectrality 
between electromagnetic and gravitational perturbations: modes of Z2 with angular 
index I coalesce with modes of with index (Z — 1) in the extremal limit. This has 
been shown to be a manifestation of supersymmetry [2591 12601 1261) . 

The resulting QNM spectrum [2411 12241 1250j looks qualitatively similar to the 
Schwarzschild spectrum of Fig. (5] The real part first decreases, approaching the 
pure- imaginary axis as the overtone index grows. A QNM seems to be located at 
w — (0, —3.0479), while Chandrasekhar's formula (|9T|) predicts a mode dX — (0, —3). 
The numerical method used so far becomes unreliable for |a;/| > 10, and better 
techniques will be needed to verify analytical predictions [1011 1107] . 
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5.3. Kerr 

The Kerr QNM spectrum has a rich and complex structure [75l [2301 ESU [2411 [262] , 
The most relevant feature of the spectrum is that rotation acts very much like an 
external magnetic field on the energy levels of an atom, causing a Zeeman splitting 
of QNM frequencies. The determination of the QNM frequencies is tangled to the 
solution of an angular equation, the spin-weighted spheroidal harmonic equation (see 
Sections 12.31 and 14. 6p . In the general case aw ^ there are no known closed- form 
solutions for the separation constant Aim or for the spheroidal harmonics. However 
the spheroidal harmonics satisfy various symmetry properties |1531 175| , namely: 

(i) Eigenvalues for negative and positive m are related: s-Aijn — S'Aj_ . 

(ii) Eigenvalues for negative and positive s are related: -sAim = sAim + 2s. 

(iii) If w and -sAim correspond to a solution for given (s, I, m), another solution can 
be obtained by the following replacements: m — m, uj —t^*, -sAim -s^;*-m- 
These symmetries are usually a source of some confusion, so we give an explicit 
numerical example. For a — 0.6M, the methods described in Section 14.61 yield the 
following eigenvalues for s = —2: 

Mlo22 = 0.49404 - 0.08376i , Muj2-2 = -0.49404 - 0.08376i , (94) 
Mlo22 = - 0.31678 - 0.08889i , Muj2-2 = 0.31678 - 0.08889i . (95) 




Figure 8. Left: "Zeeman-like" splitting of the fundamental gravitational 
mode with 1 = 2. We mark by dots the points corresponding to a/M = 
0.0, 0.1, 0.2, ... 1.0. Right: trajectory of the first eight Kerr QNM frequencies 
with m = 2 (solid lines) and m = —2 (dashed lines). Filled circles mark the 
corresponding mode in the Schwarzschild limit. An arrow indicates the small 
loop described by the "exceptional" QNM with n = 6, that does not tend to the 
critical frequency for superradiance (see also Figs. 3-4 in [230]). The data used 
to produce this figure (and more) are available online I47| . 

We illustrate the splitting of the fundamental gravitational QNM with / = 2 
in the left panel of Figure [8] Even though QNMs have both positive and negative 
frequencies, it is customary to plot only the positive-frequency part of the spectrum 
[75j : in view of the symmetry properties listed above, this yields all the necessary 
information. As the rotation parameter a/M increases, the branches with m = 2 and 
m = —2 move in opposite directions, being tangent to the the branches with m = 1 
and m = — 1 in the limit a/M 0. For low overtone numbers and small values 
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of a/M the rotation-induced splitting of the modes is roughly proportional to m, as 
physical intuition would suggest. 

The weakly damped modes of Kerr black holes In the right panel of Figure[8]we show 
the first eight gravitational QNM frequencies with m — 2 (solid lines) and m = — 2 
(dashed lines). A general feature is that almost all modes with m > cluster at the 
critical frequency for superradiance, 2Muj = m, as a/M 1. This fact was first 
observed by Detweiler |262j , and a thorough examination of the extremal limit can be 
found in Refs. [2631 12641 1265j . The mode with n = 6 (marked by an arrow) shows a 
peculiar deviation from the general trend, illustrating the fact that some positive-m 
modes do not tend to this critical frequency in the extremal limit. 




a/M a/M 

Figure 9. Frequencies and quality factors for the fundamental modes with 
i = 2, 3, 4 and different values of m. Solid lines refer to m = I,.., I (from 
top to bottom), the dotted line to m = 0, and dashed lines refer to m = — 1, .., — Z 
(from top to bottom). Quality factors for the higher overtones are lower than the 
ones we display here. 

For gravitational wave detection we are mostly interested in the frequency and 
quality factor of the dominant modes, which determine whether the mode lies in the 
sensitive frequency band of a given detector and the number of observable cycles. 
Figure [3] shows these quantities for QNMs with ^ < 5. Improving on previous results 
[9l I266j . Ref. [TOj presented accurate fits for the first three overtones with ^ = 2,3,4 
and all values of m, matching the numerical results to within 5% or better over a range 
of a/M G [0,0.99] (see Tables VIII-X in Ref. [T^ and the numerical data available 
online [47]). For instance, defining 5=1 — a/M, the frequency ujim = i^R and quality 
factor Qim = ljr/{2lui) of the fundamental I — m = 2 and ^ = 2 , to = modes are 

Mw22 ^ 1.5251 - 1.1568 S^-^^^^ , Q22 ^ 0.700 + 1.4187 b°-^^^° , (96) 
Mwao - 0.4437 - 0.0739 6°-^^^° , Q20 - 4.000 - 1.9550 S" "^" , (97) 

The highly damped modes The intermediate- and large-damping regime of the QNM 
spectrum of Kerr BHs is even more puzzling than the RN spectrum. The main 
technical difficulty in pushing the calculation to higher damping is that Leaver's 
approach requires the simultaneous solution of the radial and angular continued 
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fraction conditions. For mode order n > 50 the method becomes increasingly 
unreliable. A way around this "coupling problem" is to study the asymptotic behavior 
of the angular equation as \auj\ ~ \aujj\ — oo |225| 1153] , The leading-order behavior 
of the separation constant sJ^imio-i^) when lo ~ iuji and \auj\ — »■ oo is 

sAhn = i{2L + l)acj + 0(1) , \aLj\ oo , (98) 

with L = I — \m\ for \m\ > \s\ and L = I — \s\ otherwise [2251 1153] . By replacing 
this analytic expansion in the radial continued fraction one effectively decouples the 
angular and radial continued fractions, and the calculation of QNM frequencies with 
n 3> 50 can be performed by solving only the radial continued fraction. It turns out 
that in the highly damped limit uji^ is independent of (s , I) and proportional to m: 

2MujR = mzu{a/M) . (99) 




Figure 10. Left: asymptotic real part 2Mujn = 2ro(a/M) of the I = m = 2 
gravitational and scalar QNM frequencies extrapolated from numerical data. 
Points refer to numerical results, the line is the analytical prediction. Results 
are independent of I, s and the sign of m. Right: same for imaginary part. 



Numerical results for w{a/M) are shown in Fig. [TOl These numerical results 
have been confirmed by recent analytical calculations in an impressive tour de force 
|2671 1268] 12211 . The final result can be implicitly expressed as a contour integral, 
which in turn can be expressed as a sum of elliptic integrals j2681 1222] . The relevant 
equations are summarized in [Appendix B[ and compared against our own numerics 
in Fig. 1101 The agreement is remarkable, given how involved the numerical and 
analytical calculations are. In the Schwarzschild limit [a/M 0) it can be shown that 

1/3 

■cu{a) ~ 0.30634 (a/M) ' and 2MSloi ~ 1/2. A good fit to the analytical predictions 
for the real part (accurate to within 0.8% in the entire range) is 

ZD = O.m (a/M)^^^ -0.308 a/M+0.156 (a/Mf -0.052 {a/M f .(100) 

In the extremal limit {a/M 1) the analytical results imply m ~ 0.10341 , 2MSujj ~ 
0.51260. For any a/M, the imaginary part uji grows without bound. The spacing 
between consecutive modes Stoi is not simply given by 2ttT, but it is a monotonically 
increasing function of a/M. A power fit in a/M of the numerical results yields [225j : 

2MSlji = 1/2 + 0.0219a/M- 0.0089 (a/M)^ . (101) 
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The agreement between the analytical predictions for the mode spacing and the fit of 
Eq. (jlOip is better than 0.1% in the entire range of a/M . A generalization of these 
asymptotic results to higher-dimensional rotating BHs can be found in Ref. [222j . 

The eikonal limit The eikonal limit of Kerr QNMs is not yet fully understood. Partial 
results concern I — zLm modes, for which [73l [lOl 1208] 

LJ = ±mflc'- (n + 1/2) iricl y/SM/rcS-^ , (l = \m\ oo) , (102) 

where Tj, , flc are counter-rotating or co-rotating radius and orbital frequencies at the 
unstable circular null geodesies, and 26 = rdvc — M)/{r1 — 2Mrc + a^)- Again, this 
result can be expressed as ui ^ Iflc — i{n + 1/2)/t, with r the typical instability 
timescale of the unstable circular null geodesic [751 ITUl [208) . To our knowledge, a 
simple geometrical optics interpretation is still lacking for modes with / ^ |m|. 

5.4- Kerr- Newman 

General BH solutions of the Einstein-Maxwell system are described by the Kerr- 
Newman metric. For this metric the Klein-Gordon equation and the Dirac equation are 
still separable [269j , so QNMs can be computed using the same general methods that 
apply to Kerr BHs |256j . The scalar spectrum was analyzed by Berti and Kokkotas 
|233j , who showed in particular that the eikonal limit can still be understood in terms 
of unstable circular null geodesies. 

Unfortunately, studies of the interplay of electromagnetic and gravitational fields 
in the Kerr-Newman metric are plagued by a major technical difficulty: to date, all 
attempts to decouple the electromagnetic and gravitational perturbations have failed 
(see e.g. Section 111 of Ref. f2T). Approximations to gravitational/electromagnetic 
perturbations of the Kerr-Newman geometry either keep the geometry fixed and 
perturb the electric field, or (more interestingly) keep the electric field fixed and 
perturb the geometry. This approach should be appropriate for values of the charge Q 
at most as large as the perturbations of the spacetime metric. QNMs for gravitational 
and electromagnetic perturbations in this approximation scheme were computed in 
Refs. j 270[|233] . but a solution of the general problem is highly desirable, and it could 
shed light on many open problems in BH physics. 

5.5. Higher dimensional Schwarzschild-Tangherlini black holes 

Here we briefly discuss the QNM spectrum of the higher-dimensional analogs of the 
Schwarzschild solution, known as Schwarzschild-Tangherlini BHs. Electromagnetic 
perturbations of these BHs were considered by Crispino, Higuchi and Matsas |134| . 
The elegant work by Kodama and Ishibashi [1021 1103[ I104j laid the foundations for 
the analysis of higher-dimensional RN BHs. 

Weakly damped modes The lowest lying QNMs of a d— dimensional Schwarzschild- 
Tangherlini BH were computed in a WKB approximation in Refs. [2T6[ [27T[ I?/?, '204) . 
Leaver's method can be generalized to these higher-dimensional BHs [234, 235], 
although the number of terms in the recurrence relation rapidly grows with d. 
For instance, in c? = 5 the recurrence relation for vector-gravitational and tensor- 
gravitational perturbations has four terms, while for scalar-type gravitational QNM 
frequencies it has eight terms [235 . A naive application of Leaver's method breaks 
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down for d > 9. For large d more and more singularities, spaced uniformly on the circle 
\r\ = r/i (where rh is the horizon radius) approach the horizon at r = r^- A solution 
satisfying the outgoing wave boundary condition at the horizon must be continued 
to some mid point, and only then can the continued fraction condition be applied 
|273| . Alternative calculations of weakly damped modes for scalar and gravitational 
perturbations in d dimensions make use of time evolutions [274j . 

Highly- damped modes In the large-damping limit, the leading-order result of the 
monodromy calculation, Eq. (|86p . generalizes to the c?-dimensional Schwarzschild- 
Tangherlini metric. This was first suggested in Ref. ^lOlj , and then it was explicitly 
shown by Birmingham [275] (see also [276p . The perturbative technique of Ref. |245| 
has been extended to bosonic fields in higher-dimensional Schwarzschild-Tangherlini 
BHs, with the result [234 

uj Tln{l + 2 cos nj) + {2n+l)niT + kdu'J^'^^^^'^'^^'^^ . (103) 

Here j — Q for scalar fields and tensor- type gravitational perturbations, j — 2 for 
vector- type gravitational perturbations and j — 2/{d — 2) , 2 — 2/ (d — 2) for vector- 
and scalar- type electromagnetic perturbations, respectively. The coefficient kd can 
be determined analytically for given values of d and j [234]. For electromagnetic 
perturbations in rf = 5, Eq. (|103p is singular: this cither means that there are no 
asymptotic QNM frequencies at all (a possibility first suggested in [1011) or that the 
asymptotic frequency is zero. Numerical results support the latter possibility |234| . 

The result above illustrates a general feature of the high-damping regime, which 
concerns the damping itself. It was shown by several authors [2771 1278|, I249j that for 
spacetimes with a single horizon, for all values of I one has 

ujj — ~2'Kin/T , rt — > cxD , (104) 

where T is the Hawking temperature. General results for the oscillation frequencies 
were obtained in Ref. [279j for single-horizon spacetimes (see also |280| for BHs in 
two-dimensional dilaton gravity) . An elegant unification of these results dealing with 
a rather general class of wave equations is given in Ref. [281j (but see Ref. [282j for 
words of caution on the generality of the results). A review on the highly damped 
regime of QNMs for several spacetimes can be found in Ref. [107j . 

Eikonal limit For I ^ 1, a WKB analysis of the Schwarzschild-Tangherlini 
perturbation equations yields j2161 f208j , 

ujr+ = nc{l + d/2~3/2-i{n+l/2)Vd~^^ . (105) 

where rJc — Vd— 3/\/d— 1 [2/{d— 1)]"^ is the orbital angular velocity at the circular 
null geodesic pilSj. 

Infinite- dimensional limit (d — s- oo) To our knowledge, the interesting limit where the 
dimensionality d tends to infinity has not been studied in much detail in the literature. 
Using a WKB analysis for scalar fields or gravitational tensor modes we find 

ujr+^d/2~iK^/d/2, d ^ oo , (106) 

with K a factor of order unity depending on the perturbing field (k = 1, l/\/2 for scalar 
fields and vector gravitational perturbations, respectively). BHs in higher dimensions 
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are much better resonators, with a quality factor Q = loji/\^uji\ increasing as ^/d. As 
far as we know, no numerical studies are available to confirm this analytical prediction. 

It is important to recall that black objects in higher dimensions may have other 
topologies beside the spherical one. One family consists of the so-called black strings, 
whose gravitational perturbations were studied by Gregory and Laflamme |283[ 1284] 
and Kudoh [285j . Higher-dimensional "squashed" Kaluza-Klein BHs were considered 
in |286j . and perturbations of brane-localized BHs were studied, for instance, in 
[2871 [2881 [2891 . Apart from specific noteworthy exceptions [2901 IMl |292l [293] . 
gravitational perturbations of higher-dimensional rotating solutions are yet to be 
understood. For a review on higher-dimensional BHs, see Ref. J4J. 

6. The spectrum of asymptotically anti-de Sitter black holes 

6.1. Schwarzschild anti-de Sitter black holes 

The metric ^ of an uncharged static BH in an (asymptotically) AdS spacetime has 
two parameters, vq and L, related, respectively, to the mass of the BH and the 
cosmological constant. The horizon radius r_|_ is the largest root of the equation 
/(r) = 0. The Hawking temperature associated with the metric ^ is 

T=[d-3 + {d- l)e] I (4<L) , (107) 

where ^ = f+/L. For d > 4, this fmiction has a characteristic minimum at 
Cmin = v/(rf-3)/(d- 1) with T„,i„ = ^{d-l){d-i)/2nL (see Fig. [Til). 



1.0 
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Figure 11. The (dimensionless) Hawking temperature TL as a function of tlie 
(dimensionless) horizon radius r+/L for asymptotically AdS BHs (shown here for 
d = 5). The horizontal line is the critical temperature of the Hawking-Page phase 
transition. 



The specific heat of BHs with ^ < ^^in (^ > ^min) is negative (positive). For 
T < Tmin, no BH solution exists, and the state with minimal free energy is "thermal 
AdS" (Euclidean AdS/Z). For T > T^-nim there are two solutions with the same 
temperature, "small" and "large" BHs. Small BHs have negative specific heat and 
are thermodynamically unstable. Large BHs exist in equilibrium with the heat bath. 
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Moreover, for < 1, the minimum of the free energy still corresponds to thermal AdS, 
even for T > Tmin- Thus the stable ground state is given by large BHs with ^ > 1 
(and T > Tc — {d — 2)/27rL). The (first order) transition between thermal AdS and 
a large AdS BH is known as the Hawking-Page phase transition ^94] . As shown by 
Witten [295j , in the gauge-gravity duality this transition corresponds to a first order 
deconfinement transition in a dual thermal strongly coupled gauge theory on a sphere 
of radius L. These considerations can be extended to charged and rotating AdS BHs 
(see e.g. |296| 1297] and references therein). The role of small AdS BHs in the gauge- 
gravity duality is not well understood (see, however, Ref. [298p . The limit of extra 
large BHs with r^/L — )■ oo in the appropriate scaling leads to black branes dual to 
strongly coupled gauge theories in flat space (see Section [8]). 

BHs in AdS in the context of the gauge-string duality attract considerable interest 
because they serve as a good laboratory for studying the most acute problems of a 
theory of gravity: the information loss paradox, BH singularities and some aspects 
of quantum gravity. Although we will not venture into this fascinating field in this 
review, in Section [8] we mention some of these topics which involve the use of QNMs. 

Another peculiar feature of asymptotically AdS space is the "active role" played 
by its boundary. In AdS, null geodesies reach the boundary in finite coordinate time. 
One thus often refers to an asymptotically AdS space as a box, having in mind that 
AdS boundary conditions directly affect the bulk physics [2991 13001 1301] . This should 
be contrasted with the asymptotically flat case, where the only physically relevant 
choice for the boundary conditions of the bulk fluctuations corresponds to outgoing 
waves at spatial infinity. In the gauge-gravity duality, the choice of the AdS boundary 
conditions is dictated by a holographic prescription [98l I302| I303| I304j , see Section [51 

QNMs in asymptotically AdS backgrounds were first considered by Chan and 
Mann for a conformally coupled scalar field [305j (see also Ref. [306]). Subsequent 
interest was strongly motivated by the development of the AdS/CFT correspondence. 
An interesting albeit somewhat qualitative discussion of thermalization in AdS/CFT 
first mentioning QNMs appeared in Ref. i30Z] (see also the important work of 
Danielsson, Keski-Vakkuri and Kruczenski [3D51 I309j ). Horowitz and Hubeny f52] 
explicitly pointed out the fundamental link between QNMs of a large AdS BH, which 
describe the background's relaxation into a final state, and the dual field theory, where 
they describe the approach to thermal equilibrium. Ref. [92] computed the QNMs of 
a minimally coupled massless scalar in intermediate or large AdS BHs for dimensions 
d = 4,5, 7. After this seminal work, a series of studies in four and higher dimensions 
led to a deeper understanding of QNMs in asymptotically AdS backgrounds. 

The following is a brief overview of the QNMs of the four and higher-dimensional 
SAdS geometries. For more details we refer to the original works (see 92, 107, 130, 131] 
[3Tf)l [T941 [3TT1 [3T21 [3T31 [2231 [3T4l [3T5l [3T6l [3171 [3T8l [3T9l [3201 [27T1 132L ,322. .323. .324] 
and references therein). The results below focus mostly on the large BH regime 
r+ / L ^ 1 and on Dirichlet boundary conditions, which are more relevant for 
holography. Other boundary conditions were investigated in Refs. |162[|163Tll64[[T65] . 

The weakly damped modes of large black holes The fundamental QNM frequencies 
for scalar field perturbations were first computed by Horowitz and Hubeny |92] 
for intermediate and large BHs. In Refs. [130] I194j the analysis was extended 
to electromagnetic and gravitational perturbations. These works considered only 
d = i; the d = 5 case was analyzed in Ref. |316] . and half- integer spins were 
considered in Refs. |325[ I326j . Approximate analytical solutions have been discussed 
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in Refs. |3121 1313t 1314] with particular emphasis on the d = 5 geometry. 
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Figure 12. QNM frequencies for scalar [s = 0), electromagnetic (s = 1) and 
gravitational (s = 2) perturbations of large SAdS BH, computed for r+/L = 100 
and I = s. In the large BH regime the frequencies scale with r+/L and are 
basically independent of I for I r^/L. Furthermore, scalar-field, tensor-type 
and vector-type gravitational perturbations are nearly isospectral, except for a 
special mode belonging to the vector-type family and marked by a red diamond. 
For further details see [TM1I5TT] . 



The weakly-damped QNM spectrum of a large (r+/L = 100) SAdS BH in 
d = 4 is shown in Fig. [T^l Usually the QNM frequencies scale with the horizon 
radius, uj cx r+fL"^ [92l 1194] . so the modes of any large BH can be obtained 
by rescaling appropriately the numbers in the figure. The exception to this rule 
is marked by a red diamond and will be discussed in more detail below. The 
frequencies of different perturbations are very similar: scalar perturbations and scalar- 
type gravitational perturbations have basically the same spectra, while vector-type 
gravitational perturbations are displaced by one overtone relative to these two |194| . 
The QNM frequencies are practically independent of Z for I <^r+/L; the large-Z limit is 
discussed below. The QNM spectrum for electromagnetic perturbations has a peculiar 
structure [1301 13111 11941 13271 1328j : the real part of some modes asymptotes to zero 
when r_|_/L — > oo, and for the dominant mode uJiL'^/r^ —^ —1.5. 

The fundamental QNM for vector-type gravitational perturbations is extremely 
long-lived when compared to all other modes of other kinds of perturbations. This 
mode has an interesting interpretation in the gauge/gravity duality, discussed in 
Section [H The timescale of this long-lived mode is proportional to r+, and the mode 
itself is easily computed numerically by a straightforward application of the series 
solution [13011194] . In a general d— dimensional SAdS geometry it is well described by 

r+ w™=*°'-*yP^ = - + d- 2)/{d - 1) , r+/L » 1 . (108) 

This was first observed numerically in Refs. [1301 1194] for d = 4, and in Ref. [316j for 
d = 5. The four-dimensional result was later confirmed analytically by Miranda and 
Zanchin [315 . Eq. ([lOSp for general d was derived by Siopsis [324j . and we checked 
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that it agrees with numerical results for d = 4 up to d = 8. In d = 4 and for r+/L ^ I, 
corrections to Eq. ()108|) can be found: 

r+ w^T^*''''' = -Kl - + 2)/3 - 0.0288 ^ l^L^rl . (109) 

We estimate the uncertainty in the term to be about 5%. Analytical calculations 
of these corrections have been done in the r_|_ ^ oo limit [109|, I329|, I330j and are 
consistent with the numerical results above after proper identifications. In this limit 
the geometry becomes that of a black brane and the angular wavefunctions Yim are 
replaced by e*^^. For large l,p one has the correspondence Lp = I |331| . In particular, 
it is found that in the large l/L limit, 



o,49-91og3 + VStt iP 0.0281^2^4 

-y-^i^ 16H T (110) 



in quite good agreement with the numerical fits. For similar high-order analytical 
corrections in higher-dimensional AdS backgrounds, see Refs. [1091 13291 1330] . 

If Dirichlet boundary conditions are imposed at infinity, the scalar-type 
gravitational sector does not have such a long-lived mode. It was suggested that the 
preferred boundary conditions in the AdS/CFT framework are of Robin type, and in 
particular that the perturbations should not deform the metric on the AdS boundary 
|316| I163j . Using Robin conditions, a long-lived mode for gravitational perturbations 
was discovered in Refs. j316| I163[ I324j . 



The weakly damped modes of small black holes The series solution method, which 
works so well for large BHs, converges very slowly for small BHs |318[ 13101 I194j . 
Recent results make use of Breit-Wigner type resonances in the scattering cross-section 
to study the very small BH regime [227j . described briefly in Section 1431 Small and 
large BHs have a very different behavior from a QNM perspective. This is related 
to the qualitatively different behavior of the potential in the two regimes, which is 
shown in the left panel of Fig. [T3] for s = 0: for small BHs {r+/L < 1) the potential 
develops a well capable of sustaining trapped, long-lived modes, corresponding to 
quasi-stationary states in quantum mechanics |332[ 13311 1227] . In the right panel of 
Fig. [13] we plot the QNMs of scalar fields for different BH sizes and 1 = 0, 1. In this 
limit (and under the assumption that Mlor ^1) it is possible to prove that [333] 

lurL ~ I + 3 + 2n- kinr+/L , (111) 

ojiL- -jo{l, n){l + 3 + 2n) {r+/L) ^'+^ /n . (112) 

The constant jo{l,n) {n = 0, 1, 2, . . .) and selected values of the constants fc;o can 
be found in Ref. 227J. The analytic prediction is consistent with numerical results. 
Likewise, for gravitational perturbations with r+/L ^ 1 one finds 

ujRL^l + 2 + 2n-cinr+lL, (113) 

LOiL- -72(;, n){l + 2 + 2n) {r+/L) , (114) 

where the constants q„ and 72(^1 n) can be found in Refs. ]3341I227] . In conclusion: for 
scalar field and gravitational perturbations of small SAdS BHs the damping timescale 
ujjL ^ (r+/_L)2'"'"2, and the oscillation frequencies approach the pure AdS value in the 
limit r+/L ^ (see [3101 1227] for more details). To our knowledge, generalizations of 
these results to charged, higher-dimensional and rotating AdS BHs are still lacking. 
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Figure 13. Left: Potential for scalar field propagation of ( = modes in a 
SAdS background, for different values of the BH size r_|_/L. A local maximum 
(and a potential well) only exist for small BHs. Right: Track described by the 
fundamental scalar field QNMs with I = and / = 1 as we vary the BH size r+ /L. 
Counterclockwise starting from the top-right of the figure we marked the points 
corresponding to different decades in rj^jL {tj^/L = 10'^, 10^, lO'', 10~^, . . .). 
Modes with different I coalesce in the large BH regime, as long as Z <C rj^/L. 



The highly-damped modes The large-damping regime (n ^ 1) can be studied using 
the monodromy method for both small and large AdS BHs [1071 1223] . For n ^ 1 and 
r+/L ^ 1 (large BHs) the QNMs are given by 

^ = (d-l)sin(-^) ei 

Results for arbitrarily sized BHs can be found in Refs. [107[|223] . Eq. (|115p describes 
scalar fields, gravitational tensor and gravitational vector-type perturbations. It 
yields ujL'^/r+ « (1.299 + 2.250i)n + (0.573 + 0.419i) for d = 4 and wi7r+ w 
(2 + 2i)n + 3.22064 + 2.77936 for d = 5, in good agreement with numerical results for 
d = 4 [194j and d = 5 |316j . For large SAdS BHs the asymptotic high-damping regime 
is approached very quickly, so Eq. (jllSp describes fairly well even weakly-damped 
modes. There is some disagreement concerning the offset for scalar-type gravitational 
perturbations |107j . 

An asymptotic analysis based on the monodromy method was recently used to 
predict the existence of a new family of modes [335] . So far, this new part of the 
spectrum has not been confirmed by numerical studies. In Ref. I336J the authors 
discuss the interesting possibility of describing the highly damped QNM regime in 
terms of geodesies of the SAdS spacetime. 



1 .log2 
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(115) 



The eikonal limit of large black holes In asymptotically AdS spacetimes the eikonal 
limit is especially interesting, since large-/ modes can be very long-lived [SH I331j . 
A WKB analysis [331j and numerical investigations |337j show that for scalar field 
perturbations, r+/L':^ 1 and I 3> 1 

= / -f H„ (r+/L)^ , (116) 
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So large-? modes are very long-lived, and they could play a prominent role in the BH's 
response to generic perturbations. This is at variance with the asymptotically flat 
case, where the damping timescale is roughly constant as / varies. Notice also that the 
scaling with the BH size differs from that of the weakly-damped and highly-damped 
modes. Other types of perturbations also display a similar qualitative behavior [337j . 



The eikonal limit of small black holes For scalar and electromagnetic perturbations 
of small BHs the potential for wave propagation for I 3> 1 develops a minimum, 
potentially supporting quasi-stationary states, i.e. long-lived modes (see Fig. [T5)) . 
Define rf, > Tc to be the two real zeros (turning points) of uj^ ~P^f/r'^ = 0. Then the 
real part of a class of long-lived modes in four spacetime dimensions is given by the 
WKB condition 

^^^^'^^^*^,('2n + l + 2V ,117) 



rf V 2 

where p = I + d/2 — 3/2. Their imaginary part is given by 



uJi = ^. logr = 2i/ ^ E—L±dr. 118 

The prefactor 7, not shown in Ref. [331j . can be obtained by standard methods as 
shown in Ref. [227], where these results are also supported by numerical calculations. 



6.2. Reissner-Nordstrom and Kerr anti-de Sitter black holes 

The analysis of QNMs of large RNAdS BHs was performed in Ref. [3351 1333 for 
weakly damped modes of a massless scalar field, and later extended to charged scalar 
fields [340]. Scalar field, electromagnetic and gravitational perturbations of RNAdS 
BHs were analyzed and compared in Ref. [311J. Half-integer spins were studied in 
Refs. [3411 1342] , Gravitational perturbations of higher-dimensional charged solutions 
were considered in [343] . and the highly-damped regime was explored analytically by 
Natario and Schiappa |107j . 

Ref. 3IIJ pointed out some interesting facts: (i) a near-isospectrality of different 
classes of perturbations holds for large BHs, but it breaks down as r+/L decreases, 
i.e. in the small BH limit; (ii) the imaginary part of the purely damped modes found 
for electromagnetic and vector-type gravitational perturbations tends to zero as the 
charge Q tends to the extremal value Qext, possibly pointing to a marginal instability 
of extremal RNAdS BHs; (iii) for all kinds of perturbations, the real part of the 
fundamental QNM frequency, Lojr^ has a minimum for Q/Qcxt — 0.366, followed 
by a maximum at Lojr ~ 0.474. A reanalysis of massless scalar field perturbations 
found that the imaginary part of scalar QNMs also tends to zero in the extremal 
limit [344j . However, comparing time-evolutions of the field with results from the 
Horowitz-Hubeny series solution, Wang et al. |344| found that the BH response turns 
from a standard, oscillatory QNM-type decay below some critical value of the charge 
{Q < Qcrit — 0.3895(3ext) to a non-oscillatory behavior characterized by a purely 
imaginary QNM frequency for Q > Qcrit- Furthermore, these authors suggested that 
the potential marginal instability proposed in Ref. [311| does not pose a threat for 
extremal RNAdS BHs, because (at least for scalar perturbations) the asymptotic field 
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in the near-extremal limit is dominated by a power-law tail of the type first analyzed 
by Price [171j . The implications of these findings in the context of the AdS/CFT 
correspondence deserve further investigation. 

There are few studies of QNMs of rotating BHs in AdS backgrounds. The 
formalism to handle perturbations of these spacetimes was laid down by Chambers 
and Moss [345]. Giammatteo and Moss |346j considered axially symmetric QNMs 
of large Kerr-AdS BHs. These spacetimes behave as a BH in a box. Superradiant 
amplification of incident waves at the expense of the BH's rotational energy can 
then produce instabilities in the non-axisymmetric modes of small Kerr AdS BHs 
[3471 [3331 [3341 [3481 [3491 [292 [350|: this is an interesting example of the "black hole 
bomb" first investigated by Press and Teukolsky [351j . 

6.3. Toroidal, cylindrical and plane- symmetric anti-de Sitter black holes 

The uniqueness theorems that apply for asymptotically flat BHs (see e.g. ^39j) can be 
evaded when the cosmological constant is non-zero: solutions have been found with 
the topology of a cylinder, of a plane or of a doughnut [122 (HSl [TM [1251 [1261 [T27] . 
Defining /(r) = r^/i^ — AM/r, these spacetimes are described by 

ds^ ^ fdt^ - r^dr^ - d?? - r^dcf)"^. (119) 

The range of the coordinates z and (j) dictates the topology of the BH spacetime. For 
a BH with toroidal topology the coordinate z is compactified such that z/L ranges 
from to 27r, and cj) ranges from to 27r as well. For the cylindrical BH, or black 
string, the coordinate z has range —oo < z < oo, and < < 27r. For the planar BH, 
or black membrane, the coordinate (j) is further decompactified {—oo < L cj) < oo). For 
the torus, M is related to the system's ADM mass; for the cylinder, to the mass per 
unit length of a constant-z line; and for the plane, to the mass per unit area of the 
{(j),z) plane [l22l [T23l [T24l [1251 FIM [T27] . 

A formalism to handle electromagnetic and gravitational perturbations of these 
spacetimes was developed and used to investigate numerically the QNMs of scalar, 
electromagnetic and gravitational perturbations in Ref . [131) . A thorough analysis by 
Miranda and Zanchin [315[ I327 j confirmed and extended the results of Ref. ^131] . 

7. The spectrum of asymptotically de Sitter and other black holes 

7.1. Asymptotically de Sitter black holes 

BHs in a dS background have a rich structure. Consider an uncharged spherically 
symmetric BH, the SdS^ solution, 

ds^ = fdt^ ~ f-^dr^ - r'^dn}^^ , (120) 

where /(r) — 1 — r^/i^ — rp^'^/r'^^^, dr2^_2 the metric of the (d — 2)— sphere, and tq 
is related to the mass M of the spacetime by M = {d— 2) Ad- 2rQ~^ / (IGtt). Depending 
on the value of ro/L'^, in general this solution has two horizons, the event and the 
cosmological horizon. When the cosmological and event horizons coalesce one has a 
so-called extremal SdS^ solution, which is related to a topologically different solution, 
known as the Nariai spacetime [TM [MH [^ [TUS] . 
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Weakly damped modes The first calculations of the fundamental gravitational modes 
were done by Mellor and Moss |129j for Reissner-Nordstrom-de Sitter (RNdS) BHs, 
using numerical techniques. Brady et al. [352j complemented this study through 
a numerical time evolution of scalar fields in the SdS spacetime. Time evolutions 
were also performed in Ref. [353j . where the results were compared against WKB 
predictions (see also Ref. |274j ). Approximations to the correct, lowest order QNMs 
were considered in Refs. [I62] (s = 2, SdS), 354J (s = 0, SdS), [355] (s = 1/2, RNdS) 
and [356j (s = 2, RNdS), where the true potential was approximated by a Poschl- 
Teller potential. Cardoso and Lemos [191] showed that the true potential reduces 
to the Poschl- Teller potential for near extremal SdS geometries, which explains why 
previous results based on the Poschl- Teller approximation gave accurate predictions 
for the QNMs. The results in [191j were later generalized to d-dimensional RNdS BHs 
|197j . Other analytical results used the WKB approximation for spins s — 0, 1, 2, 1/2 
in the vicinity of SdS BHs [557j . 

Highly- damped modes The highly-damped limit, where the imaginary part is much 
larger than the real part, was studied numerically by Yoshida and Futamase |238j for 
near-extremal uncharged BHs (s = 1, 2) and in the general case in Ref. [35811239] . The 
results were confirmed analytically in Refs. [223j (s = 0, 1, 2 SdS) and [107] (s = 0, 1, 2 
RNdS). 

Eikonal limit In the large-Z limit, WKB techniques [190|, 1357] yield 

sVsMlj = Vl - 27M2/i2 ^1^1/2- i{n + 1/2)] . (121) 

QNMs in the eikonal limit can be interpreted as perturbations of unstable circular 
null geodesies [208]. For studies on charged and rotating BHs in dS backgrounds, see 
|359] and references therein. 

7.2. Black holes in higher- derivative gravity 

Among all possible theories of gravity with higher derivative terms, theories modified 
by the addition of a Gauss-Bonnet term TZgb = Rp.upaR'^'''"' - '^Rf^uR^" + R^ , are 
particularly attractive and have been considered by many authors ^360] . There are 
simple BH solutions in these theories [MD IMl [M3 ISMl [Ml [Mi [MH [Ml [3M1 
1370) . whose perturbations were studied in a series of works by Dotti and Gleiser 
[5711 [371 [M], Moura and Schiappa [M] and Takahashi and Soda [371]. In four 
dimensions the Gauss-Bonnet term is a total divergence, and yields upon integration 
a topological invariant, being therefore equivalent to Einstein's theory. Therefore 
these theories are interesting in higher dimensions only. Non-trivial four-dimensional 
scenarios can be accomplished by coupling the Gauss-Bonnet term to a dilaton 
P67": the resulting theory also arises from the low-energy limit of certain string 
theories [360] . Perturbations of BH spacetimes in these scenarios were considered 
in Refs. [3671I3751I365] . 

QNMs of BHs in these theories were first investigated by Iyer et al. [376] , who 
used a WKB approach to study scalar perturbations. QNMs of these spacetimes were 
also studied through a WKB approach [3771 ISlSl [379], and numerically [3801 l377] 
for the low-lying modes. The highly-damped regime was analytically explored in 
Ref. [282j . The eikonal limit was considered by Konoplya |380] and interpreted in 
terms of circular null geodesies in Ref. [208j . 
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1.3. Braneworlds 

Braneworld scenarios, where the standard model lives in a four-dimensional brane 
embedded in a higher dimensional spacetime, have been a popular research topic 
in the last decade. The extra dimensions can be compact [38 1[ I382j or even 
infinite, fiat [3831 or curved [384) 1385] . BH solutions in these theories are extremely 
difficult to find (see for instance Ref. |386|, 1387] for a discussion), some solutions are 
known perturbatively in some regimes. For instance, in the case of flat compact 
extra dimensions, Tangherlini BHs should be a good approximation to a static BH 
solution as long as the horizon radius is much smaller than the size of the extra 
dimension [386]. In these scenarios the standard model is localized on the brane and 
therefore BH oscillations are non-trivial, especially for QNMs of standard-model fields 
[M81IM91 [3901 [39T1 [3921 [393] . 

7.4- Black holes interacting with matter 

Astrophysical BHs are not expected to be in complete isolation, so it is important to 
understand how QNMs change when BHs interact with the surrounding environment. 
Leung et al. }394[ I395j investigated how the low-order QNMs of a BH are affected by 
a small amount of matter. For such a BH, in the static, spherically-symmetric case, 
Medved et al. |277U278j proved that highly-damped QNMs depend only on the surface 
gravity. Specific models for BHs interacting with matter were constructed by several 
authors. A popular cosmological scenario invokes the existence of dynamical vacuum 
energy ("quintessence", see e.g. j396[l397] ) or phantom fields |398ll399j to explain the 
acceleration of the universe. The QNMs of BHs with quintessence or phantom fields 
have been investigated in Refs. [1001110111101110111011105] and gOO], respectively. 

8. Quasinormal modes and the gauge-gravity duality 

In this Section, we review a particular entry in the gauge-gravity duality dictionary 
directly related to QNMs. It turns out that quasinormal spectra of asymptotically 
AdSd-i-i and more general backgrounds correspond to poles of the (retarded) thermal 
correlators of dual d-dimensional strongly interacting quantum gauge theories. The 
lowest quasinormal frequencies of black branes have a direct interpretation as 
dispersion relations of hydrodynamic excitations in the dual theory. More information 
on near-equilibrium properties of BHs and black branes and their holographic 
interpretation can be found in the reviews d 11171 [5S1 [571 [55] . 

8.1. The duality 

The discovery |407] and subsequent studies of _D-branes in string theory led to 
the concept of gauge-string duality. In the original example of the duality, 
known as the AdS/CFT correspondence [6j, the full type IIB string theory on 
the background AdSs x (five-dimensional AdS space times a five-sphere) was 
conjecturec|§| to be equivalent to the specific supersymmetric gauge theory, M = A 
SU [Nc) supersymmetric Yang- Mills theory (SYM) in fiat four-dimensional spacetime. 
The equivalence is understood as an equivalence of quantum partition functions. In 

§ During the last decade, the AdS/CFT correspondence and, more generally, the gauge-string duality, 
survived numerous, often very non-trivial tests of validity. At the moment, there is very little doubt, 
if any, that the conjecture is valid. 
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quantum field theory, there is very strong evidence supporting the claim that M ~ A 
SYM is a conformal field theory (CFT), with the beta- function identically equal to 
zero and coupling constant qym being independent of the energy scale. The two 
parameters characterizing A/" = 4 SYM are the 't Hooft coupling A = gy^-^c and the 
number of colors Nc- In the AdS/CFT correspondence, these parameters are mapped 
into the string theory parameters L and gs'- 

9ym Nc ~ gsNc - L^lt , Nc ^ 1^1% , (122) 

where L is the parameter of AdSs and the radius of the five-sphere, gs is the string 
coupling and IstIp are the string and Planck lengths, respectively. The full quantum 
string theory on AdSs x is poorly understood. However, its low-energy limit, type 
IIB supergravity, has been extensively studied since the 1980s. Restricting duality to 
the supergravity limit of the full string theory restricts the values of the gauge theory 
parameters to gyn ^ 1, -/Vc ^ 1. Thus the gauge theory at large values of coupling 
and large Nc is effectively described by classical gravity in the AdS background. 

Following the original example of AdS/CFT correspondence, many more dual 
pairs have been discovered, including those involving non-supersymmetric and non- 
conformal theories. Gauge-string duality thus includes the original AdS/CFT 
correspondence and all its "non-conformal" and "non-AdS" generalizations, often 
commonly referred to as "AdS/CFT". The gauge-string duality in the supergravity 
approximation is known as the gauge-gravity duality. The duality provides a 
quantitative correspondence between classical gravity in ten (or five) dimensions and 
a gauge theory (in the limit gyM -^c ^ 1, Nc ^ 1) in flat four-dimensional spacetime. 
Such a correspondence between higher-dimensional gravity and lower-dimensional non- 
gravitational theory is often referred to as "holography" . The gauge-gravity duality 
serves as a quantitative example of the "holographic principle" proposed by 't Hooft 
and Susskind [40811409] . 

Since classical higher-dimensional gravity is holographically encoded into the dual 
gauge theory's properties, one may wonder about the gauge theory interpretation 
of the QNM spectrum. The short answer, conjectured in [94], established in [98], 
and further generalized in |109] and many subsequent publications, is that the 
QNM spectrum of the fluctuation 5(j) of a higher- dimensional gravitational background 
coincides with the location of the poles of the retarded correlation function of the gauge 
theory operator O dual to the fluctuation 5<j). In the rest of this Section, we elaborate 
on this statement and provide some explicit examples. 

The main ingredient of the gauge-gravity duality is the ten-dimensional 
(super) gravity background characterized by the values of the metric and other 
supergravity fields such as the dilaton, the axion and various tensor fields. The 
background fields must satisfy supergravity equations of motion. (In most cases, only 
bosonic supergravity fields are considered, thus eliminating the need for the prefix 
"super".) For example, the near-horizon limit of the black three-brane background, 
which is the basic ingredient of the AdS/CFT duality at finite temperature, consists 
of the metric 



ds^ - — 



-fdt^ + dx^ + dy^ + dz' 



t2 

-jf-'^dr^ +L^dnl, (123) 



where f{r) = 1 — rg/r^, and the Ramond-Ramond five-form field, 

F5^- — {l + *)dtAdxAdyAdzAdr, (124) 

1j 
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with all other fields vanishing. According to the gauge/gravity correspondence, the 
background (fT23l) - (fT24|) with non-extremality parameter tq and Hawking temperature 
T = tqIt:!? is dual to TV 4 SU{Nc) SYM at finite temperature T in the limit of 
Nc — > oo, gy^Nc — > oo. The A/" = 4 SYM is defined in Minkowski space with 
coordinates t,x,y,z. The fifth (radial) coordinate r of the dual metric p23p plays 
the role of the energy scale in the gauge theory (with the boundary at r — > oo 
corresponding to the ultraviolet in the gauge theory), and the five-sphere describes 
internal degrees of freedom associated with the i?-symmetry group SU{A) specific to 
theories with A/" = 4 supersymmetries. 

Quite often, the internal degrees of freedom are of less interest, and the 
background can be dimensionally reduced from ten to five dimensions. For the metric 
(jl23|) , the result of such a reduction is the five-dimensional Schwarzschild- AdS metric 
with translationally invariant horizon 



J.2 

d4 = j, 



'fdr + dx^ + dy^ + dz" 



^r'dr^ , (125) 



obeying Einstein's equations i?^^ — 2A/35p,y in a five-dimensional space with 
cosmological constaniliJ A = — G/L^. Thus the gauge-gravity duality often appears 
as a correspondence involving the five-dimensional (super)gravity bulk and the four- 
dimensional boundary gauge theory. One should always remember, however, that all 
five-dimensional fields and their fluctuations have ten-dimensional origin. 

As a remark, we note that the metric (jl23p with translationally invariant horizon 
can be obtained from the Schwarzschild- AdS BH metric 

dsl^-fdt' + — +r'dQl, f = i + — --l- (126) 

by rescaling r — > A^/^r, tq — > A^/^rg, t A~^/'^i, and taking the limit A — s- cx) while 
simultaneously blowing up the sphere 

dnj A"i/2 {dx^ + dy^ + dz^) . (127) 

Similar rescalings for more general metrics can be found in [296[ 1410] . This difference 
between black hole and black brane metrics leads to the fact that black brane QNM 
spectra are functions w = uj{q) of the continuous parameter q (rather than a discrete 
parameter Z), where q is the momentum along the translationally invariant directions. 



8. 2. Dual quasinormal frequencies as poles of the retarded correlators 

Quantitatively, the gauge-string (gauge-gravity) duality is the equivalence of the 
partition functions 

Zym [J] - (e- ^ ^<^'^'^)yM = ^string [J] ^ e~^--^''^ , (128) 

where the semiclassical approximation on the right-hand side corresponds to passing 
from the gauge-string to the gauge-gravity duality in the appropriate limit (e.g. in 
the limit Nc ^ oo, gyM^c ^ oo for A/" = 4 SYM). The equivalence (|128p means that 
the classical gravity action effectively serves as a generating functional for correlation 
functions of gauge-invariant operators O in the dual gauge theory. On the gravity 
side, the role of J for a given operator O is played by the boundary value 50o of 

II The cosmological constant in five dimensions arises as a result of the dimensional reduction of the 
five-form II124I I on 5^. 
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the background fluctuation 5(j) (for a moment, we are ignoring all indices the field 5(j) 
might have). For example, the boundary value of the background metric fluctuation 
hfiv plays the role of J in computing the correlators of the energy-momentum tensor 
in a four-dimensional gauge theory (see Table [5]) . 

The recipe for applying the equivalence (|128p is the following. To compute the 
correlators of a gauge-invariant operator O, one has to 

• identify the dual fluctuation fleld 5(j) associated with O 

• solve the linearized bulk equations of motion satisfled by (50 with the boundary 
conditioij^ Scj) ScjjQ = J 

• using this solution, compute the on-shell supergravity action S'grav['^] as a 
functional of 6(I)q = J 

• compute the correlators in the usual field theory sense by taking functional 
derivatives of exp {—SgmviJ]) with respect to J. 

The recipe given above is sufficient for computing Euclidean correlation functions 
from dual gravity. For Minkowski space correlators, there are subtleties resolved in 
[Ml [302 (recent work on the Lorentzian AdS/CFT includes gTIl WM)- R-ef. [M] 
contains an exact prescription on how to compute the Minkowski space two-point 
functions from fiuctuations of a dual gravity background. As a byproduct, the 
prescription establishes a one-to-one correspondence between poles of the retarded 
quantum field theory correlators and QNM spectra of the dual background. Indeed, 
the dual gravity fluctuation fleld 5(j){r, t, x, y, z) associated with the operator O, whose 
retarded two-point function we are interested in, satisfies an ordinary linear second- 
order differential equation with respect to the radial coordinate r. The fluctuation's 
dependence on the "usual" four-dimensional space-time coordinates t, x, y, z in the 
bulk is typically trivial, allowing one to Fourier transform with respect to them: 

<50(r, t, X, y, ^) = / 1^ ^-^^t+^^^ ^^(^^ q) _ (129) 

Note that at finite temperature the Lorentz invariance is broken, thus the components 
w, q of the four- momentum are independent variables. For theories with unbroken 
rotation invariance, the fluctuation will depend on the magnitude of the three- 
momentum q = |q|, and one can conveniently choose q along the direction z of the 
brane, with the four- momentum being given by [uj^ 0, 0, q). The ODE satisfled by the 
fluctuation S(f){r,LU,q) typically has singular points at the horizon r = vq and at the 
boundary r — > oo. The solution S(j){r,LU,q) satisfying the incoming wave boundary 
condition at the horizon can be written near the boundary as 

50(r, LO, q) = A{lu, q) r'^- (1 + • • •) + Biu, q) r-^+ (1 + • • •) , (130) 

where A+, A_ are the exponents of the ODE at r = cx) (these exponents are related 
to the conformal dimension of the operator O), and ellipses denote higher powers of r. 
In most cases, flelds can be redefined so that A-|_ > 0, A_ = 0. Applying the gauge- 
gravity duality recipe for Minkowski correlators, for the retarded two-point function 
of the operator O one finds 

B 

G^{uj, q) ^ —7 + contact terms . (131) 

% The second boundary condition on the fluctuation 54> is either a regularity condition (e.g. for zero- 
temperature global AdS space) or the incoming (outgoing) wave boundary condition. For metrics with 
horizons the incoming (outgoing) wave condition corresponds to computing the retarded (advanced) 
correlators in the boundary theory. 
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Table 2. The correspondence between the boundary gauge theory operators and 
the dual fivc-dimensional gravity bulk fields. 



Gauge theory operator O 


Dual gravity fluctuation 


energy-momentum tensor T^i, 
conserved current 
Tr 


metric fluctuation h^^, 
Maxwell field A^, 
minimally coupled massless scalar Lp 



Zeros of the coefficient A correspond to the poles of G^{uj,q). On the other hand, 
from the general relativity point of view, the condition ^ = (for A_ = this is 
just the Dirichlet boundary condition) defines the QNM spectrum of the fluctuation 
S(j). Thus, all the information about the poles of the retarded correlators of a quantum 
field theory with a gravity dual description is encoded in the QNM spectra of the dual 
gravity fluctuations. This statement is a useful entry in the gauge-gravity duality 
dictionary, since the poles of thermal correlators carry important information about 
transport properties and excitation spectra of the theory. Consider, for illustration, 
the relatively simple case of a two-dimensional CFT at finite temperature dual to the 
(2 -|- l)-dimensional BTZ BH background. The retarded two-point function of the 
operator of conformal dimension A = 2 at finite temperature is given by [98j 



9 9 

uj — 
47r2 



V- 1 



AttT 



+ ip\ l-i- 



(132) 



where tpi^) = r'(z)/r(z); we put — Tji and ignored the constant prefactor for 
simplicity. The correlator (jl32p has infinitely many poles in the complex frequency 
plane, located at 

cj„ ==±g-i47rT(n+l), n = 0,l,2,.... (133) 

These poles are precisely the BTZ quasinormal frequencies j305|, [95l I412j . 

The role of QNM spectra as frequencies defining relaxation times in the dual 
thermal field theory was realized and discussed in early publications [307[ [92l I309j . 
Birmingham, Sachs and Solodukhin [94j were the first to note explicitly that the QNM 
spectrum of the BTZ BH coincides with the poles of the retarded correlators of the 
dual (l-l-l)-dimensional CFT. With the appearance of the Minkowski AdS/CFT recipe 
for the retarded correlators ^SSj , the relationship between QNMs and the poles of the 
correlators has been established quantitatively first for the scalar |98], [99j and later 
for general fluctuations '100', '109]. In the case of non-scalar fluctuations, considering 
gauge-invariant combinations of fluctuating fields is especially useful [109,, although 
this is not the only possible approach [4131 1414j . 

A word of caution is necessary. As we have seen, asymptotically AdS spacetimes 
offer a variety of choices for the boundary conditions at spatial inflnity. Not all such 
choices produce QNM spectra which have a meaningful interpretation in the dual 
quantum fleld theory. In order to say that a computed QNM corresponds to a pole 
of a dual theory correlator, one has to analyse the bulk fluctuation along the lines 
leading to Eq. (|13ip to establish the precise form of the boundary condition. 



8.3. The hydrodynamic limit 

The most interesting results for QNM spectra with a dual field theory interpretation 
are obtained for five-dimensional gravitational backgrounds (note, however, the 
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growing body of work on the AdS-Condensed Matter Theory correspondence [37 l I38 j . 
where, for the purposes of studying (2 + l)-dimensional condensed matter systems, 
one is interested in (3+ l)-dimensional gravitational backgrounds). For example, the 
poles of the retarded thermal two-point function 

G';,,,Au;,q) = -t [ d^xe~^^-^eit) ([T^.(x), T;,.(0)])t (134) 



of the stress-energy tensor in four-dimensional A/" = 4 SU{Nc) SYM in the limit 
Nc ^ oo, gyj^jNc oo are given by the QNM frequencies of the gravitational 
perturbation /i^^ of the metric (|125p . By symmetry, the perturbations are divided into 
three groups |415[ I109j . Indeed, since the dual gauge theory is spatially isotropic, we 
are free to choose the momentum of the perturbation along, say, the z direction, leaving 
the 0(2) rotational symmetry of the (x, y) plane intact. The perturbations h^i,{r, t, z) 
are thus classified according to their transformation properties under 0{2). Following 
|4151 1109] . we call them the scalar (hxy), shear {htxj hzx or hty, hzy) and sound {hu, 
htz, hzz, hxx + hyy) channels. (Here we partially fixed the gauge by requiring /i^^^O.) 
The 0{2) symmetry ensures that the equations of motion for perturbations of different 
symmetry channels decouple. Linear combinations of perturbations invariant under the 
(residual) gauge transformations /i^^ — > — V^i^,y — V^/^^ form the gauge-invariant 
variables 

Zi = qHfx + ujHzx , (135) 
Z2 - q^fHtt + 2LuqHtz + w'iJ.z + (2 - /) - OJ^] H , (136) 

Z3 = Hxy (137) 

in the shear, sound and scalar channels, respectively (here Htt — L^htt/r'^ f , Hij = 
L^hij/r^ t), H = L^{hxx + hyy)/2r'^). From the equations of motion satisfied 

by the fluctuations one obtains three independent second-order ODEs for the gauge- 
invariant variables Zi, Z2, Z3. 

The QNM spectra in all three channels share a characteristic feature: an infinite 
sequence of (asymptotically) equidistant frequencies approximated (for g = 0) by a 
simple formula [Ml [Tonl [1(77] 

LOn = 2TiTn (±1 — i) + LjQ , n — !■ cx) . (138) 

Each frequency has a non-trivial dependence on q [99l llOOj . Spectra in the shear 
and sound channels are shown in Fig. (I14p . In addition to the sequence mentioned 
above, they contain the so called hydrodynamic frequencies shown in Fig.[T4]by hollow 
red dots. The hydrodynamic frequencies are remarkable, in that their existence and 
dependence on q are predicted by the hydrodynamics of the dual field theory. For 
example, low-frequency, small momenta fluctuations of the stress-energy tensor T^^^ 
of any d-dimensional theory are characterized by the dispersion relations (see e.g. 

[miiiioiiiiB]) 

' 0{q') , (139) 



UJ = — i- 



p 



UJ = ± Csq - 



P 



+ O , (140) 



where e, P, ry, C are, respectively, energy density, pressure, shear and bulk viscosities, 
and Cs is the speed of sound. For conformal theories, some of the higher order terms 
are also known (see e.g. [416] ). In quantum field theory, the dispersion relations 
()139p , (|140p appear as poles of the retarded correlation functions of the stress-energy 
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tensor. Therefore, according to the holographic dictionary discussed in Section [821 
Eqs. ()139|) . (|140p are precisely the lowest quasinormal frequencies in the spectrum of 
gravitational perturbations. For A/" = 4 SYM in c? = 4 with zero chemical potential 
we have ^ = 0, Cg = l/VS, e + P ~ sT, where s is the (volume) entropy density, and 
the dispersion relations (|139p . (|140p become 

-i-S^q' + 0{q*) , (141) 

On the other hand, the lowest QNMs of the fluctuations Zi, Z2 of the dual black 
brane background (|125p can be computed analytically [4151 1413| : 

0,= -i^<i' + 0(<,') , (143) 

Comparing Eqs. pH]) and (fTIHl) . (fH^ and (HHj), one finds that 1) the real part of 
the mode (jl44p is correctly predicted by hydrodynamics (yet another piece of evidence 
in favor of the AdS/CFT conjecture) and 2) assuming the validity of AdS/CFT, the 
ratio rj/s is equal to l/47r in TV = 4 SYM theory (in the limit of infinite coupling and 
infinite Nc, where the dual gravity description is valid). 
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Figure 14. Quasinormal spectrum of black three-brane gravitational fluctuations 
in the "shear" (left) and "sound" (right) channels, shown in the plane of 
complex frequency XO = u)/2TrT, for fixed spatial momentum q = q/2-KT = 1. 
Hydrodynamic frequencies are marked with hollow red dots (adapted from | 109| ). 

This simple example illustrates a general method of extracting physical quantities 
of strongly coupled quantum field theories from the QNM spectra of their gravity duals. 
Bulk fluctuations corresponding to operators of conserved currents in the dual field 
theory in Minkowski spacctime are guaranteed to have hydrodynamic frequencies in 
their QNM spectra. All such frequencies have the generic property that cj for 
g ^ 0, characteristic of long wavelength, small frequency fluctuations in flat space. 
For M — A SYM, such operators are (in addition to T^^) the i?-current and the 
supercurrent. Their corresponding bulk fluctuations are, respectively, the C/(l) and 
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the Rarita-Schwinger fluctuations in the background (|125p . The hydrodynamic QNMs 
of these fluctuations were first computed in Refs. [Mil HOOl HTT]. The full QNM 
spectrum of electromagnetic/gravitational fluctuations in the background p25p was 
computed in Refs. |1001 1109] , The full QNM spectrum of the Rarita-Schwinger fleld 
has not yet been determined. Note that for theories in a flnite volume (dual to black 
holes rather than branes) hydrodynamic QNMs, strictly speaking, do not exist as the 
momentum q is discrete. However, for the large-radius asymptotically AdS BHs, the 
emergence of the hydrodynamic behavior in the limit r+ /L — s- oo is easily detected 
(see e.g. Eq. ([TMl) ). 

The interpretation of non-hydrodynamic QNMs is obscured by the fact that 
very little is known about thermal correlation functions at strong coupling. Typical 
singularities of such correlators at weak coupling appear to be cuts rather than poles 
(this problem is discussed in Ref. ^418J). Poles associated with the QNM sequence 
(|138p . for example, cannot be interpreted as quasiparticles of A/" — 4 SYM, since their 
imaginary part is large. 

For backgrounds more complicated than the one of the black brane (|125p . 
the gauge-invariant variables Zi may involve flelds other than metric fluctuations. 
Generically, these variables form a system of coupled ODEs (simple examples can be 
found in pTOl [TTTl [4T(1] ) . 

A separate problem is the computation of the corrections to QNM spectra coming 
from higher-derivative terms in the relevant supergravity actions. In the dual field 
theory such corrections correspond, in particular, to coupling constant corrections 
to transport coefficients. For instance, higher-derivative corrections to the type IIB 
supergravity action result in 't Hooft coupling corrections to the shear viscosity- 
entropy density ratio in A/" = 4 SYM gH 1^ HH] 




A>1, (145) 



where C(3) ~ 1.202 is Apery's constant. Corrections to other transport coefficients 
have also been considered l423 l [4241 1425] . 

8.4- Universality of the shear mode and other developments 

The existence of "hydrodynamic" QNMs is a generic feature of black branes: these 
QNMs must appear in the spectra of fluctuations dual to conserved currents in 
translationally invariant backgrounds. In some cases, the lowest QNMs can be 
computed analytically even for rather general metrics. For example, for a black p- 
brane metric of the form 

p 

ds^ = gtt{r)dt^ + grr{r)dr'^ + gxxir) Y.{dx^)\ (146) 

1=1 

with all other background flelds vanishing, one flnds that the component of an electric 
field fluctuation parallel to the brane possesses a "hydrodynamic" frequency 

uj = -iD + O {q^) , (147) 

where the coefficient D is given in terms of the metric components: 



^ ^ ^-g{rH) f ^^-gu{r)9rr{r) ^^^^^ 
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and is the position of the horizon. If there exists a quantum field theory 
(QFT) dual to this background, the electromagnetic fluctuation naturally couples 
to a conserved U{\) current in the theory, and D is interpreted as a diffusion constant 
of a corresponding U{1) charge. 

Similarly, the lowest gravitational QNM in the shear channel l|135[) for the metric 
(|146|) can be computed analytically [114j. The result is an expression similar to 
Eqs. p47p - (|148p . but it can be further simplified using either Buchel-Liu's theorem 
|426j or the alternative proof given in [7]. Surprisingly, it turns out that this mode is 
universal: its frequency is given by Eq. p43p for any background with a metric of the 
form (|146p . (It is important to note that the shear mode fluctuation decouples from 
the fluctuations of other fields [427].) 

The universality of the gravitational quasinormal shear mode has an important 
consequence for a dual QFT: the ratio of shear viscosity to entropy density has a 
universal value rj/s — I/A-k for all theories with gravity duals, in the limit where the 
gravity dual description is valid (e.g., in the limit of infinite Nc and infinite 't Hooft 
coupling for AT = 4 SU{Nc) SYM in four dimensions). Earlier proofs of the shear 
viscosity universality not involving QNMs were given in [4271 1428j . Charged black 
brane backgrounds were considered in [HOl |429l |430l |43T1 |432l HSU |434l HSSl [424] . 

The viscosity bound conjecture |436j prompted further investigation of the black 
brane QNM spectrum in higher-derivative gravity [1371 iSHl SSS 11101 SHI SH SSI 
144411445] . The conjecture states that the ratio rj/s is bounded from below by l/47r in all 
physical systems. In the language of QNMs, this would mean that a correction to the 
g^-coefficient in Eq. (|143p coming from higher-derivative terms in the action describing 
gravity duals of such systems is always positive, as in Eq. (|145p . In some gravity models 
with higher-derivative terms this appears not to be the case [4371 14381 14391 14401 1425| . 
It would be highly desirable to have a comprehensive understanding of the influence 
of higher-derivative terms on the shear mode and other QNMs. 

QNMs of Dp-branes and more complicated backgrounds in the holographic 
context were computed in |446| I447|. 11121 I448j . Time-dependent backgrounds 
especially relevant for modeling the behavior of the quark-gluon plasma in heavy ion 
collisions were investigated in [4491 13161 1163] . QNMs in models with holographic 
mesons, Sakai-Sugimoto model and other QCD-like models were studied in [4501 14511 
14521 14531 14M 14551 14561 14571 [458] . 

The gauge-gravity duality is primarily used to investigate strongly coupled gauge 
theories with the help of a dual classical gravity theory. However, at least in principle, 
the holographic dictionary can be used to explore aspects of quantum gravity with the 
help of a weakly coupled gauge theory. QNMs in this and similar contexts of duality 
have been studied in [1051 (Ml 13311 [^ . 

Recently, holographic methods were extended to include gravity backgrounds dual 
to non-relativistic theories as well as systems with spontaneous symmetry breaking. 
Exploration of QNM spectra in these models is an active area of research [371 12H] • 

9. Quasinormal modes of astrophysical black holes 

In the context of the Einstein-Maxwell equations, the most general solution describing 
stationary axisymmetric BHs is the Kerr- Newman metric [39] . For astrophysical BHs 
the electric charge Q is likely to be negligible, being shorted out by the surrounding 
plasma [4591 1460] . Astrophysical BHs are effectively the simplest of all macroscopical 
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objects, characterized only by their mass M and angular momentum parameter a. As 
a consequence, their structure and their oscillation spectra are remarkably simple. 

The complex frequency of a QNM yields two observables: the actual oscillation 
frequency and the damping time of the oscillation. For each given mode, these 
observables depend only on M and a. Therefore a measurement of the frequency 
and damping time of a QNM can be used to infer the mass and angular momentum of 
the BH with potentially high accuracy [9l [79l [TOl [11] . Since the whole QNM spectrum 
depends solely on M and a, the measurement of two or more QNM frequencies 
provides a stringent observational test of the no-hair theorem of general relativity 
|461[ I462[ [TTl 1463] . The prospects for detecting the signature of BH osciUations in 
gravitational waves are the main topic of this section. 

There is strong and growing observational evidence for the existence of at least 
two different classes of astrophysical BHs. Solar-mass BHs with A/ ^ 5 — 20 Mq 
are usually found in X-ray binaries [464] . and SMBHs with M ~ 10^ — 10^'^ Mq are 
believed to harbor most Active Galactic Nuclei (AGN) '465]. At the moment there is 
only tentative evidence for intermediate-mass BHs (IMBHs) of mass M ~ 10^ — 10^ Mq 
|466| I467|. I468j . Ringdown detectability from these different classes of astrophysical 
BHs depends on several factors, the first of which is the sensitivity of any given 
gravitational wave detector. Earth-based detectors, such as LIGO and Virgo, have 
an optimal sensitivity band corresponding to stellar-mass BHs and IMBHs, while the 
planned space-based detector LISA is most sensitive to ringdowns from high-mass 
IMBHs and SMBHs (see Section l9T] below) . Another important factor is the extent 
to which QNMs are excited in astrophysical settings. The most promising scenarios 
to excite ringdown to a detectable level are reviewed in Section 19.21 Whatever the 
source of the excitation, the frequency and amplitude of ringdown waves depend on 
the mass and spin of the BH. The evidence for astrophysical BHs and the present 
understanding of their mass and spin distributions are reviewed in Section 19.31 
Measuring ringdown waves will allow us to extract interesting information, ranging 
from accurate measurements of the mass and spin of BHs to tests of the no-hair 
theorem of general relativity. These applications of gravitational wave detection are 
reviewed in Sections 19.6119.81 



9.1. Physical parameters affecting ringdown detectability 



Present and planned gravitational wave detectors are located at large distance from 
astrophysical BHs. Therefore, for all practical purposes, a QNM as seen by a 
detector is well approximated by the asymptotic behavior of the wave equation at 
infinity, Eq. (|5D|) . We can express the waveform measured at the detector as a linear 
superposition of the gauge-invariant polarization amplitudes h^, where, for a given 
mode m, n). 



M 

—Re 

r 

hx — — Im 
r 



4+ p 



(149) 
(150) 



Here and 0/^„ are the (real) amplitude and phase of the wave, and S';m„(t, /?) 

denotes spin- weighted spheroidal harmonics of spin-weight —2 [10]. The angles {l, (3) 
are adapted to the source, so that the 2— axis is aligned with the spin of the BH. 
Interferometric detectors are sensitive to the effective strain 

h = h+F+{9s, cl)s,i^s) + hxFx {Os, cl)s,i^s) , (151) 
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[AtS + i-^Lnf] ■ (153) 



where f+,x are pattern functions that depend on the orientation of the detector and 
the direction of the source (specified by the polar angles 6s, 4>s) and on a polarization 
angle ipS |153| I469j . A crucial quantity for gravitational wave detection is the signal- 
to-noise ratio (SNR) p, defined as 

where h{f) is the Fourier transform of the waveform, and Sh{f) is the noise spectral 
density of the detector [79j. In discussing the SNR we will usually average over 
source direction, detector and BH orientations, making use of the sky averages: 
{FD = = 1/5, (F+Fy,) = 0, and {\Simn\^) = l/Air. An analysis taking into 

account different sky locations and orientations of the source probably requires Monte 
Carlo methods. At the moment, such an analysis is still lacking. 

Our chances of detecting and measuring ringdown waves are mainly determined 
by the BH's mass M, by the spin parameter a, and by the ringdown efficiency eid. 
The latter quantity is defined as the fraction of the total mass-energy of the system 
radiated in ringdown waves, and it is well approximated by |1Q[ I89j 

QlmnMuJlmn 

These three parameters (mass, spin and efficiency) affect the detectability and 
measurability of the signal in different ways. 

The BH mass sets the frequency scale and damping time of the emitted radiation. 
For a Schwarzschild BH, the fundamental QNM with / = 2 (that dominates the 
radiation in most cases |89j ) has frequency / and damping time r given by 

/ = 1.207 • 10-^(10^MQ/Af) Hz , (154) 
T = 55.37(M/1O^M0) s . (155) 

Earth-based detectors are limited at low frequency by a seismic cutoff fs (a 
plausible estimate for second-generation detectors being fs ^ 10 Hz for the Einstein 
Gravitational Wave Telescope, and fs ^ 20 Hz for Advanced LIGO). Therefore 
they can detect the fundamental I = m ^ 2 QNM as long as the BH mass 
M < 1.2 X lO''(Hz//,)M0 if the BH is non rotating, and M < 2.7 x 10'^(Hz//,)Mq 
if the BH is rotating near the Kerr limit (see Table I in Ref. [II]). LISA is limited at 
high masses (low frequencies) by acceleration noise, and at low masses by the condition 
that the damping time r should be longer than the light-travel time TJight — 16.7 s 
corresponding to the planned armlength (i ~ 5 • 10^ m). Thus LISA can detect 
ringdown waves from BHs in the range 10^ Mq < M < lO^M© [TU]. To summarize: 
Earth-based detectors are sensitive to the ringdown of stellar-mass BHs and of 
relatively low-mass IMBHs, and LISA can observe mergers of IMBHs and SMBHs 
throughout the whole universe. There is a chance that the ongoing ringdown searches 
in data from Earth-based gravitational wave detectors [12, 470, 471, 472. 473^ .4741 113] 
may provide the first incontrovertible evidence of the existence of IMBHs. 

As discussed in Section 15.31 the BH spin (for a given BH mass) determines all 
frequencies of the Kerr QNM spectrum. For QNMs with m > the quality factor 
increases with spin (see Fig. [9]) . Since the detectability of a gravitational wave signal by 
matched filtering scales with the square root of the number of cycles, highly spinning 
BHs could be the best candidates for detection [5S]. However, exciting QNMs of fast- 
rotating BHs seems to be harder, as the excitation factors tends to zero as a/M 1 
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(see [70 l 1167] and Section [3T2|) . In hindsight, this is not surprising: the build-up of 
energy in a long-lived resonant mode usually takes place on a timescale similar to the 
eventual mode damping, so it should be difficult to excite a QNM with characteristic 
damping several times longer than the dynamical timescale of the excitation process 
|263| . Numerical simulations of the merger of comparable-mass BH binaries suggest 
that QNM excitation is mildly dependent on the initial spin of the components (see 
[475] and Sect ion l9.2p . but further investigation is required to clarify this issue. 

There are different ways of quantifying the excitation of QNMs by generic initial 
data [77] 11731 11661 [31 1167| . One can operationally define a ringdown efficiency e-^d, 
which is directly related to the gravitational wave amplitude as illustrated by Eq. (|153p , 
and hence to the SNR of the signal for a given detector [89l ]T0] . Different empirical 
notions for the "ringdown starting time" (which is intrinsically ill-defined |169] ) yield 
very different ringdown efficiencies [476] [44] 1477] . From the point of view of detection, 
a suitable definition of the ringdown starting time is the one proposed by NoUert, 
using an "energy-maximized orthogonal projection" of a given numerical waveform 
onto QNMs O ]TT]. The relative excitation of different modes is even harder to 
determine than the overall ringdown efficiency, but it is particularly relevant for tests 
of the no-hair theorem using ringdown waves |4611 fTO], [TT] . Ref . [28] discusses this 
issue in a general context, and Berti et al. discuss this issue in a general context [28], 
and give preliminary estimates of the excitation of different multipolcs in binary BH 
mergers [4411475] . 

In Section 19.21 we review promising astrophysical scenarios that could produce 
detectable ringdown signals (i.e., large efficiencies): accretion, stellar collapse leading 
to BH formation and compact object mergers. In Section r9.3l we summarize the current 
theoretical and experimental understanding of BH masses and spins. 

9.2. Excitation of black hole ringdown in astrophysical settings 

In principle, most dynamical processes involving BHs excite QNMs to some degree. 
For the purpose of gravitational wave detection from astrophysical BHs, the question is 
not whether QNMs are excited, but whether they are excited to a detectable level. BH 
QNMs can be excited in a variety of astrophysical settings, such as accretion, collapse 
and compact binary mergers. As we will see, analytical estimates and numerical 
calculations show that the most promising source of detectable ringdown waves is the 
merger of two compact objects leading to BH formation. 

Ringdown excitation by accretion An early study highlighting the importance of 
QNM ringing is the classic analysis of the gravitational radiation emitted by particles 
falling radially into a Schwarzschild BH ^43]. Unlike stellar oscillation modes (which 
play an important role in the orbital dynamics of compact binaries [4781 14791 14801 |48T] ) 
BH QNMs are hard to excite by matter orbiting around the BH. The reason is that 
weakly damped QNMs are associated with unstable geodesies at the light ring (see 
[208] and references therein), and for Kerr BHs the light-ring frequency is always 
larger than the frequency of the innermost stable circular orbit (ISCO), as illustrated 
in the left panel of Fig. 1151 Higher overtones may have lower frequencies, but they are 
harder to excite because their quality factor is too small. According to this intuitive 
description, QNM excitation requires the accreting matter to cross the light ring. 
Indeed, QNMs are always excited by particles falling along generic geodesies into Kerr 
BHs (see [35] and Appendix C of [33] for comprehensive lists of references) . Lumps of 
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Figure 15. Left: ISCO frequency, light-ring frequency and radiation efficiency 
for Kerr BHs (positive values of a/M refer to corotating orbits, negative to 
counterrotating orbits). Right: bounds on the final spin expected to result from 
a binary BH merger as a function of the binary's mass ratio q (from Ref. [482) ). 

matter accreting onto a BH at appropriate rates could potentially excite QNMs to a 
detectable amplitude. Unfortunately, even the most optimistic estimates suggest that 
the wave amplitude is too small [266,. For generic accretion flows, simple analytical 
arguments [28j show that destructive interference reduces the ringdown amplitude 
P51 14831 1484j . This conclusion is confirmed by numerical simulations of neutron star 
(NS)-BH mergers: if the NS is tidally disrupted well before merger, accretion of the 
NS material onto the BH proceeds incoherently, and the ringdown signal is replaced by 
an abrupt cutoff in Fourier space at the tidal disruption frequency (see Section [9. 2 1 and 
Ref. [485j ). If it occurs in nature, a highly intermittent hyperaccretion scenario |486j 
may excite QNM ringing in a nearly resonant fashion for large spin parameters, and 
potentially lead to ringdown radiation strong enough to be detectable by Advanced 
LIGO. 

Collapse to a black hole There are several excellent reviews on gravitational wave 
signatures from supernova core collapse and on the underlying physics [23l 14871 14881 
14891 I490j . When the core of a massive star collapses, it produces a nonspherical 
protoneutron star. Depending on details of the supernova explosion, the protoneutron 
star may collapse to a BH, emitting a burst of gravitational waves due to the rapidly 
shriking mass-quadrupole moment of the protoneutron star and to the QNM ringing 
of the nascent BH. 

Perturbative calculations of gravitational wave emission from rotating gravita- 
tional collapse to a BH were first carried out in the 1970s by Cunningham, Price and 
Moncrief plElliS], improved upon by Seidel and collaborators [4911 . l492 t l493 t l494] 
and more recently by Harada, Iguchi and Shibata [495] . These studies suggest that 
gravitational waves are mainly generated in the region where the Zerilli potential 
is large, and that the signal is usually dominated by QNM ringing of the finally 
formed BH (see |495j for exceptions). Simplified simulations based on a free-fall 
(Oppenheimer-Snyder) collapse model yield a small energy output, with a typical 
core collapse radiating up to ~ 10~^Af in gravitational waves. Most of the radiation 
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is quadrupolar (1 — 2), radiation in Z = 3 being typically two to three orders of mag- 
nitude smaller (see Fig. 9 in [H^)- Recent perturbative studies suggest that magnetic 
fields could increase the energy output by several orders of magnitude |4961 1497j . 

The first numerical simulation of collapse in two dimensions was carried out by 
Nakamura [498] , but numerical problems prevented gravitational radiation extraction. 
For a long time the seminal 1985 work of Stark and Piran [71] has been the only 
nonperturbative, axisymmetric calculation of gravitational wave emission from stellar 
collapse. The waveform resembles that emitted by a point particle falling into a 
BH, but with a reduced amplitude. The total energy emitted increases with the 
rotation rate, ranging from - 10"®M for a/M = to ~ 7 x 10"^M as a/M 1. 
Rotational effects halt the collapse for some critical value of (a/M)crit which is very 
close to unity, and depends on the artificial pressure reduction used to trigger the 
collapse. Stark and Piran find that the energy emitted E/M ~ 1.4 x 10^'^(a/M)'* 
for < a/M < (a/M)cnt, and for larger spins it saturates to a maximum value 
Cmax ^ 10^"*. Note however that the maximum energy radiated is very sensitive 
to the amount of artificial pressure reduction used to trigger the collapse: the 99% 
pressure reduction used in the simulations of Stark and Piran essentially produced a 
free-fall collapse, presumably overestimating the radiation efficiency. 

This calculation has recently been improved using a three-dimensional code 
|4991 15001 [7^ 1501] . Ref. [500 choose as the initial configuration the most rapidly 
rotating, dynamically unstable model described by a polytropic equation of state 
(EOS) with r = 2 and K = 100, having a dimensionless rotation rate ~ 0.54, 
and trigger collapse by reducing the pressure by < 2%. The polarization is 
essentially a superposition of modes with I — 2 and I — 4, and the "x" polarization 
is a superposition of modes with I = 3 and 1 — 5. The energy lost to gravitational 
waves according to these simulations is at most ~ 1.45 x 1O~^(M/M0), two orders 
of magnitude smaller than the estimate by Stark and Piran for the same value of the 
angular momentum, but larger than the energy losses found in recent calculations of 
rotating stellar core collapse to protoneutron stars |502j . Ref. |501j confirms the basic 
scaling E/M Ki {a/M)'^ for a/M < 0.54 (the largest rotation rates yielding equihbrium 
models in uniform rotation), but at variance with Stark and Piran's simulations it 
finds that the efficiency has a local maximum for large rotation rates. If the collapse 
is triggered only by pressure depletion, the overall efficiency for uniformly rotating 
models is quite low {E/M « 10"'' — 10~^), the mass quadrupole does not change 
very rapidly, and higher multipoles are not significantly excited. However the total 
energy radiated can increase by up to two orders of magnitude if velocity perturbations 
are present in the collapsing star. A rapidly rotating polytropic star at 10 kpc can 
produce maximum characteristic amplitudes /ic ~ 5 x 10"^^ (M/Mq) at characteristic 
frequencies ranging between fc « 500 Hz and fc ~ 900 Hz, and could be detectable 
with SNRs as large as about 30 by Earth-based interferometers. 

In conclusion, the ringdown efficiency is very sensitive to poorly known details of 
the mechanism triggering the collapse. Simulations with more realistic microphysics 
are required to study the complex signal preceding the ringdown phase, as well as 
details of QNM excitation by matter accreting onto the newly formed BH. 

Mergers of compact objects leading to black hole formation Numerical relativity 
simulations of compact binary mergers made enormous progress since 1999, when 
the two classic reviews on QNMs |3J [51 were written. This section is an attempt to 
summarize aspects of this progress of interest for the detection of ringdown waves from 
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astrophysical BHs. 

Since NSs cannot have masses larger than about SMq, Eq. (f84|) implies that 
NS-NS mergers are potentially relevant only for the detection of BH ringdowns by 
Earth-based interferometers, such as LIGO and Virgo. On the other hand, NS-BH 
and BH-BH mergers are plausible targets for ringdown detection by both Earth-based 
and future space-based detectors, so we will discuss them in more detail. 

(1) NS-NS mergers - In 1999 Shibata and collaborators carried out the first successful 
equal- mass NS-NS merger simulation for a polytropic index F = 2 }503|, [93] . The 
first reasonably accurate calculation of gravitational waveforms was possible a few 
years later, in 2002 [504j . Unequal-mass binaries with a F = 2 polytropic EOS were 
studied in [505] . where it was found that a BH forms when the total rest-mass of the 
system is larger than ~ 1.7 times the maximum allowed rest-mass of spherical NSs, 
irrespective of mass ratio (which however affects significantly the waveforms and the 
mass of the disk forming around the newly born BH). More realistic EOSs and larger 
parameter spaces were considered in Refs. |5061 1507] . Ref. |507] uses stiff EOSs and 
models binary NSs of ADM mass M > 2.6Mq. For all mass ratios 0.65 < 9 < 1, a 
BH forms whenever the mass M > Mth- The threshold value Mth — 1.3 — 1.35Afniax 
(where M„iax is the maximum mass allowed by the given EOS for cold, spherical NSs) 
depends on the EOS. If M < Mth the merger results in a hypermassive NS of large 
ellipticity, emitting quasiperiodic gravitational waves at frequencies ~ 3 — 4 kHz for 
< 100 ms. After this phase the NS may or may not collapse to a BH. For total binary 
masses in the range M ^ 2.7 — 2.9Mq, the I — m — 2 ringdown frequency emitted 
in the collapse to a BH is « 6.5 — 7 kHz, with amplitude ^ 10~^^ at a distance of 
50 Mpc. Therefore BH ringdown from NS-NS mergers is unlikely to be detected: the 
amplitude is too low and (most importantly) the frequency is too high for present and 
planned Earth-based gravitational wave detectors. It is interesting to note that NS-NS 
merger simulations typically lead to final BH spins a/M « 0.8 j507[ 1508] . not very 
different from the value a/M ~ 0.69 predicted by equal-mass BH merger simulations. 
Various groups have recently carried out NS-NS merger simulations with magnetic 
fields, finding that aligned poloidal fields can delay the merger and strongly affect the 
gravitational wave signal [5091 15101 1511] . 

(2) NS-BH mergers - NS-BH binaries are potentially among the most interesting 
BH ringdown sources. From Eq. (|154|) . the QNM frequency for non-spinning BHs 
with M > IOM0 is within the sensitivity window of LIGO and Virgo. Theory and 
observations [5121 1464[ 1513] suggest that binaries containing BHs with masses in this 
range should be common. Furthermore, high mass-ratio systems do not tidally strip 
the NS, producing a merger with a clean ringdown signal. The simulations described in 
the rest of this section suggest that NS-BH mergers are promising sources of ringdown 
waves and (when the BH is spinning) excellent candidates as central engines for 
gamma-ray bursts (GRBs). 

NS-BH merger simulations became possible once the moving puncture approach 
[15l[T6] proved successful in evolving BH-BH binaries. Shibata and Uryu [514il515j first 
evolved two binaries consisting of (1) a non-rotating BH of mass « 3.2Mq or ^ 4Mq, 
respectively, and (2) a NS of rest mass « 1.4M0, modeled by a F-law EOS with F = 2. 
A larger set of simulations with non-spinning BHs of masses M ~ 3.3 — 4.6M0 and 
NSs of mass M ~ I.4M0 was performed in Ref. [516] . The NS is tidally disrupted, and 
the system results in a BH with spin a/M 0.5 — 0.6. Non-spinning NS-BH mergers 
can only produce massive tori and fuel short-hard GRBs if the NS has compactness 
M/R < 0.145. Furthermore, when tidal disruption occurs the QNM amplitude quickly 
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decreases, because the incoherent accretion of material is ineffective at resonantly 
exciting the mode (see Fig. 6 of Ref. |516j ). 

The first code capable of handling NS-NS, NS-BH and BH-BH binary evolutions 
comprising more than 4 orbits was recently described in Ref. [508] . Gravitational 
waveforms from NS-BH mergers were presented in Ref. [485) and classified into three 
classes as follows: 1) For small mass ratio q and small NS compactness (e.g., if 9 ^ 3 
for M/R — 0.145 and F = 2) tidal disruption occurs outside the ISCO, and there 
is no QNM excitation; 2) For some systems, mass shedding occurs before the binary 
reaches the ISCO. Most of the NS is swallowed by the BH before tidal disruption is 
completed and QNMs are excited, but only to a low amplitude; 3) If tidal effects do 
not play an important role the waveforms shows significant QNM excitation, as they 
always do for BH-BH mergers. The latter class of NS-BH mergers is clearly the most 
promising for BH ringdown detection. Whenever ringdown is not significantly excited, 
kicks are also suppressed: this confirms the crucial role played by the merger/ringdown 
phase in determining the magnitude of the kick resulting from compact binary mergers 
(see Section |9^ below) . The results of Ref. [485] are in good agreement with NS-BH 
codes developed by other groups [517[ I518[ I519j . Ref. [517j deals with systems where 
the NS is irrotational, the BH is non-rotating and the mass ratio q = 1/3. These 
simulations lead to the formation of BHs with a/M k, 0.5 — 0.8. Most of the NS 
material is prompty accreted, and no more than 3% of the NS mass is ejected into a 
gravitationally bound disk. This disk mass, while larger than the typical values found 
in Ref. [485], is probably not enough to trigger short-hard GRBs. 

Results from non-spinning NS-BH mergers show that the disk mass is typically 
too low to fuel GRBs. For example, Fig. 7 of Ref. [485j shows that, if the BH is 
non-spinning, the formation of a disk requires a "fat" NS with radius i? > 14 km. 
The few simulations presently allowing for (aligned or anti-aligned) BH spins have a 
drastically different outcome, indicating that spin plays a crucial role in fueling GRBs 
[518] . In these simulations the number of orbits before merger increases as the spin is 
varied between a/M = —0.5 (anti-aligned), 0.0 and 0.75 (for fixed mass ratio q ~ 1/3). 
In the latter case the final BH spin is a/M ~ 0.88, and the tidal disruption of the 
NS leads to the formation of a massive disk of about 0.2Mq, potentially capable 
of driving GRBs. The production of GRBs by NS-BH mergers and the possibility 
to detect the final spin orientation using observations of the merger/ringdown phase 
are promising areas of research for multi-messenger astronomy using a network of 
ground-based detectors in conjunction with traditional electromagnetic observations 
(see e.g. [5^ 15^ 15^ 15^ 15^ 1. 

(3) BH-BH mergers - After 40 years of developments in numerical relativity, recent 
breakthroughs [Mj [151 [16] finally allowed simulations of the merger and ringdown 
of BH binaries. Extensive collaborations to use numerical merger waveforms in 
gravitational wave searches have just started [525] . A discussion of the accuracy 
of numerical simulations and of their rapid progress in the last few years would take 
us too far (see Refs. [Mj 1526) for reviews). From the point of view of this review, 
the main result of these merger simulations is that QNM ringing is observed in all 
binary BH merger simulations. This is quite unlike NS-BH and NS-NS mergers, where 
tidal effects can sometimes suppress ringdown excitation because of the incoherent 
accretion of material onto the newly formed BH. Another important difference is that 
the gravitational waveform from a BH-BH merger depends on the total mass of the 
system via a trivial rescaling, so BH-BH systems are interesting for ringdown detection 
by both Earth-based and space-based detectors. 



Quasinormal modes of black holes and black branes 



67 



The mass-ratio dependence of the ringdown efficiency can be estimated by simple 
arguments [89l |44j . The quadrupole moment of a body of mass M with a " bump" of 
mass n M is Q ^ fiM'^ . The oscillation frequency of the system / ~ l/M; hence, 
the radiated power dE/dt ~ ((PQ/dt^) ^ [f^QY ^ [i^ jM"^. For a binary with mass 
ratio ^/A/ the inspiral lasts ^ {M/ ji) cycles times the orbital time scale T ^ M, 
so the total energy loss during the inspiral is i?insp ~ (M//i)(/i^/Af^)M ^ ^. By 
contrast, a typical ringdown waveform lasts only a few cycles, so the ringdown energy 
loss -Eringdown ~ M {fj? / M"^) ^ [i} j M (comparc with the classic result for infalling 
particles in Ref. f43|) and -EiingdownZ-Einsp ^ [ijM: ringdown is negligible with respect 
to inspiral for extreme-mass ratio binaries. If we naively extrapolate these estimates 
to bodies of comparable masses (interpreting M as the total mass of the binary and 
/J, — > mim2/M as the reduced mass) we find that -Eringdown/Af ~ 77^, where r] = ^i/M 
is the so-called symmetric mass ratio (77 1/4 in the comparable-mass limit). 

Physical arguments to estimate the prefactors suggest that the merger/ringdown 
waveform should actually dominate over inspiral for binaries with mass ratio q > 1/10 
[891 I10| - and numerical simulations of quasicircular inspirals of comparable- mass 
mergers have borne out this expectation. For non-spinning binary BH mergers, the 
fraction of energy radiated (Af — Affin)/A'/ (where M denotes the total mass of two 
BHs in isolation and Man is the mass of the final BH), as well as the final spin a/A/, 
have been extensively studied [4771 HI ISMl MM- Ref. [527] fits data from the 
simulations by the Goddard group by a relation of the form 



This result is consistent with the fitting formula given in ^44j using a different 
normalization. Frequencies and damping times of different multipolar components can 
be estimated using either a standard least-squares algorithm [477] or Prony methods, 
which are in many ways optimal to estimate the parameters of damped exponentials 
in noise [529j . By monitoring the frequencies and damping times after merger we can 
monitor inaccuracies in the higher multipolar components of numerical simulations, 
and possibly explore non-linear effects (see e.g. Section IV in Ref. 44J). 

Quantifying the fraction of energy radiated in ringdown is inherently ambiguous. 
The reason for this difficulty can essentially be traced back to the non-completeness of 
QNMs [169|, [2|[3ll44j. An operational viewpoint to isolate the ringdown contribution is 
given by Nollert's energy-maximized orthogonal projection (EMOP) criterion [169|, l44j. 
The idea is to determine the starting time of a ringdown waveform by assuming 
that the frequency of the ringdown waveform is known, and performing matched- 
filtering (in white noise) of the numerical waveform, using a damped sinusoid as the 
"detection template" . Since ringdown is essentially monochromatic, this should give 
a reasonably good, frequency-independent and detector-independent estimate of the 
fraction of energy that we can expect to detect by a ringdown search. For more details 
on ringdown search techniques, see Section [9. 41 b clow . 

An application of the EMOP criterion to numerical simulations shows that 
ringdown typically contributes ~ 42% of the energy radiated by the last two cycles 
of BH-BH mergers with mass ratio g > 1/4 (see Tables VI and VII of Ref. p4]). 
By combining Prony methods for estimating frequencies and the EMOP approach to 
estimate the ringdown starting time, we find the following estimates for the energy 
emitted in ringdown as a result of the merger of non-spinning, quasi-circular BH 



Mfi„/Af = 1 + (78/9 - l)7j - 0.498(±0.027)7y2 
a/Mfin = VUt] - 2.900(±0.065)772 . 



(156) 
(157) 
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binaries [4j 

^'rngdown ^ ^ 0.104^^^i^i]^ . (158) 

For non-spinning, equal-mass binaries, the merger/ringdown signal as the 
system's orbital angular momentum is reduced (so a quasicircular merger slowly turns 
into an head-on collision) has been studied in Refs. [53011531] . In the head-on limit one 
gets a radiated energy ^ 0.1%M. In the intermediate regime, where two equal-mass 
BHs merge along orbits with large residual eccentricity, the radiated energy decays 
(roughly) exponentially (see Table I and Fig. 7 in Ref. |530j ) and the final BH spin has 
a local maximum jfin — 0.724zfc0.13 |530l[53T] . For prehminary studies of the ringdown 
efficiency in the merger of quasicircular, spinning binaries, see Refs. |4751 15321 1533| . 

Predictions of the spin of a BH resulting from a merger are very interesting from 
the point of view of ringdown detection. For example, if we measure the masses and 
spins during the inspiral we may be able to predict the final spin and reduce the 
errors in parameter estimation [5341 1535j . For a summary of semianalytic models and 
fitting formulas to predict the final spin from generic mergers, we refer the reader to 
|5361 15371 15381 1539] . Some insight into the general outcome of a spinning merger can 
be obtained by looking at the right panel of Fig. [151 For simplicity, consider three 
different merger scenarios: i) in the isotropic scenario, both BH spins are distributed 
isotropically; ii) in the aligned spin scenario, the individual BH spins in the binary 
are assumed to be aligned (for example, in "wet mergers" the alignment could be 
caused by torques from accreting gas, as suggested in [540| ): iii) in the equatorial 
merger scenario, the smaller BH is supposed to orbit in the equatorial plane of the 
larger hole (e.g. because of Newtonian dynamical friction in a flattened system), but 
the spin orientation of the smaller BH is distributed isotropically. The right panel of 
Fig. [m shows the maximum and minimum spin resulting from a merger in these three 
scenarios. These curves are obtained by (1) fixing some value of the mass ratio g, (2) 
averaging over angles according to the three different assumptions listed above, and 
(3) maximizing or minimizing the final average spin resulting from a merger (where the 
average is computed using the fitting formulas of Ref. |537| ) in the (|ji|, jjjl) plane, 
where |jj| is the spin magnitude of BH i = 1, 2. Not surprisingly, the minimum 
average final spin always corresponds to the case where both BHs are non-spinning 
(dashed black line). The maximum average spin in the three cases is shown by the 
continuous black (isotropic), red (equatorial) and blue (aligned) fines. The dashed 
blue line shows the (modulus of) the minimum spin that could be achieved if we allow 
for antialignment of both spins with respect to the orbital angular momentum (a spin 
flip becomes possible when the mass ratio q « 1/3). The most interesting prediction 
of this plot is the existence of a narrow funnel between the solid black and dashed 
black lines: on average, isotropic major mergers (with q > 0.2 or so) always produce 
a final spin which is very close to the spin resulting from equal-mass non-spinning BH 
binaries, i.e. a/M ~ 0.69. Furthermore, in all three scenarios the most likely spin 
resulting from "major" mergers with q > 0.1 is quite close to |jg„| ~ 0.69. 

A critical assessment of the available predictions on the final spin and on the final 
kick is outside the scope of this review; the interested reader is referred to Ref. [541j . 

9.3. Astrophysical black holes: mass and spin estimates 

From the discussion in the previous section it should be clear that ringdown 
detectability depends crucially on the physical parameters of astrophysical BHs. 
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What are the most promising scenarios leading to the formation of massive stellar- 
mass BHs that would predominantly be seen in merger/ringdown by Earth-based 
interferometers? What are the event rates for mergers to be observed by LISA? Given 
that ringdown detectability scales strongly with mass ratio, what is the most likely 
mass ratio for binaries whose merger is detectable by LIGO and LISA? What is the 
most likely spin magnitude of BHs in binaries, and what are the odds that we can 
unveil connections between merger/ringdown waveforms and the central engine of 
GRBs? In this section we briefly survey theoretical expectations and state-of-the- 
art measurements of the mass and spin distribution of BHs. We summarize some of 
the most relevant information that astronomers have collected by working as busy 
bees over the past thirty years, paying particular attention to the implications for the 
detection of gravitational waves from the merger/ringdown of BH binaries. 

Our focus here is on testing the BH nature of astrophysical objects by 
gravitational wave observations, but there are excellent reviews on measuring BH 
masses, spins and (possibly) providing evidence of an event horizon by "traditional" 
electromagnetic astronomy. Narayan reviews the status of BH astrophysics, focusing 
on observational progress in measuring mass and spin and on (circumstantial) 
observational evidence for the defining property of a BH: the event horizon [8]. Psaltis 
discusses how electromagnetic observations of BHs and neutron stars can be used to 
probe strong-field gravity; in the process he describes various ways of identifying BHs 
and measuring their properties, including continuum spectroscopy, line spectroscopy, 
and attempts at imaging the vicinity of BHs to constrain their angular size |32j . 



Stellar-mass black hole candidates The most accurate mass measurements for stellar- 
mass BH candidates are made via dynamical methods, that is, by looking at how the 
unseen BH affects the orbit of a companion star. Consider a test particle in circular 
orbit around the BH. If the orbit is wide enough for Newtonian physics to apply, then 
the mass M = uj'^r^ = v^r = f^/w, where r is the orbital separation, v is the orbital 
velocity and uj — 2tt/T with T the orbital period (simple modifications can account 
for orbital eccentricity). By measuring any two of f , r and w, we may estimate the 
BH mass M. In the case of BH X-ray binaries it is relatively easy to measure cu and 
the maximum line-of-sight Doppler velocity Kc — v sin t of the companion star. From 
these quantities one can compute the "mass function" 

./,rN Kc Msin^t 

/ M = — = ■ — -r , (159) 

where M, Mc are the masses of the BH candidate and of the companion, respectively. 
The inclination angle i of the orbit can be estimated from the light curve of the binary, 
and sometimes it's even possible to estimate Mc- By combining measurements of uj 
and Kc with estimates of l and Mc, one can in principle determine the masses of both 
binary members. However, to identify BH candidates the essential point is to note that 
the mass function /(M), which depends only on to and provides a strict lower 
hound on M. Since NSs cannot be more massive than about SM© j542j . all X-ray 
binaries for which f{M) > 3Mq should contain a BH. Remillard and McClintock 
[464j review the phenomenology of 20 X-ray binaries with dynamically confirmed 
BHs, presenting a census of BH candidates and a critique of different methods for 
measuring spins. Their Table I provides a list of (lower bounds on) the mass of about 
20 BH candidates. The most massive stellar-mass BH candidate to date is IC 10 
X-1, with a minimum mass M = 23.1 ± 2.1Mq if the companion's mass Mc = 17Mq 
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(M = 32.7 ± 2.6Afo if one trusts an estimate of — SSM© for the companion) 
|513j . This system is particularly interesting for gravitational wave detection, because 
it should become a close double BH binary with coalescence time of ~ 2 — 3 Gyr [543j . 

The relative importance of ringdown with respect to inspiral waves decreases 
for extreme mass ratios. For the stellar-mass binaries of interest for Earth-based 
interferometers, population synthesis codes suggest that q should always be very close 
to unity [512j . A rather speculative kind of source (in the absence of solid evidence 
for IMBHs) consists of the intermediate mass ratio inspiral (IMRI) of a compact 
object, such as a NS or BH, into an IMBH. The relatively low energy content in 
ringdown waves is compensated, in this case, by the fact that the ringdown frequency 
is close to the minimum of the Advanced LIGO noise power spectral density (see 
Appendix B of Ref. |544j ) . Another promising ringdown source for advanced Earth- 
based interferometers (albeit with highly uncertain event rates) are IMBH-IMBH 
inspirals. These binaries, if they are numerous enough to be detectable, present an 
interesting data analysis challenge: the initial inspiral phase could be detected by 
LISA, while the ringdown phase is in the optimal bandwidth for Advanced LIGO, 
that could therefore be used for "follow-up" ringdown searches [545] . 

Supermassive black hole candidates A good review of SMBH observations from a 
historical perspective can be found in Ref. [465 . The first quasar was identified in 
1963 [ni], when the Kerr solution had just been discovered fSU" and its astrophysical 
relevance was unclear. In the intervening years astronomers gathered strong 
observational evidence for the presence of SMBHs in the bulges of nearly all local, 
massive galaxies j546| 15471 1548j . Reliable mass estimates are available for many of these 
systems. The most precise measurement comes from observations of stellar proper 
motion at the center of our own galaxy, indicating the presence of a "dark object" of 
mass M ~ (4.1 ±0.6) x 1O^M0 [54911550] and size smaller than about one astronomical 
unit |551| . A Schwarzschild BH of the given mass has radius R ~ 0.081 astronomical 
units, compatible with the observations. Any distribution of individual objects within 
such a small region (with the possible exceptions of dark matter particles or asteroids, 
which however should be kicked out by three-body interactions with stars) would 
be gravitationally unstable [55211553] . Theoretical alternatives to SMBHs (e.g., boson 
stars and gravastars) have been proposed by various authors, but the formation process 
of these hypothetical objects is unclear, and many of these exotic alternatives can be 
shown to be unstable [554[ 1555] . Another accurate mass measurement comes from 
the motion of the gas disk at the centre of the nearby galaxy NGC 4258 [556], as 
monitored by radio interferometry of the waves emitted from water molecules via 
maser emission. The observations imply the presence of an object of mass 3.5 x IO^Mq 
within ~ 4 X 10^^ m. Other techniques include applications of the virial theorem to the 
velocity dispersion of stars near the galactic center [546] and reverberation mapping to 
obtain more crude estimates for distant, variable AGNs [557] . The reader interested 
in a SMBH mass census can consult Graham's survey [558] . listing 76 galaxies with 
direct SMBH mass measurements and (when available) their host bulge's central 
velocity dispersions. Graham also lists 8 stellar systems that could potentially host 
intermediate-mass BHs. For our purposes, it suffices to note that SMBHs have masses 
in the range M ~ 10^ — IO^Mq, approximately proportional to the mass of the host 
galaxies, M ^ 10^"^ -Mgaiaxy [559] , There is an almost-linear relation between the 
mass of a SMBH and the mass of the galactic bulge hosting it [5461 15471 13i5] . The BH 
mass is also tightly correlated with other properties of the galactic bulge, such as the 
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central stellar velocity dispersion a, the bulge light concentration and the near-infrared 
bulge luminosity [560| I56H I558j . These correlations clearly indicate that SMBHs are 
causally linked to the surrounding galactic environment. The growth of galaxies and 
SMBHs must be intertwined, and observations of the merger and ringdown of SMBHs 
with LISA holds great promise to clarify their formation history. 

Supermassive black hole binary candidates The SMBHs harboring nearby galactic 
cores are expected to grow via a combination of mergers and accretion, and hierarchical 
merger models of galaxy formation predict that binary SMBHs should be common 
in the Universe [5621 15631 1564] . The merger of SMBH binaries is one of the most 
luminous gravitational wave events in the universe, and it is the strongest conceivable 
astrophysical source of ringdown waves. Observational evidence for close SMBH 
binaries started emerging very recently, and it is one of the most exciting frontiers 
of relativistic astrophysics. 

The formation of SMBHs during galaxy mergers is a challenging problem in 
theoretical astrophysics. The general scenario was outlined in a pioneering analysis 
by Begelman, Blandford and Rees [565j (see Ref. [562] for a more recent review). The 
evolution of a SMBH binary can be roughly divided in three phases: i) as the galaxies 
merge, SMBHs sink to the center via dynamical friction; ii) the binary's binding energy 
increases because of gravitational slingshot interactions, i.e. the ejection of stars on 
orbits intersecting the binary (these stars' angular momentum must be in a region of 
phase space called the "loss cone" ) ; iii) if the binary separation becomes small enough, 
gravitational radiation carries away the remaining angular momentum. Notice that 
the gravitational wave coalescence time is shorter for more eccentric binaries j566j . 
so high-eccentricity binaries are slightly more likely to coalesce within a Hubble time 
(see e.g. [ISJJ. The transition from phase ii) to phase iii) is a field of active research, 
that has been referred to as the "final parsec problem" |562| . Since the binary will 
quickly eject all stars through gravitational slingshot interaction, the problem is to find 
some mechanism (such as gas accretion, star-star encounters and triaxial distortions of 
galactic nuclei) to refill the loss cone. It is generally believed, based on both theoretical 
and observational arguments, that efficient coalescence should be the norm [567 . The 
main point here is that only SMBH binaries with separations ^ 1 pc can merge within 
a Hubble time under the sole influence of gravitational radiation. 

The mass ratio distribution is an important variable from the point of view of 
ringdown detection. The impact of different SMBH assembly models on the mass and 
mass ratio distribution of detectable binaries has been discussed by various authors. 
The general consensus is that mass ratios q < 1/10 (and down to g w 10^^ — 10^**) 
should be common [568[ 15691 15701 1571[ I572[ 1573] . This may not be true for mergers 
between BHs in separate dark matter halos, because the smaller halo could get tidally 
stripped and it would not be able to sink efficiently toward the center of the main 
halo [574j . In any case, SMBH mergers will be strong LISA ringdown sources even for 
modest mass ratios. 

Possible observational smoking guns of SMBH binaries (such as X-shaped radio 
galaxies, double-double radio galaxies, helical radio jet patterns, semi-periodic signals 
in lightcurves, double-peaked emission-line profiles and galaxies which lack central 
cusps) are reviewed by Komossa |575j . At the time of her review, the most spectacular 
example of a SMBH binary was the ultraluminous infrared galaxy NGC 6240 [576j . 
containing two active SMBHs separated by a relatively short projected distance 
1 kpc. Since then, more quasars have been identified as promising hosts of SMBH 
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binariefl Even more interestingly, in the last year three systems have been proposed 
to host binary SMBHs at separations smaller than one parsec [5831 I584j . These 
observations could be extremely important for our understanding of SMBH binary 
mergers and their rates |585j . Since they are still controversial and they would take 
us too far from the main topic of this review, we briefly summarize them for the 
interested reader in [Appendix C[ 

Spins: theoretical expectations Estimating BH spins is an exciting observational 
frontier for the next decade of observational astronomy [586|, I587[ [588]. Present 
spin estimates based on electromagnetic observations are all, to some extent, model- 
dependent. Few dependable and accurate measurements are available, and there 
are different opinions on the expected spin magnitude of both stellar-mass BHs and 
SMBHs. The punchline of the theoretical work summarized in this section is as follows: 
1) stellar-mass BHs should typically retain their natal spin, and so they can be used 
to infer information about the mechanism triggering the collapse; 2) SMBH spins 
encode the history of the hole, and particularly the relative importance of mergers 
and accretion in the hierarchical formation process responsible for growing the holes. 
The rest of this section gives arguments in support of these conclusions. 
1) Stellar-mass black holes: Theoretical arguments suggest that stellar-mass BHs in 
binaries retain the spin they had at birth: neither accretion nor angular momentum 
extraction are likely to change significantly their mass or spin. A BH must accrete an 
appreciable fraction of its original mass in order to significantly change its spin. For 
BHs with low-mass companions, even the accretion of the entire companion star will 
only change the spin by a small fraction; for BHs with high-mass companions, even 
Eddington-limited accretion will only grow the BH spin by a small amount before the 
high-mass companion explodes [589| . Therefore the spin of stellar-mass BHs should 
depend mainly on their formation process. 

Detailed studies of spin evolution in compact binaries have been carried out by the 
Northwestern group, focusing mostly on NS-BH binary systems and using progressive 
improvements of the STARTRACK stellar evolution code (see j512j and references 
therein). The evolution prior to the supernova explosion involves mass-transfer phases, 
which are expected to align the spins of both the BH and the NS progenitor, but 
a significant natal kick of the NS at birth is required to form a coalescing NS-BH 
binary |590j . In general, the plane of the post-supernova orbit is tilted with respect to 
the pre-supernova plane, and hence tilted with respect to the BH spin axis by some 
angle i, inducing precession of the binary's orbital plane. Preliminary results suggest 
that precession only marginally impacts the detection of gravitational waves from the 
inspiral waves, but it should be significant for parameter estimation |591[ 1512"! I592j . 
Belczynski et al. [512 consider both BH-BH and NS-BH binaries. For NS-BH binaries 

+ A binary AGN with separation ~ 4.6 kpc lias been claimed in Arp 299 15771 . and a system with 
projected distance ~ 10.5 kpc has been found in the galactic pair ESO 509-IG066 |578) . Evans et 
al. |579l reveal the AGN nature of the companion of the FRII radio source 3C 321. Bianchi et 
al. identify a binary AGN in Mrk 463 with projected separation ^ 3.8 kpc |580l . Finally, using 
multifrequency observations with the Very Long Baseline Array (VLB A), Rodriguez et al. I581II582| 
report the discovery of a SMBH binary in the radio galaxy 04024-379 with a total estimated mass of 
1.5x IO^Mq and a projected orbital separation of just 7.3 pc. This is the smallest orbital separation by 
more than two orders of magnitude, but even for this relatively close binary the emitted gravitational 
waves have frequency ~ 2 X 10"^'' Hz, way too low to be observed by LISA, and a merger time 
~ 10^® yr (much longer than the age of the Universe) i/ gravitational radiation is the only dissipative 
mechanism. 
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they confirm the qualitative predictions of |589j : BHs cannot be significantly spun up 
by accretion in the common envelope phase. For example, only 20% of initially non- 
spinning BHs spin up to a/M > 0.1, and no BHs spin up to a/M > 0.5. The spin-up 
is even smaller for BH-BH binaries, with the highest attainable spins being very close 
to the initial spins of the individual BHs (see Fig. 5 and Fig. 9 in |512| ). Furthermore, 
the kick-induced tilt angle t < 45° for ^ 50% of NS-BH systems. The fraction of 
events that can potentially produce short-hard GRBs, and which therefore is relevant 
for merger/ringdown searches in association with gamma-ray bursts, is significant (of 
order ~ 40%) only if the initial BH spin a/M > 0.6 [H^ . 

These studies suggest that Advanced-LIGO measurements of BH spins in the 
inspiral of binaries containing stellar-mass BHs should be an excellent probe of the 
collapse mechanism that produced the BH in the first place, because the BH essentially 
retains the spin it had at birth. Furthermore, in view of the NS-BH simulations 
discussed in Section Wl^ NS-BH binaries where the BH is rapidly spinning are good 
candidates as central engines of GRBs, and the resulting gravitational wave signal 
may have a significant ringdown component. This could potentially allow a ringdown- 
based measurement of the final BH spin magnitude and direction, with interesting 
implications for coincident electromagnetic/gravitational observations of GRB events. 
2) Supermassive black holes: Since SMBHs are expected to grow by a combination 
of mergers and accretion, their spin will depend on three main ingredients [593j : i) 
the spin of "seed" BHs at birth, that in some sense defines the initial conditions for 
the problem; ii) the spin resulting from a binary BH merger; iii) the maximum spin 
attainable by accretion. We discuss these points, in turn, below. 

i) Natal spins - Little is known observationally about the formation of the first BHs in 
the universe. A popular formation scenario involves the collapse of primordial, massive 
(M ~ 30 — 300 Mq), metal-free Population HI stars at cosmological redshift 2; ~ 20 to 
form primordial BHs with M ~ 10^ Mq, clustering in the cores of massive dark-matter 
halos [594j , but details of the collapse are uncertain |595| I596j . Shibata and Shapiro 
simulated the collapse of uniformly rotating stars supported by radiation pressure and 
spinning at the mass-shedding limit, finding numerically [597j that the final BH spin 
(independently of the progenitor mass) is a/M « 0.75, and supporting this result 
by analytical arguments [598] . Alternative scenarios suggests that BH seeds would 
form at z~12 from low-angular momentum material in protogalactic discs |599| 1600] : 
these seeds would have larger mass M ~ 10^ Mq, but their angular momenta depend 
on the dynamics of the collapsing material. Whether this "initial value problem" is 
relevant for the overall spin distribution of SMBHs is a matter of debate. Merger tree 
simulations where the accretion disk orientation is chosen randomly show that the spin 
distribution does indeed retain memory of the initial conditions |5691 1482] . However, 
a recent model where the BH spin directions are "linked with the galaxies" suggests 
that the spin distribution could be largely independent of the initial conditions |601j . 

ii) Spin from mergers - A pioneering attempt to study massive BH spin evolution 
by repeated mergers, predating the 2005 numerical relativity breakthrough, is due to 
Hughes and Blandford [602j . Extrapolating results from BH perturbation theory they 
found that "minor mergers" (g < 1/10) of a large BH with an isotropic distribution 
of small objects tend to spin down the hole. An implementation of the results shown 
in the right panel of Fig. [TS] within merger tree scenarios confirms these results. 
Furthermore, it shows that SMBH spins cluster around a/M ~ 0.7 if alignment and 
accretion are inefficient, so that the BH spin growth is dominated by mergers [482]. 
However, Volonteri et al. |569|I482] argued that on average accretion should dominate 
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over mergers in determining the spin evolution in hierarchical SMBH formation 
scenarios. 

Hi) Spin from accretion - The details of SMBH growth by accretion are very 
uncertain. It is usually believed that prolonged accretion should lead to large spins 
|569| . However, King et al. [603[ I604[ 16051 1606] suggested that gas accretion may 
occur through a series of chaotically oriented episodes, leading to moderate spins 
a/M ^ 0.1 — 0.3. Since numerical relativity suggests that comparable-mass mergers 
should not be efficient at spinning up BHs, and such comparable-mass mergers are 
expected to be common j540| . a few measurements of spins a/M > 0.9 (such as the 
value of a/M = 0.989^^;^^^ claimed by Brenneman and Reynolds [607] for the Seyfert 
1.2 galaxy MCG-06-30-15) would imply that chaotic accretion is not the norm. 




Inefficient Prolonged Chaotic Inefficient Prolonged Chaotic 



Figure 16. Left panel: spin distributions for merging BHs. Red dashed 
histograms show the spin distribution of inspiralling SMBHs in binaries, while 
black histograms show the spin distribution of SMBHs post-merger. Right panel: 
spin distribution of the whole MBH population in the mass range 10^ Mq < M < 
3 X IO^Mq. For simplicity, we consider on ly the isotropic merger scenario. Left 
to right in each panel, accretion is assumed to be inefficient, prolonged or chaotic, 
respectively (see Ref. 14821 and the text for details). 

Given the present uncertainty on the physical agents responsible for SMBH 
growth, it would be extremely valuable to find observational signatures of 
different formation scenarios. Figure (TH] (adapted from Ref. [482] ) shows that 
electromagnetic spin measurements and gravitational wave measurements from the 
inspiral and ringdown may provide an excellent way of discriminating between different 
mechanisms of BH growth. The spins are assumed to be isotropically distributed 
(but efficient alignment, as suggested in Ref. |540j . would only marginally alter the 
picture: see Ref. [482]) and the seed BHs are assumed to be non-spinning. In the left 
panel, the red dashed histograms show theoretical estimates of the spin distribution 
that could be measured by observing gravitational radiation from inspiralling SMBH 
binaries; the black histograms show the distribution of SMBH spins post-merger, as 
would be measured by observations of ringdown waves. Finally, the histograms in 
the right panel show the distribution of spins of the whole BH population, that can 
be probed by electromagnetic observations of the kind described in [Appendix D[ In 
each panel time runs upward (the histograms correspond to different redshift cuts, as 
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indicated in the inset) and the three "columns", from left to right, correspond to: (1) 
spin growth being determined only by mergers (inefficient accretion), (2) spin growth 
being driven by mergers and prolonged accretion [569], (3) spin growth being driven 
by mergers and chaotic accretion .603, .6041 1605j . The spin distributions are obviously 
very different. According to the chaotic accretion scenario, ringdown measurements 
would never observe BH spins larger than 0.7 or so, and most spins would be very 
low. If prolonged accretion dominates, then most BHs should be rapidly spinning, and 
if accretion is inefficient the spin distribution should have an attractor around j ~ 0.7. 
The spin distribution clearly encodes information on the SMBH merger history. 

It is important to keep in mind that the expected spins of SMBHs may well depend 
strongly on their masses. It has been suggested that SMBHs with M < 2 x IO^Mq 
may grow primarily by disruption of stars (see e.g. Fig. 9 of [608) ). which would 
then lead to low spins. These BHs would be in the optimal sensitivity window for 
LISA, but they are more difhcult to observe electromagnetically. In contrast, the 
Soltan [609j argument that SMBHs grow mainly be accretion (see [Appendix D ) really 
applies only to M > lO^M©, because only a small fraction of SMBH mass is in BHs 
with M < IO^Mq . As a result, higher-mass SMBHs could grow mainly by accretion, 
but lower-mass BHs would grow mainly by stellar disruption, mergers with stellar- 
mass compact objects and comparable-mass mergers. Therefore, the observation of 
spin as a function of mass could be a powerful diagnostic of SMBH evolution. 



Spins: observational estimates Spin estimates based on electromagnetic observations 
made enormous progress in the last three years. A summary of estimates available 
in the literature is provided in Table [3l Notice that in some cases (most notably 
for 4U 1543-475 and GRO J1655-40) different methods yield sensibly different spin 
estimates (see e.g. |612) ). The main methods used so far to estimate spins are 
continuum spectroscopy of accretion disks, spectroscopy of relativistically broadened 
Fe Ka fluorescence lines, and energetic arguments based on the radiative efficiency 
of quasars. A discussion of these topics would take us too far, but it is essential to 
appreciate the statistical and systematic errors involved in the numbers listed in Table 
[3l To improve readability, we review recent literature in this field in [Appendix D| 

A glance at the Table shows that the estimated spin magnitudes of stellar mass 
BHs cover the whole range from zero to one. This seems to confirm the doctrine that 
stellar- mass BHs in X-ray binaries essentially retain the spin they had at birth 589]. 
The situation is even more unclear for SMBHs, where (as shown in the last two lines 
of Table [3|) uncertainties in observational data on the mean efficiency of quasars lead 
to very different conclusions on the average values of SMBH spins. Gravitational wave 
detection could be instrumental in resolving this controversy. 



9.4- Detection range for Earth-based and space-based detectors 

We usually say that a gravitational wave signal is detectable when the SNR, as defined 
in Eq. (I152|) . is larger than some threshold, typically p > 10. Since the gravitational 
wave amplitude decreases linearly with the (luminosity) distance from the source, 
the distance corresponding to p = 10 is sometimes called the detector range. In the 
left panel of Fig. [T7| we show LISA's SNR for equal-mass inspirals (solid lines) and 
the ensuing ringdown (dashed lines) as a function of the total mass of the system 
M", where the superscript "0" means that the mass is measured in the source frame 
(see Ref. [TU] for details of the assumptions going into this calculation). Each line 
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Table 3. A list of spin estimates available in the literature, along with the method 
used for the estimate (see |Appcndix Dt and the relevant references. Tables 1 and 2 
of Ref. |610 | list spin estimates for 19 powerful FRII radio sources (FRIIb) and for 
29 central dominant galaxies (CDGs). For MS0735. 6+7421 Ref. [5T0) estimates a 
spin of 0.83 it 0.39, consistent with Ref. 611 . 



system 


Estimated spin 


Method 


Reierence 




Stellar- 


-mass BHs 




Cygnus A-1 


0.05 ± 0.01 


Line spectroscopy 


1612| 


T T\ /T/^ V O 


« 0.2 — 0.4 


Continuum 


[6131 


All /ITC 

4U 154o-4(^o 


O.O ±0.1 


Line spectroscopy 


[612) 




0.75 — 0.85 


Continuum 


[6141 


SAX J 171 1.6-3808 


0.6^0.4 


Line spectroscopy 


1612| 


XTE J1550-564 


« 0.1 - 0.8 


Continuum 


I6l3| 




0.76 ± 0.01 


Line spectroscopy 


|612| 


SWIFT J1753. 5-0127 


0-76Io.i5 


Line spectroscopy 


1615| 


~\ /TOO "V T 

Moo A- 7 


0.77 ± 0.05 


Continuum 


[6161 


XTE J1908H-094 


0.75 ±0.09 


Line spectroscopy 


1612| 


XI E J 1650-500 


0.79 ±0.01 


Line spectroscopy 


1612| 


GRS 1915-1-105 


0.7-0.8 


Continuum 


IMII 




0.98 - 1 


Continuum 


[61811619] 






Continuum 


[620j 


KjX ooy-4 


0.94 ±0.02 


Line spectroscopy 


[6211 1612] 


GRO J 1655-40 


> 0.25 


QPOs 


[622j 




0.65-0.75 


Continuum 


[618U614] 




w 0.1 -0.8 


Continuum 






0.98 ±0.01 


Line spectroscopy 


|612| 


XTE J1655-40 


K. 1 


Line spectroscopy 


|623| 


XTE J1550-564 


« 1 


Line spectroscopy 


(623] 




SMBHs 




29 CDGs 


0.1-0.8 


Energetics 


|610| 


19 FRIIb 


0.7-1 


Energetics 


|610| 


SWIFT J2127.4-t5654 


0.6 ±0.2 


Line spectroscopy 


|624| 


MCG-06-30-15 


U.t)Ot)_Q QQ2 


Line spectroscopy 


|607j 


1H0419-577 


K. 1 


Line spectroscopy 




MS0735.6-f7421 


« 1 


Energetics 


I6llj 




Large 


Average efficiency 


(6261 16271 16281, 16291 1630] 




Small 


Average efficiency 


(631il632, 633, 606, 



corresponds to a different luminosity distance or, which is the same, to a different 
cosmological redshift z: for example, a redshift z ~ 0.5 corresponds to a luminosity 
distance of about 3 Gpc in a standard ACDM cosmology. The plot illustrates a 
number of important points: (1) LISA can detect the last year of the inspiral of equal- 
mass binaries with total mass IQ^Mq < M° < IO^Mq out to cosmological distances 
(z > 10); (2) the ringdown phase has (under reasonable assumptions, which are 
confirmed by numerical simulations: see e.g. Fig. 14 in Ref. (634j ) a larger maximum 
SNR than the inspiral phase, and this maximum is achieved at larger values of the 
SMBH's mass. This is important, because it implies that ringdown searches are better 
suited for the detection of binaries with 10^ Mq < A/" < 10'^ Mq, which is closer to 
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Figure 17. Left: typical SNR for equal-mass inspirals and ringdown signals 
detected by LISA, assuming a final spin j = 0.8 and a conservative ringdown 
efficiency e^^ = 1% (from 1101 ). Top to bottom, the lines correspond to sources 
at redshift 2 = 0.54, 1, 2, 5, 10. Right: reach of Earth-based interferometers for 
a ringdown producing a BH with final spin j = 0.6 and e^d = 3%, as a function 
of the mass of the final hole (from |llj ) . 



the typical mass range estimated for SMBHs at galactic centres (see Section 1^75)) . 

In the right panel of Fig. [17] we show the ringdown detection range corresponding 
to p = 10 for Earth-based detectors. The plot shows that Advanced LIGO and the 
Einstein Gravitational Wave Telescope (ET) have the potential to detect ringdown 
from IMBH-IMBH systems of mass up to ^ 10"^ out to a luminosity distance of 
a few Gpc ^545j . In fact, third-generation Earth-based interferometers could probe 
the first generation of "light" seed BHs of M ^ 10^ — 10^ Mq, providing information 
complementary to LISA on the earliest BHs in the universe [635] . A more speculative 
source is the intermediate-mass ratio inspiral (IMRI) of stellar-mass BHs into IMBHs. 
For these systems, ringdown would be suppressed relative to inspiral because of the 
small mass ratio. However, what is lost in terms of number of cycles is gained in terms 
of detector sensitivity: IMRI ringdowns would happen in the optimal frequency band 
of second- and third-generation detectors, and therefore they could be detectable by 
a network of Advanced LIGOs [544j . 

9.5. Event rates 

Gravitational wave interferometers (unlike traditional electromagnetic observatories) 
respond to the waves' amplitude, not to their energy. Since the wave amplitude decays 
linearly with distance, a modest increase of (say) a factor two in sensitivity means that 
the detectable volume increases by a factor eight. Given that merging compact binaries 
are the most promising ringdown source, the relevant question for ringdown detection 
is then: how many merging compact binary events can we expect in a given volume? 
The issue of estimating event rates is one of the most pressing in gravitational wave 
detection. The uncertainties involved are so large, and progress in the field is so rapid, 
that any estimates we quote are likely to become rapidly obsolete. For this reason we 
dedicate little space to event rate estimates, providing a few references as a guide for 
the interested reader. 
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Stellar-mass and intermediate-mass black hole ringdowns In the simplest models, 
compact binary coalescence rates should be proportional to the stellar birth rate in 
nearby spiral galaxies, which can be estimated from their blue luminosity. Therefore 
the coalescence rates are usually given in units of Myr~^ L^q , where Lio is 10^" times 
the blue solar luminosity. To convert these numbers into detection rates, one must 
take into account the fact that detection ranges for ringdowns are different from those 
for inspirals. If the distance to an event is above ~ 50 Mpc the local over-density 
of galaxies can be ignored, and the number of galaxies containing possible sources 
is TV = [(47r)/3](D/Mpc)3(2.26)-30.0117 Milky-way equivalent galaxies, where 2.26 
is a correction factor to include averaging over all sky locations/orientations, and 
0.0117Mpc^'^ is the density of Milky- Way equivalent galaxies. For shorter distances 
one should use a sky catalog, such as Ref. [636j . 

For NS-NS binaries, early and conservative estimates were made by Phinney 
|637j . At present, the most reliable NS-NS merging rate estimates are obtained 
by extrapolating from observed binary pulsars [6381 I639j . Expected rates are « 
50MyT~^ L^Q , but they could be an order of magnitude lower or larger. For NS-BH 
and BH-BH rates we have to rely mostly on population-synthesis models |6401 1641) . 
Plausible rate estimates are « 2Myr~"'^Lj~Q^ for NS-BH binaries and « 0.4Myr~^Lj~Q'^ 
for BH-BH binaries, but they could be roughly two orders of magnitude larger or lower. 
These rates translate into tens to thousands of inspiral events per year in Advanced 
LIGO. The typical end-product of these mergers are BHs of mass ~ lOM©, and the 
range for Advanced LIGO detection of these BH ringdowns is more than an order 
of magnitude less than the inspiral range, so ringdown rates should be ~ 10^ — 10'^ 
times smaller than inspiral rates. Rates for ringdowns involving IMBHs are even more 
uncertain. In optimistic scenarios. Advanced LIGO could see ~ 10 IMBH binary 
mergers per year }545) and perhaps ~ 20 ringdowns from the merger of stellar-mass 
BHs into IMBHs [Mj (see also p2ll643] ). 

These predictions rapidly change as our understanding of the underlying physics 
and compact binary observations improve. For example, Ref. }644j argues that 
potential BH-BH binary progenitors may undergo a common envelope phase while 
the donor is evolving through the Hertzsprung gap. This would probably lead to a 
merger, thus shutting off a channel for BH-BH production and sensibly reducing BH- 
BH merger rates. On the other hand, based on observations of a very massive BH 
binary, Ref. |543j estimates an initial LIGO rate of order 0.5/yr for relatively massive 
BH binaries (that therefore would be observed mostly in the merger/ringdown phase). 

Particularly interesting for ringdown detection are compact binary mergers from 
dense star clusters. Ref. |645) points out that a "high-mass" selection occurs because, 
in nuclear star clusters at the centers of low-mass galaxies, three-body interactions 
usually "pair up" the two heaviest members of a triple system. This typically produces 
binary mergers with M > 20Mq , which means that Earth-based interferometers would 
usually observe the merger/ringdown phase. These findings are consistent with work 
by other authors [6461 [6471 [648] . 

Supermassive black hole ringdowns The LISA noise curve determines the optimal 
mass and redshift range where binary inspiral and ringdown events have large SNR, 
allowing a precise measurement of the source parameters. Reliable estimates of the 
number of events detectable during the mission's lifetime, and of their mass spectrum 
as a function of redshift, will play a key role in the planning of LISA data analysis. 
For this reason, over the last few years the calculation of SMBH merger event rates 
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and of their mass spectrum has become an active field of research. 



Table 4. SMBH binary rates (events/year) predicted by different models 
(adapted and updated from Ref. |567| ). 



Reference 


Rate 


Rcdshift range 


Haehnelt 2003 |649j 


0.1-1 
10-100 


< z < 5 (gas collapse only) 
z > 5 (hierarchical buildup) 


Enoki et al. 2004 [650 


1 


z > 2 


Mcnou et al. 2001 [651j 


10 


z < 5 


Rhook and Wvithe 2005 [6521 1653j 


15 


z-3-4 


Volonteri, Haardt and Madau '568, 570, 654 


- 30 


z - 4 - 16 


Begelman. Volonteri and Rees 600, |570| i654j 


- 20 


z-3-10 


Koushiappas and Zentncr 2005 [571] 


> 10^ 


mostly z ^ 10, down to z ^ 1 


Islam et al. 2003 [655J 


10^ - 10'^ 


z ~ 4- 6 



A discussion of rate estimates is out of the scope of this review (see e.g. |567[ 
I654[ 1656] ). but the large uncertainties in SMBH binary formation models and in the 
predicted event rates are quite evident from Table ID The numbers we list should be 
interpreted with caution. Each prediction depends on a large number of poorly known 
physical processes, and the notion of "detectability" of a merger event is defined in 
different ways: some authors define detectability setting a threshold on the SNR, 
others set a threshold on the gravitational wave effective amplitude. Furthermore, 
different authors use different LISA noise curves. A tentative bottom line is that we 
could face one of the following two scenarios. According to a class of models, we should 
observe ~ 10 events/year at redshifts (say) 2~z~6. However, we cannot exclude the 
possibility that hundreds or thousands of SMBH binaries will produce a large (and 
perhaps even stochastic) background in the LISA data. Clearly, the detection strategy 
to use strongly depends on which of the two scenarios actually occurs in nature. 

Besides being able to observe SMBH mergers throughout the universe, LISA 
should also be able to detect IMBH-IMBH binary mergers (that is, binaries containing 
a 10 - 100 M© BH orbiting a 100 - 1000 M© BH). Rates for IMBH binary detections 
were first estimated by Miller [65 7j and then revised by Will [658| . The revised 
estimates are very pessimistic, predicting ^ 10~^ events/year for typical values of 
the parameters. A more promising scenario involves an IMBH spiralling into a SMBH 
|659| . For these systems, Ref. |660| estimates a detection rate of a few events/year, 
suggesting that mergers of a lO'^ Mq IMBH into a 10^ Mq SMBH could be observed 
out to z '--^ 20 with an SNR of 10 in a one-year integration. Ref. ^661, predicts an even 
more optimistic rate of ~ 10^ events/year throughout the universe. These estimates 
are very preliminary and even more uncertain than the corresponding estimates for 
SMBH binaries, but they should be taken into account in design choices concerning 
(for example) the optimal armlength of LISA. 

9. 6. Inferring black hole mass and spin from ringdown measurements 

We have seen in Section 19.41 that the prospects for detection of ringdown radiation 
by LISA and advanced Earth-based detectors are quite encouraging. Interesting 
physics can be extracted from the observation of BH ringdowns [lOl [79l 146 Ij . Since 
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astrophysical BHs in general relativity are fully characterized by their mass and 
angular momentum, the detection of a single mode is in principle sufficient to estimate 
the mass and spin parameter of the hole, by inverting the experimentally determined 
Wim„(M, a/M),Qimn{a/M). Indeed, one finds that accurate measurements of SMBH 
mass and angular momentum can be made. 




Figure 18. Left: errors (multiplied by the signal-to-noise ratio p) in 
measurements of different parameters for the fundamental I = m = 2 mode as 
function of the rotation parameter a/M. Solid (black) lines give p(Ja/M i dashed 
(red) lines paj^i/M, dot-dashed (green) lines per a/ A, dot-dot-dashed (blue) lines 
pcr^, where crj. denotes the estimated rms error for variable k, M denotes the mass 
of the Btl, and A and denote the amplitude and phase of the wave (adapted 
from 1101 .) Right: actual errors in a LISA measurement for a source located ad 
Dl = 3 Gpc, with ringdown efficiency t^^ = 3%. 

For example, the left panel of Fig. [TSl shows the estimated error (multiplied by the 
SNR p) in measuring the SMBH mass Af , angular momentum parameter a/M, QNM 
amplitude A, and phase (j) for circularly polarized radiation from the fundamental 
I = m = 2 bar mode (cf. Eq. pSOp for definitions). If an energy ^ lO^^M is radiated 
into the fundamental mode of a IO^Mq SMBH with a/M = 0.8 at 3 Gpc {p ~ 200), M 
and a/M could be measured to 1%; if the energy deposition is only 10^^, they could 
still be measured to 10%. These numbers were computed for LISA in Ref. [TO], but 
they carry over to other detectors through a simple rescaling by p. Generalizing to 
multi-mode detection (and specifically to the detection of two modes with a range of 
relative amplitudes) one finds similar results ^iDJ. Gravitational wave detectors will be 
able to determine the mass and spin of BHs with excellent precision from observations 
of the ringdown phase. 

Event loss and bias in parameter estimation using single-mode templates Current 
ringdown searches are performed using matched filtering and single-mode templates, 
consisting of a single exponentially damped sinusoid. These are the simplest possible 
templates, and they are expected to capture the relevant physics of the problem 
when one ringdown mode dominates over the others. Unfortunately, this notion of 
dominance must be precisely formulated. Consider for illustration the case of non- 
spinning binary BH mergers, and suppose for simplicity that there are only two modes 
in the signal, say I = m = 2 and Z = m = 3. In particular, assume that the strain h 
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as seen by a detector, Eq. (I151|) . is of the form 

h = Aie-'' sin(27r/i t-(/)i) + y^2e"''^'*/'^' sm(27r/2i-(/)2) .(160) 

Estimates of the multipolar energy distribution give a relative amplitude A2/A1 — 
^33//i22 ~ 0.3 — 0.4(1 — 1/q), where q is the mass ratio and /122 (h^^) are the amplitudes 
of the / = m = 2 (3) modes of the radiation, respectively 44]. For g > 3, which includes 
most likely merger scenarios, we get /133//122 > 0.2 — 0.3. It is now natural to ask: 
given a relative amplitude of this order, how many events would we miss in a search 
with single-mode ringdown templates? 




<t', *i 



Figure 19. Left: contour plots of the fitting factor as a function of tiie phase 
angles of the signal. For this illustrative calculation we assume that the Kerr 
parameter of the final BH is j = 0.6, and that the relative amplitude of the 
second mode is ^ = 0.3. Right: quality factor estimated by a single-mode filter. 
The "true" frequencies and quality factors for a/M = 0.6 are Mu>[ii = 0.4940, 
Qi = 2.9490, MujR2 = 0.7862, Q2 = 4.5507. See [U] for more details. 

The answer is quantified in the left panel of Fig. [HI where we show Apostolatos' 
fitting factor (FF) [6621 1663| resulting from searching a two-mode signal with single- 
mode templates from a 200Mq, a/M = 0.6 BH, with relative amplitude A2/A1 = 0.3, 
in Advanced LIGO data. Contour plots of the FF are shown as a function of the 
(unknown) phases in Eq. (jl60p . The FF is essentially the ratio p/pmaxj where p is the 
actual SNR achieved by matched filtering and Pmax is the maximum possible SNR, 
attained when the template and waveform coincide. FFs larger than 0.97, leading to 
a loss of less than about 10% in events, are typically considered acceptable. The FF 
achievable by a detection is usually assumed to correspond to the minimum of the FF 
in the space of the unknown parameters (0i, 02 ) (this is known as the "minimax" 
criterion). Therefore, according to this simplified calculation, single- mode templates 
may well lead to unacceptable SNR degradation in a search. The situation can be 
even worse for larger masses and other detectors [TT]; similar conclusions apply also 
to LISA, albeit in a completely different mass range. 

Besides leading to a loss in the number of events, single-mode templates lead to 
a large bias in the estimation of the BH's mass and spin. The estimated frequency 
has relatively small bias, and it always corresponds to the dominant (least-damped) 
mode in the pair. Results are more interesting for the estimated quality factor, shown 
in the right panel of Fig. [THIas a function of the phase angles. For our chosen value 
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of the Kerr parameter the quality factor of the I — m — 2 and / = m = 3 modes 
are Qi — 2.9490 and Q2 — 4.5507, respectively. Comparing with the left panel, we 
see that relative minima in the FF (white "islands" in the left panel) occur, roughly 
speaking, when the quality factor "best fits" the subdominant, I = m = S mode. This 
is a rather remarkable result: the minimax filter "best fits" the subdominant mode 
in the pair, leading to significant bias in the estimation of the quality factor (and 
hence of the spin). Unfortunately, maxima in the FF do not correspond to the filter 
being optimally adapted to the I — m = 2 mode. As the filter tries to maximize the 
SNR, the estimated value of the quality factor becomes significantly biased, and it 
deviates quite sensibly from the value expected for the dominant {I = m = 2) mode. 
The bottom line is, again, that single-mode filters may be useful for detection, but a 
multi-mode post-processing will be necessary for accurate spin measurements. 

There is, of course, a price to pay when using multi-mode templates. About 
A/" ~ 6(5inaxlog/max//min 500 single-mode templates are enough to cover the 
parameter space for Earth-based detector searches, if we assume an event loss of 
no more than 10% (i.e., a minimal match larger than 0.97 [663] ). For two-mode 
templates, rough estimates suggest that this number may increase up to A/" ^ 5 x 10^, 
with 6 ~ 1 a detector-dependent constant [TT]. A more detailed data analysis study 
(e.g. using better template placing techniques, along the fines of [472] I473| I474|, [T3]) 
will be needed to reduce computational requirements. 



. 7. Tests of the no-hair theorem 



The fact that all information is radiated away in the process leading to BH formation, 
so that astrophysical BHs in Einstein's theory are characterized completely by their 
mass and angular momentum, is known as the "no-hair theorem". To test this 
theorem, it is necessary (but not sufficient) to resolve two QNMs [TO l 146 Ij . Roughly 
speaking, one mode is used to measure M and a, and the other to test consistency 
with the Kerr solution. 

Can we tell if there really are two or more modes in the signal, and can we 
resolve their parameters? Physical intuition suggests that if the noise is large and the 
amplitude of the weaker signal is very low (or if the two signals have almost identical 
frequencies) the two modes could be difficult to resolve. A crude lower limit on the 
SNR required to resolve frequencies and damping times was presented in Refs. [lOlITT]. 
The analysis uses the statistical uncertainty in the determination of each frequency 
and damping time, which a standard Fisher-matrix calculation estimates to be 10, 
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These errors refer to mode "1" in a pair; errors on /2 and T2 are simply obtained by 
exchanging indices (1 ^ 2) [10l|Tl]. The expressions above further assume white noise 
for the detector, and that modes "1" and "2" correspond to different values of I or m. 

A natural criterion {d la Rayleigh) to resolve frequencies and damping times 
is that l/i — /2I > max(tT/j , cTy^ ), |ti — T2I > max(CT7.j, CT-a). In interferometry this 
would mean that two objects are (barely) resolvable if the maximum of the diffraction 
pattern of object 1 is located at the minimum of the diffraction pattern of object 2. 
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We can introduce two "critical" SNRs required to resolve frequencies and damping 
times, p(^■^ = max(/9cr/i , pa/J/l/i - /2I ,Pcrit = niax(p(T^i , pcr^J/|Ti - T2I, and recast 
our resolvability conditions as 

P > Pent min(pf^it, , (163) 

P > Pboth = max(p4.it> Pcrit) ■ (164) 

The first condition implies resolvability of either the frequency or the damping time, 
and the second implies resolvability of both. 

Now let us consider how to resolve amplitudes, i.e. the minimum SNR needed to 
determine whether two or more modes are present in a given ringdown signal. Suppose 
again, for simplicity, that the true signal is a two-mode superposition. Then we expect 
the weaker signal to be hard to resolve if its amplitude is low and/or if the detector's 
noise is large. Appendix B of Ref. [TT] quantifies this statement by deriving a critical 
SNR for amplitude resolvability pglrt based on the generalized likelihood ratio test. 
The derivation of this critical SNR is based on the following simplifying assumptions: 
(i) using other criteria one has already decided for the presence of at least one damped 
exponential in the signal, and (ii) the parameters of the ringdown signal (frequencies 
and damping times), as well as the amplitude of the dominant mode, are known. In 
practice the latter assumption is not valid, so these estimates of the minimum SNR 
should be considered optimistic. Fig. [20l compares the critical SNR pglrt (Eq- (B12) 
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Figure 20. Minimum SNR required to resolve two QNMs from a BH resulting 
from the inspiralling of a BH binary with mass ratio q. The dominant mode (mode 
"1") is assumed to be the fundamental QNM with I = m = 2; mode "2" is either 
the fundamental I = m = 3 QNM (solid lines) or the fundamental I = m = A QNM 
(dashed lines). The relative mode amplitude A{q) is estimated from numerical 
simulations 1441 . If p > PGLRT we can tell the presence of a second mode in the 
waveform, if p > Pcrit we can resolve either the frequency or the damping time, 
and if p > Pboth we can resolve both. 

of Ref. [TT]) and the two different criteria for frequency resolvability, Eqs. p63p and 
(|164p . The plot shows that pcrit < Pglrt < Pboth for all values of q. Therefore, given a 
detection, the most important criterion to determine whether we can carry out no-hair 
tests is the GLRT criterion. If p > pqlrt we can decide for the presence of a second 
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mode in the signal. Whenever the second mode is present, we also have p > Pcrit^ 
that is, we can resolve at least the frequencies (if not also the damping times) of the 
two modes. A SNR p ~ 30 — 40 is typically enough to perform the GLRT test on the 
Z = TO = 3 mode, as long as q > 1.5. From Figs. [1711^ we conclude that not only 
LISA, but also advanced Earth-based detectors (Advanced LIGO and FT) have the 
potential to identify Kerr BHs as the vacuum solutions of Finstein's general relativity. 

In conclusion, ringdown radiation can be used to distinguish BHs from exotic 
alternatives, such as boson stars [463j or gravastars [664] • Ringdown tests of the Kerr 
nature of astrophysical BHs are independent from (and complementary to) proposed 
tests using a multipolar mapping of the Kerr spacetime, as encoded in FMRI signals 
according to "Ryan's theorem" and its generalizations [6651 16661 1667] . 

9.8. Matching inspiral and ringdown: problems and applications 

Models of gravitational waveforms from inspiralling compact binaries usually rely on 
the post-Newtonian approximation to general relativity [668^. Ideally, for matched- 
filtering detection we would like to have detection template banks with "complete" 
waveforms encompassing the inspiral, merger and ringdown. Phenomcnological 
template families based on physically motivated fits of numerical merger waveforms 
or on effective-one-body models are now available [6691 16341 16701 1671j . However, it 
would be desirable to have a full understanding of the merger process, connecting the 
post-Newtonian approximation during the inspiral to a description of the ringdown as 
a superposition of QNMs. 

QNM fits of numerical relativity waveforms are routinely used to check that the 
total angular momentum of the system is conserved during BH merger simulations. 
This is usually achieved by computing three independent quantities: the angular 
momentum radiated at infinity, the angular momentum obtained from isolated horizon 
calculations, and the angular momentum obtained by inverting the frequencies and 
damping times resulting from a QNM fit of the ringdown of the final Kerr BH (see 
e.g. t477i |44» £72] ) . Buonanno, Cook and Pretorius noticed that, as more and more 
overtones are included, a QNM expansion gets in better agreement with numerical 
waveforms for equal-mass binaries j477j . 

So far, efforts to "stitch" the post-Newtonian approximation (or one of its 
effective-one-body variants) to the ringdown have used rather crude models for 
ringdown excitation. The original effective-one-body approach simply attached the 
plunge waveform to the fundamental I = m = 2 QNM of a Kerr BH [673J. In 
their study of recoils, Damour and Gopakumar attached the plunge waveform to the 
fundamental Schwarzschild QNMs with I = 2 and I = 3 [674]. A comprehensive 
analysis of Pretorius' numerical waveforms for equal-mass BH binaries |477| clearly 
illustrated the importance of higher-order overtones. Damour and Nagar obtained a 
good match to numerical relativity waveforms by requiring continuity of the plunge 
waveform with a (Schwarzschild) ringdown waveform including five overtones on a grid 
of points, or "comb" 176J. A "hybrid comb" procedure, imposing the continuity of 
the waveform and its derivatives, was later introduced in Ref. [174j . These stitching 
methods have great phenomcnological interest, and they can do remarkably well at 
reproducing numerical waveforms. For example, the authors of Ref. |174j extend an 
effective-one-body waveform through merger by stitching the inspiral waveform to 
QNMs, finding striking agreement (at the 0.1% level) between the numerical QNM 
frequency and the perturbative prediction for the same frequency computed from the 
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final mass and spin of the numerical simulation. A slightly different approach was 
adopted in Ref. |670j . where the frequency-domain waveform was truncated abruptly 
at some cutoff frequency close to the QNM frequency. The truncated waveforms 
display a Gibbs phenomenon in the time domain, and they are (surprisingly) quite 
effective at detecting the full signal. The reason is that a cutoff in the frequency 
domain around the fundamental QNM frequency corresponds, in the time domain, to 
a damped sinusoid with frequency close to the cutoff frequency. 

These phenomenological "stitching recipes" may well be enough for the purpose 
of gravitational wave detection, but it is desirable to have a complete understanding 
of QNM excitation in a compact object merger within perturbation theory, based 
on the concept of excitation coefficients. Ref. T6T carried out a first step in this 
direction by computing the excitation factors for Kerr BHs (see Section 13. 2p , but 
more work is required to compute the excitation coefficients for generic initial data, 
and to understand the validity of the linear approximation in a BH merger. A 
correct matching of effective-one-body waveforms with the ringdown signal may have 
astrophysical implications because QNMs play an important role in the recoil of 
merging BH binaries, through a process sometimes called "ringdown braking" (see 
|675| 16741 1676] and references therein). A complete analytical description of the 
merger of spinning, precessing binaries would also prove useful for statistical studies 
of astrophysical relevance [677j . 

Even if phenomenological waveforms may be good enough at detecting 
gravitational wave signals, a complete description of the waveform (including both 
inspiral and ringdown) has been shown to improve parameter estimation by effectively 
decorrelating the source parameters (see |534j for a study predating the numerical 
relativity breakthrough, and [6781 16791 16691 16341 I535j for recent efforts in this 
direction). Large-scale efforts (dubbed the NINJA and Samurai projects, respectively) 
are now attempting to use complete numerical waveforms in LIGO data analysis, and 
to cross- validate numerical waveforms produced by different codes |525| 1680) . 

An important point to keep in mind is that inspiral and merger/ringdown 
really probe different BH populations, so they provide complementary information. 
This should be quite clear by inspecting the red and black histograms in Fig. (THl 
which show that the spin distribution of SMBHs in binaries is usually different from 
the spin distribution of the remnant BHs formed as a result of the merger. This 
complementarity can be used for interesting physical applications. For example, if we 
can determine accurately enough the masses and spins of BHs before and after merger 
for the same system, we could use this information to test Hawking's area theorem in 
astrophysical settings [681] . 

10. Other recent developments 

10.1. Black hole area quantization: in search of a log 

One of the main driving forces behind the development of new tools to study QNMs 
is the possibility that classical BH oscillations could yield insights into their quantum 
behavior. First suggested by York [68j and Hod [17], this idea was further explored 
by Dreyer [97] in the context of Loop Quantum Gravity, and subsequently revisited 
by many authors (see e.g. [551 1271 [Ml MM MB MB [MZl MB [Ml MB [MI])- 

The idea can be traced back to arguments by Bekenstein and collaborators 
|6921 1693j . A semi-classical reasoning suggests that the BH area spectrum is quantized 
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according to 

An^jII, iV=l,2,... , (165) 

where Ip is the Planck length (but see |694j for criticism). Statistical physics 
arguments impose a constraint on 7 |6931 1695j : 

7 = 4 log fc , (166) 

where k is an integer, which in general is left undetermined (but see [696] ) . 
Inspired by Bekenstein's ideas, Hod [17] proposed to determine k via a version of 
Bohr's correspondence principle in which the highly-damped QNM frequencies play 
a fundamental role. At the time, the only available exploration of highly-damped 
QNMs was the numerical study by NoUert [90J, which indicated that gravitational 
highly-damped QNMs in the Schwarzschild geometry asymptote to 

Afcj = 0.0437123- i(2n+ l)/8. (167) 

While looking for classical oscillation frequencies proportional to the logarithm of an 
integer. Hod realized that 0.0437123 ~ ln3/(87r), and went on to suggest that the 
emission of a quantum with frequency ujr = ln3/(87r) corresponds to the smallest 
energy a BH can emit. The corresponding change in surface area would then be 

AA = 327rMdM = i2TiMhw = Ahhii . (168) 

and by comparing with Eq. (|165p we get fc = 3, thus fixing the area spectrum to 

An =4/|ln3n; n = 1,2,... (169) 

A few years later Dreyer [97J used similar arguments to fix the Barbero-Immirzi 
parameter of Loop Quantum Gravity. Shortly afterwards, Motl [Qlj and Motl and 
Neitzke |101| showed analytically that the highly damped QNM frequencies of a 
Schwarzschild BH are indeed given by Eq. (j86p . containing the desired logarithm 
of an integer. 

Motl and Neitzke's work lent some support to Hod's ideas. It was followed by a 
flurry of activity to explore the highly-damped QNMs of several BH spacetimes, whose 
outcome has been described and summarized in previous sections. Unfortunately 
the conjectured relation between QNMs and area quantization stumbled before the 
charged and rotating four-dimensional geometries, for which the highly-damped regime 
is not as simple as suggested by Hod's original argument (see Sections 15. 3115. 2p . 
Recently Maggiore |688] observed that Hod's prescription was based on experience 
with "normal" quantum systems, for which the relevant frequency is ljr. For highly- 
damped systems, Maggiore noted that one should rather consider the imaginary part 
of the QNMs lji, solving a number of puzzles and obtaining a different expression 
for the quantum of area. A similar prescription has been extended also to Kerr BHs 
|689| 1690] and other geometries |691] . It is not clear whether Maggiore's suggestion 
can be extended consistently to all geometries. Whether a relation between QNMs 
and the quantum behavior of BHs exists or not. Hod's suggestion was at the very least 
an important thrust to complete our understanding of classical BH oscillation spectra. 

10.2. Thermodynamics and phase transitions in black hole systems 

In the last few years, remarkable relations between classical and thermodynamical 
properties of black objects have been uncovered. For instance, a correspondence 
between classical and thermodynamical instabilities of a large class of black branes 
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conjectured by Gubser and Mitra [697| I698j was proved by Reall |699j (see |700j for a 
review) . Manifestations of this duality are expected to appear in the QNM spectra. 
Indeed, some indications that phase transitions correspond to changes in the QNM 
spectrum were provided in specific cases by various authors [7011 17021 17031 17041 17051 
I706j . However, at present there seems to be no obvious correspondence between 
thermodynamical phase transitions of the kind suggested by Davies |707| I708 j and 
QNM spectra (see |709j and references therein). 

Another connection between QNMs and BH thermodynamics may follow from 
Hod's proposed "universal relaxation bound" [7101 1711[ I712j . Hod's proposal asserts 
that the relaxation time t of any thermodynamic system is bounded by r > h/ijrT), 
where T is the temperature of the system. For BHs this implies the existence of 
(at least one) QNM frequency with imaginary part uji < ttT, where T is now the 
Hawking temperature of the BH. The bound seems to be valid for various kinds of 
BHs, and it may be saturated by extremal BHs. The significance of Hod's bound is 
not completely clear to us. The bound is trivially satisfied by any physical system 
exhibiting hydrodynamic behavior, since such a system always possesses sufficiently 
long-lived modes (and, correspondingly, QNMs with imaginary part sufhciently close 
or even infinitely close to zero: these hydrodynamic frequencies are discussed in Section 
18. 3p . Finite- volume systems might be more interesting, but then the concept of the 
"relaxation time" used in the bound needs a proper definition. In any case, this is an 
interesting idea which might require better understanding. The relation between Hod's 
bound and the viscosity-entropy bound (see Section 18. 4p was discussed in Ref. [713j . 

10.3. Non-linear quasinormal modes 

QNMs are usually defined and studied by considering only first-order perturbations. 
Being an intrinsically non-linear theory, general relativity is expected to display non- 
linear effects, which might conceivably be captured by going to higher orders in 
perturbation theory. The second-order formalism laid down several years ago by 
various authors [85l I714( I715j has recently been used to compute corrections to the 
QNM frequencies, their detectability and their influence on the late-time behavior of 
the system |7161 17171 1718) . The encouraging outcome of these studies is that non- 
linear effects may well be observable by future gravitational wave interferometers. 
Favata [719j explores the interesting possibility of detecting another important non- 
linear effect of general relativity (the so-called gravitational wave memory) through 
gravitational wave observations of merging binaries. 

Higher-order perturbations of BH spacetimes have been explored systematically 
by Brizuela et al. [7201 1721 ] , using the Gerlach-Sengupta formalism and the computer 
algebra methods described in [722J (see also Kol |723j ). A complete gauge-invariant 
formalism for second-order perturbations of Schwarzschild BHs was recently reported 
[724] . 

10.4. Quasinormal modes and analogue black holes 

Strong-field effects of general relativity are very small in Earth-bound experiments. 
For this reason BH physics is most easily studied via observations of astrophysical 
phenomena such as accretion. X-ray spectra, and hopefully gravitational wave 
emission. The possibility to devise gravitational experiments probing strong-field 
general relativity in the lab, as appealing as it sounds, may seem out of reach. However, 
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Unruh, Visser and others [7251 [3T] ingeniously showed that some defining properties 
of BHs can be reproduced and studied by "analogue BHs" . These systems display at 
least a subset of the properties traditionally associated with BHs and event horizons. 
Unruh's analogue BHs do not carry information about the dynamics of Einstein's 
equations, but share many kinematical features with true general-relativistic BHs. 

The basic idea behind these analogue BHs is quite simple. Let us focus on a 
particular analogue BH, the acoustic or "dumb" hole [725\ I726j . Consider a fluid 
moving with a space-dependent velocity Vq{x^), for example water flowing through a 
variable-section tube. Suppose the fluid velocity increases downstream, and that there 
is a point where the fluid velocity exceeds the local sound velocity c{x^) = y^dp/dp, 
in a certain frame. At this point, in that frame, we get the equivalent of an apparent 
horizon for sound waves. In fact, no (sonic) information generated downstream of 
this point can ever be communicated upstream (for the velocity of any perturbation 
is always directed downstream, by simple velocity addition). This is the acoustic 
analogue of a BH, sometimes referred to as a dumb hole. These are not true BHs, 
because the acoustic metric satisfies the equations of fluid dynamics and not Einstein's 
equations. However, sound waves propagate according to the usual curved space Klein- 
Gordon equation with the effective metric [7251 17261 [3T| , 
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Analogue BHs should Hawking-radiate, though an experimental verification of 
Hawking radiation in the lab is not an easy feat [727[ I31j . Furthermore, sound wave 
propagation in these metrics should reproduce many features of wave propagation 
in curved spacetimes. Most importantly for this review, acoustic BHs have a (this 
time literally!) "characteristic sound" [5] encoded in their QNM spectrum. QNMs 
of acoustic BHs, which may be important in experimental realizations of the idea, 
were computed in Refs. [7281 1240[ 17291 [75D] for a simple (2 -I- 1)— dimensional acoustic 
hole, the "draining bathtub" [726] . for which vq — (— + B(t>)/r. The modes of 
(3 -|- l)-dimensional acoustic holes with a "sink" at the origin were computed in 
Refs. [7281 1731[ I732j . Scattering from these holes is discussed in Ref. |733j . 

The experimental confirmation of these predictions is an interesting topic for 
future research. In practice, one may need a device to accelerate the fluid up to 
supersonic velocities, such as a Laval nozzle [734[ 1735] (see also Unruh's discussion in 
[727] ). Laval nozzles were first used in steam turbines, but they find applications in 
other contexts, including rocket engines and nozzles in supersonic wind tunnels. They 
consist of a converging pipe where the fluid is accelerated, followed by a throat (the 
narrowest part of the tube) where the flow undergoes a sonic transition, and finally 
a diverging pipe where the fluid continues to accelerate. QNMs of flows in de Laval 
nozzles were discussed and computed in Refs. [7361 1737] . 

Following on Unruh's "dumb hole" proposal, many different kinds of analogue 
BHs have been devised, based on condensed matter physics, slow light etcetera. We 
refer the reader to Refs. [727] [31] for thorough reviews on the subject. QNMs of 
condensed-matter analogue BHs were recently computed in Refs. [738] I739j . Finally, 
we should mention that accretion of material onto astrophysical BHs can give rise 
to supersonic walls, i.e., acoustic horizons outside the event horizon. This process 
is particularly interesting: astrophysical BHs may provide a nature-given setting for 
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producing analogue BHs |740|, I741j . 
11. Outlook 

The investigations of the last decade show that quasinormal spectra encode a wealth of 
information on the classical and quantum properties of BHs and black branes. From 
a holographic, high-energy perspective, BH QNMs yield important information on 
the quasiparticle spectra and transport coefficients of the dual theory, and have an 
intriguing hydrodynamic description. From an astrophysical viewpoint, gravitational 
wave measurements of QNMs may allow us to accurately measure BH masses and spins 
with unprecedented accuracy and to test the no-hair theorem of general relativity. The 
suggestion that BH oscillation frequencies may be related to their quantum properties 
is controversial, but at the very least it has stimulated tremendous technical progress in 
the calculation of previously unexplored regimes of the QNM spectrum. In this review 
we tried to summarize these recent developments. There is no doubt that further 
interesting connections between QNM research, fundamental physics and astrophysics 
will be unveiled in the future. In closing, we think it might be useful to list some of 
the outstanding, unresolved issues. 

General mathematical problems, asymptotic expansions and related issues - It would 
probably be fair to say that the level of mathematical rigor in QNM studies deserves 
an improvement. Whether or not boundary value problems for dissipative systems 
can be treated in the framework of a sufficiently general mathematical theory appears 
to be an open question. Less ambitiously, a unified treatment of resonances (QNMs 
included) at the level of mathematical physics would be desirable. Reliable regular and 
asymptotic expansions based on systematic procedures rather than ad hoc recipes have 
not yet been developed in many cases. Improving numerical algorithms, in particular 
for QNMs with large imaginary parts, also remains an avenue of research. 
Black hole perturbation theory in asymptotically-flat spacetimes - Despite a fifty-year 
long history of investigation, BH perturbation theory in asymptotically flat spacetimes 
has a number of unresolved issues. The problem of metric reconstruction from the 
Teukolsky formalism is still open. QNM excitation for generic perturbations of Kerr 
BHs have not been fully explored yet. The eikonal limit of Kerr QNMs is still 
poorly understood, as is the relation of QNMs with generic (/, m) to unstable circular 
geodesies in the Kerr metric. Finally, an outstanding problem concerns the decoupling 
of linear perturbations of the Kerr-Newman metric. 

Gauge- gravity duality - At the moment, computing QNMs in the context of the gauge- 
gravity duality is a thriving industry, naturally expanding to include models with 
higher-derivative gravity, holographic models with spontaneous symmetry breaking, 
or backgrounds dual to non-relativistic theories. As discussed in Section [51 one 
distinguishes hydrodynamic-type modes from the rest of the QNM spectrum. On 
a technical level, these modes are characterized by small parameters and thus can 
often be determined analytically. A full analysis of hydrodynamic QNMs for generic 
backgrounds is currently lacking. In particular, it would be helpful to obtain an 
expression for the parameters of the sound mode similar to Eq. (|148|) for the diffusion 
constant. The spectra of fermionic fluctuations (e.g. of the Rarita-Schwinger field, 
see j417j ) in black brane backgrounds have not yet been fully computed. It would be 
very interesting to explore further the connection between holography and the black 
hole membrane paradigm [117[ I742[ I114[ 1115] . 

Gravitational wave astronomy - An important problem in gravitational wave physics 
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concerns the matching of inspiral waveforms (as computed by post-Newtonian theory) 
to the ringdown phase. This is necessary to produce complete phenomenological 
templates to be used in gravitational wave detection. Most attempts at solving 
this problem have adopted a purely phenomenological approach, but a study of the 
excitation coefficients induced by generic initial data (possibly considering non-linear 
corrections) could provide a more consistent and systematic solution to the matching 
problem. A study of the QNM contribution to the Green function has recently been 
shown to hold promise as a computational approach to solve the self-force problem 
|172| . which is of fundamental importance to model extreme- and intermediate-mass 
ratio inspirals for LISA data analysis. Other open issues concern ringdown data 
analysis. Simple calculations show that single-mode templates only produce moderate 
losses in SNR when detecting multi-mode signals. However, the loss in terms of 
parameter estimation accuracy is much more significant. The problem of optimal 
template placing for detection by multi-mode templates needs to be addressed. The 
potential of higher multipoles in estimating the source parameters without angle- 
averaging should also be explored: the different angular dependence of subdominant 
multipoles may provide information on the spin direction of the final BH. If BH spins 
are linked to jets, this information could be used for coincident searches of gamma-ray 
bursts or other electromagnetic counterparts to compact binary mergers j522i 1523) . 

In the last few years, unexpected connections between QNMs and seemingly 
unrelated phenomena (such as analogue BHs or thermodynamical BH phase 
transitions) have been uncovered or proposed. More intriguing connections will surely 
emerge in the coming decades. 
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Appendix A. 

Isospectrality, algebraically special modes and naked singularities 

A number of striking relations among gravitational perturbations of the Kerr geometry 
were revealed by Chandrasekhar and colleagues 11411 [5^ I142[ I743j . Here we review 
these relations specializing most of the discussion to the Schwarzschild geometry. The 
reader is referred to the original works for more details and for extensions to the Kerr 
spacetime [22l [Mil [Ml EM ■ 
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By direct substitution, it can be checked that the gravitational potentials 



± 



=2 



in Eqs. (|14p and p5|) can be rewritten in the form 
where 



3M(6Mr2 + 2AL2(2M-r)) A(A + 1) 3M , 

W ^ ^ — . (A. 2) 

Ar2 (6Af + 2Ar) i2 3^ ^ ' 

(these equations correct some typos in Ref. |130j ). Eq. (IA.1|) emerged from 
Chandrasekhar's investigations [62l 11411 I142j of the nature of the gravitational 
potentials. Potentials of this form are called superpartner potentials [744], and they 
imply the following relation between the corresponding wavefunctions and ^'^ 

= -J—(^W^^ + ^) . (A.3) 



(3 — uj^ \ dr^ J 

Eq. (jA.ip seems to be unique to four-dimensional spacetimes, and does not generalize 
to higher dimensions [102[ I103j . The potential for electromagnetic- type and 
gravitational-type perturbations of extremal RN BHs can also be expressed in the 
form (jA.ip [2591 1260] . This justifies the fact that electromagnetic perturbations with 
angular index / are isospectral to gravitational with index l — l, as discussed in Section 
15.21 The isospectrality is a manifestation of supersymmetry between electromagnetic 
and gravitational perturbations for extremal charged BHs [2591 12601 1261] . 
Isospectrality and asymptotic flatness Suppose that u; is a QNM frequency for i.e. 
that ^+ A+e"*"''* when r — > r+ and '5+ — > a^^jC*'^''* when r 00. Here a^J^^ = 1 
for the Schwarzschild metric and, if we impose Dirichlet boundary conditions, a^J^j ~ 
for the SAdS metric. Substituting into Eq. (jA.3D we see that 

*----^(T4-(r+)-zc.)e-^'^'^*, r ^ r+ , (A.4) 
W{oo)e"^''' + —— , r^oo. (A.5) 



— a;^ \ °" dr 

The key point is that d'^~ /dr^ ~ g^^^r, infinity for the Schwarzschild geometry. 
Therefore, if w is a QNM frequency for it is also a QNM frequency for VP", 
and the two types of gravitational perturbations are isospectral [BU I141j . In general 
d^~ /dr^ ^ for SAdS, so the isospectrality is broken. The above relations are valid 
also for dS backgrounds: since the outer boundary conditions are imposed at the 
cosmological horizon, it is easy to see that gravitational perturbations of both parities 
are again isospectral. A specialized analysis is needed at the points where P — uj'^ = 0: 
this condition defines the so-called algebraically special modes, discussed below. 
Algebraically special modes A class of special modes can be found analytically using 
Eq. (|A.1[) . In the Petrov classification, this condition corresponds to a change in the 
algebraic character of the spacetime. For this reason, Chandrasekhar called these 
modes "algebraically special" (AS) [1421 Defining ±W — -£-logx±, the wave 
equation can be written as 

1 d^^± , / 2 1 rf'x± 



[uj- 



^± drl x± dr\ 
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where we have used the identity ^^P^ ~ niodes have frequencies 

oj = (li such that [3 — Clf — 0. In this case Eq. (|A.6|1 can be easily integrated: 

*± = X± y dr,/xl = X± {c±i + C±2 y dr^/xl^ , (A.7) 

where x± = exp [± / VFdr*] . The special relation between the two gravitational 
potentials extends to the Kerr geometry, where AS modes correspond to a vanishing 
Teukolsky-Starobinsky constant [22]. The nature of the boundary conditions at the 
Schwarzschild AS frequency is extremely subtle, and has been studied in detail by 
Maassen van den Brink [745] . Let us introduce some terminology [745| : 

1) "standard" QNMs have outgoing-wave boundary conditions at both sides (that is, 
they are outgoing at infinity and "outgoing into the horizon"); 

2) total transmission modes from the left (TTM^'s) are incoming from the left (the 
BH horizon) and outgoing to the right (spatial infinity); 

3) total transmission modes from the right (TTM^'s) are incoming from the right and 
outgoing to the left. 

The Regge- Wheeler equation and the Zerilli equation must be treated on different 
footing at the AS frequency, since the supersymmetry transformation leading to the 
proof of isospectrality is singular there. It turns out that the Regge- Wheeler equation 
has no modes at all, while the Zerilli equation has both a QNM and a TTM^ at the 
AS frequency [745] . 

Numerical calculations of AS modes have yielded some puzzling results. Studying 
the Regge- Wheeler equation, Leaver [75] found a QNM very close (but not exactly 
at) the AS frequencies of Eq. ([83ll . Namely, he found QNMs at frequencies 2MVl'2 — 
0.000000 -3.998000i and 2MVt'^ = -0.000259 - 20. 015653i. Similarly, in the extremal 
RN case one finds a QNM frequency very close to, but not exactly equal to, the 
AS frequency [241] . Maassen van den Brink ^745| speculated that the numerical 
calculations may be inaccurate and that no conclusion can be drawn on the coincidence 
of and "if the latter does exist at all". 

An independent calculation was carried out by Andersson [746], who found 
that the Regge- Wheeler equation has pure imaginary TTM^'s very close to f2; for 
2 < ^ < 6. He therefore suggested that the modes he found coincide with Vli, which 
would then be a TTM. Maassen van den Brink 745] observed again that, if all figures 
in the computed modes are significant, the coincidence of TTM's and QNMs is not 
confirmed by this calculation, since (l\ and (li are numerically (slightly) different. 

Onozawa [230] calculated the (TTM) AS mode for Kerr BHs, improving upon the 
accuracy of the Kerr AS frequencies computed in Ref. ^142^ . He showed that the Kerr 
QNM with overtone n = 9 tends to the AS frequency (as defined by the Teukolsky- 
Starobinsky identity) when a ^ and suggested that QNMs approaching Q.i from the 
left and from the right may cancel at a = 0, leaving only a special (TTM) mode. The 
situation concerning Kerr QNMs branching from the AS Schwarzchild mode is still far 
from clear. Using slow-rotation expansions of the perturbation equations, Maassen 
van den Brink [745J drew two basic conclusions on these modes. The first is that, for 
a > 0, the so-called Kerr special modes are all TTM's (left or right, depending on 
the sign of s). The TTMj^'s should not survive as a ^ 0, since they do not exist in 
the Schwarzschild limit. In particular, in this limit, the special Kerr mode becomes a 
TTMl for s — —2; furthermore, the special mode and the TTM^j cancel each other 
for s — 2. Studying the limit a — > in detail, Maassen van den Brink reached a 
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second conclusion: the Schwarzschild special frequency f2; should be a limit point for 
a multiplet of "standard" Kerr QNMs, which for small a are well approximated by 

, 33078176 ma 3492608 

700009 2M 41177 4M2 ^ ' v a ; 
when I — 2, and by his Eq. (7.33) when I > 2. Numerical studies found QNMs 
close to the AS frequency, but not in agreement with this analytical prediction 
|231| . It was suggested (note [46] in Ref. [745 J) that QNMs corresponding to the 
AS frequency with to > may have one of the following three behaviors in the 
Schwarzschild limit: they may merge with those having to < at a frequency fij 
such that < (but > for I > 4) and disappear, as suggested by 
Onozawa [230] : they may terminate at some (finite) small a; or, finally, they may 
disappear towards uj — — ioo. Recently another alternative was suggested [747j, that 
in the Schwarzschild geometry a pair of "unconventional damped modes" should exist. 
For I = 2 these modes were identified by a fitting procedure to be located on the 
unphysical sheet lying behind the branch cut (hence the name "unconventional") at 
2Mw± = =f0.027+ (0.0033 — 4)i. An approximate analytical calculation confirmed the 
presence of these modes, yielding 2Mu!+ ~ —0.03248 + (0.003436 — 4)i, in reasonable 
agreement with the above prediction. If the prediction is true, an additional QNM 
multiplet should emerge near fli as a increases. This multiplet "may well be due to 
uj± splitting (since spherical symmetry is broken) and moving through the negative 
imaginary axis as a is tuned" 747 . The emergence of such multiplets was shown in 
Ref. |231| , but these do not seem to behave exactly as predicted |747j . 
Instability of naked singularities AS modes play an important role in the stability 
analysis of certain spacetimes containing naked singularities. An example of such a 
spacetime is the negative-mass Schwarzschild geometry, Eq. ^ with L — > oo and 
M < 0. From the general solution (jA.7p with C_2 = it follows that 

= C_ix- = C_ir (6M +{l- 1){1 + 2)ry^ e^^^''* , (A.9) 

where r, = (r + 2Aflog[— 2A/ + r]), is a regular Zerilli or scalar-type gravitational 
mode in the entire spacetime with frequency uj — iX{X + 1)/(3A/). Since ujj > 0, the 
spacetime is unstable. This instability was first found numerically by Gleiser and Dotti 
[748j and recognized to be an "algebraically special instability" in Ref. [749]. In fact, 
the calculation above can be extended to (negative-mass) charged BHs in de Sitter 
spacetimes [749] . The AS mode is not unstable for negative-mass Kerr geometries, but 
numerical results show that both negative-mass and overspinning "Kerr" geometries 
are unstable [750117511 1551 1555] . 
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Appendix B. 

Highly damped modes of Kerr black holes 

Here we reproduce the formulas from Kao and Tomino [222] , that allow one to compute 
analytically the highly damped modes of the Kerr spacetime, and which have been used 
to compare against numerical results in Fig. [TOl Define ui — a(A^/^ — A~^/^)/2\/3, 
VI = a(Ai/3 + A~i/3)/2 and A = 3\/3M/a + ^/l + 27AP/a^. Define also 

/o(ro) ^ ^rJ3ul + vf + 2^u,v,E [Tp] - %=^ilgi±ilx [To] 



\/'Sul + vl + 2iuiVi 



i2ro{3ui + ivi){ro — ui + ivi) 

\/3ul + vl + 2iuivi 
i{3ui + ivi){3ul -vl- 2rg) 



\/3u1 + vf + 2iuiVi 



n 



n[T_i,To] 
-2ivi 



3ui 



IVl 



To 



(B.l) 



with T_i 



-2ivi (ro+2Mi) 
/m(''o) 



, To 



4mi vi 



and 



{F[sin-iTi,T2] -i^[T2]} 



ro(ro + 2ul)^y3ul + vf - 2iuiVi 

X {n[T3,sin-iTi,T2] -n[T3,T2]} 



(B.2) 



Here Ti 



3u(+v'i+2tuivi ' -'-2 — 3u(+v'i-2iuivi J- 3 — (ro+2ui )(«i -iwi ) ' ^ '^^ 

functions E{m), E{ip^m)^ K{m), F[ip,rn),\\{n^m) and n(n, 93, m) are the elliptical 
integrals: 



E{m) 
K{m) 
n(n, m) 



Vl - m sin^ 6 dd, E{ip, m) = / Vl - m sin^ 6* ^6* 



VI 



■.dd,F{(p,m) = 



V 1 — ™ sin^ 6* 



[l-nsin^ 6l]-i[l -msin^ 0]-^/^ d9, 



rv 

n{n,ip,m)= / [l-nsm^9]-^[l~msin^d]-^^^d0. 
Jo 

Because of the pole at r = ro and branch cuts, there is an ambiguity in Eqs. (jB.ip 
and (|B.2p . The correct analytic expression is 

' fo{r+) - fo{r-) , iA7rMr+^ 



So 



-2i 



-8(3Mir+ — — Vi) 



(B.3) 



— 2ivi (r^+2ui) iiuivi 
{r^—u\—iv\){3u\~iv\) ' 3uj^-\-v'f+2iu\v\ 



where Q{x) is the step function. The term with a step function is introduced to com 
pensate the discontinuity caused by the term H 
Similarly, 

6,n = 4iMa [r+/„i(r+,ui,t)i) - r_/„i(r_, ui, wi)] /(r+ - r_). (B.4) 

Note that here there is no discontinuity and we do not need to introduce any term 
with step function. The quantities tu, 5loi, defined in Eqs. and (jlOip . are related 
to 5o, 5m via vj = 5m/ 5q , 5ivi = 2tt/6o- 
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Appendix C. 

Supermassive black hole binary candidates 

At the moment of writing there are three plausible SMBH binary candidates at close 
separation. Due to the relevance of these observations for SMBH mergers and their 
rates (see e.g. j585j ). here we briefly summarize the observations and their current 
interpretations . 

SDSSJ092712. 65+294344-0 - This quasar shows two sets of narrow emission lines, 
only one of which has associated broad lines, separated in velocity by 2650 kms~^, 
as well as additional emission and absorption lines at intermediate redshift. Komossa 
and collaborators |752j interpret the velocity separation between the two sets of lines 
as evidence for a recoiling BH with mass M « IQ^-^Mq- This large recoil speed can 
only be achieved if the final BH formed from the merger of two large-spin BHs in the 
so-called "superkick" configuration |7531 17541 1755j . The recoiling BH interpretation 
has been criticized by various authors. Bogdanovic et al. [756] propose a model where 
the blueshifted narrow lines originate from an accretion stream within the inner rim 
of the circumbinary disk of a massive BH binary with mass ratio g « 0.1 and mass 
M2 ~ IO^Mq for the secondary. Dotti et al. [75 7j also propose a model with a 
massive BH binary embedded in a circumbinary disk, where the blueshifted lines 
originate from gas swirling around the secondary BH. The BH binary has mass ratio 
q « 0.3, Ml Ri 2 X 1O^M0, a semi-major axis a « 0.34 pc and an orbital period 
P « 370 yrs. More detailed observations seem to favor the idea that the system 
represents the superposition of two AGN, rather than a recoiling SMBH |758j . Other 
interpretations suggest that the system is a more distant analog of NGC1275 [759], a 
large and a small galaxy interacting near the center of a rich cluster ^760] . 
SDSSJ153636. 22+044127.0 - This quasar shows two broad-line emission systems, 
separated in velocity by 3500 km/s, and unresolved absorption lines with intermediate 
velocity. Boroson and Lauer 584 interpret this quasar as a binary system of two 
SMBHs with masses Mi « lO*^-^Af0 and M2 ~ IO'^-^A/q (hence q « 40) separated by 
0.1 pc, with an orbital period of 100 yrs. Depending on unknown geometrical 
factors of the orbit, this system could coalesce either in 3 x 10^^ yrs or in 7 x 10^ yrs. 
Several alternative interpretations of the observations have been proposed. Gaskell 
|761j points out that the blueshift/redshift of the broad emission lines in this system 
can be interpreted in terms of normal line emission from a disk in an AGN. Radio 
observations by Wrobel and Laor |762j identify two faint, spatially distinct radio 
sources within the quasar's optical localization region, and suggest that the system 
could consist of a quasar pair separated by ~ 5 kpc. Chornock et al. [763[ I764j find a 
third broad emission component in the system, note a lack in velocity shift between the 
blueshifted and redshifted components, and argue that the system may be an unusual 
member of a class of AGNs known as "double-peaked emitters," rather than a SMBH 
binary or a quasar pair. Lauer and Boroson |765j remark that the lack of amplitude 
variations is unusual for a double-peaked emitter, and that longer temporal baselines 
are required to rule out the binary hypothesis. 

OJ287 - This BL Lac object shows quasi-periodical optical outbursts at 12-year 
intervals, with two outburst peaks (separated by one to two years) per interval. 
Optical observations of this source date back to 1890. Valtonen and collaborators 
interpret the two outbursts as happening when a smaller BH pierces the accretion 
disk of the primary BH, producing two impact flashes per period (see [766) for 
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details of the model). A model with non-spinning BHs of mass Mi — 1.8 x 
M2 = 1 X IO^Mq, semi-major axis a ~ 0.045 pc and eccentricity e ~ 0.66 was able 
to predict the date of the latest outburst. Remarkably, this model only reproduces 
the observations if gravitational radiation reaction is included |583j . Indeed, if 
the interpretation is correct, the system would inspiral very quickly and merge in 
~ 10^ years. Alternative models attribute the observed behavior to (1) oscillations in 
the accretion disk or jet of a single BH, or (2) variations of the accretion rate in a disk 
or a wobble of a jet in a binary BH (see |576| I583j for more details). 

If the binary BH interpretation is correct, the orbital parameters of the two 
candidate BH binaries from Refs. |583[I584] may not be too dissimilar: in gravitational 
wave lingo, both systems would be IMRIs (intermediate mass ratio inspirals). 

Appendix D. 

Black hole spin estimates 

Estimating BH masses is relatively easy because mass has a measurable effect even at 
large radii, where Newtonian gravity applies. Spin, on the contrary, is only measurable 
by looking at orbits close to the BH, where relativistic effects are important. The 
reason is that the gravitational potential around a rotating object can be expanded as 

$(r,6l) = ~{M /r) - q{M/rf P2icos9) + Oir-^) . (D.l) 

The parameter q = Q/M^ is a measure of the mass quadrupole moment, and to lowest 
order it scales quadratically with the spins |767[I768] : for BHs, Q — ~Ma?. Since the 
leading-order spin contribution scales like (M/r)^, spin measurements must rely on 
observations of test particles on orbits with very small radii. 

Luckily, we do have a chance to observe such orbits by looking at accreting gas 
close to the IS CO allowed by general relativity. All methods to measure spin using 
electromagnetic observations are based on variants of this idea, which is illustrated 
in the left panel of Fig. [151 There we show the variation of the ISCO frequency 
-^^isco and of the light-ring frequency Af IIlr with the dimcnsionlcss spin parameter 
a/M [160j . Positive values of a/M correspond to corotating orbits, and negative 
values correspond to counterrotating orbits. The orbital frequency at the ISCO 
has a maximum when a/M = 1 (i?isco = M) and a minimum when a/M = —1 
(^isco = 9M); for a Schwarzschild BH, i?isco = 6M and Mf^isco = 6-3/2. j.^^ 
gas in an accretion disk spirals in through a sequence of quasi-circular orbits as it 
viscously loses angular momentum; when it reaches the ISCO, it accelerates radially 
and falls into the BH, so the ISCO can be thought of as the inner edge of the accretion 
disk. The radiation efficiency 77, also plotted in Fig. [151 yields the energy radiated by 
the accretion disk per unit accreted mass. In principle 77 is determined by the binding 
energy of gas at the ISCO, which depends only on a/M. 

Most spin estimates are based on the observation of electromagnetic emission 
from accreting gas close to the ISCO. Unfortunately, all of these estimates are to some 
degree model-dependent. Some corrections to the highly idealized scenario described 
above are discussed in a number of papers [618[l769Tl770| . For example, Beckwith and 
Hawley |770j claim that an accurate modeling of magnetic fields would significantly 
displace the "inner edge" of the accretion disk from the standard general relativistic 
ISCO, and that this would lower the spin estimates inferred from Fe Ka observations. 
More in general, systematic errors originate from our limited understanding of i) the 
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physical processes producing the X-ray spectra , and ii) accretion physics in strong- 
field gravity regions. Because of these uncertainties, different models are sometimes 
capable of fitting the observations while predicting different values of the spin (see e.g. 
Ref. [35] for examples). A discussion of accretion disk physics as applied to a study 
of the strong-field gravity region around BHs is outside the scope of this paper, and 
we refer the reader to the many excellent review articles on this topic Q. 

Schematically, we can classify the main techniques used so far to estimate BH 
spins as belonging to four groups [SI I464|, [32] : 

(a) Continuum spectroscopy of accretion disks has been applied to various stellar-mass 
BH candidates. This method was pioneered by Zhang et al. j618] to suggest that 
two galactic superluminal jet sources, GRO J1655-40 and GRS 1915-1-105, should 
harbor rapidly spinning BHs. Narayan, McClintock and collaborators have embarked 
in a program to estimate spins for about a dozen BH candidates [614, ,619] I778|, 1779] 
WM mH inn] (but other groups are also very active [ST71 jWl HH [ZS2])- They 
also proposed that accretion can be used to provide hints of the presence of an event 
horizon for galactic stellar-mass BH binaries [32l [8] (but see [783] for criticism) . 

The basic idea is that when BHs have a large mass accretion rate the accreting gas 
tends to be optically thick, radiating approximately as a blackbody. In this spectral 
state (known as the "high soft state") the flux of radiation F{R) emitted by the 
accretion disk can be computed [784] and used to obtain an effective temperature 
profile T{r) ^ [F{R)/a]^^^, where a is the Stefan-Boltzmann constant. By comparing 
the blackbody radiation with the spectral flux received at Earth one can estimate 
the quantity Ri^cosi/D, where i?in is the inner edge of the accretion disk, D is the 
distance to the source and l the inclination angle. Unfortunately the method relies 
on independent estimates of t and D. Perhaps the major weakness of the method 
consists in the fact that, in practice, the observed spectrum deviates significantly 
from a blackbody. These deviations are usually modeled by artificially increasing the 
temperature of the emitted radiation by a poorly known "spectral hardening factor" , 
which is usually (and roughly) approximated by a constant /coi — 1.7. Perhaps the 
most precise spin measurements to date using this technique regard the eclipsing X- 
ray binary M33 X-7, with a claimed value a/M — 0.77 ± 0.05 [616] . and the first 
extragalactic X-ray binary LMC X-1, for which a/M = 0.90^°;°^ [BID]. 

(b) Spectroscopy of relativistically broadened Fe Ka fluorescence lines has been 
proposed as a promising alternative to continuum fitting. This method originated 
from the discovery of a strong, broad spectral line in the X-ray spectrum of the 
Seyfert 1.2 galaxy MCG-6-30-15 [785]. Brenneman and Reynolds [^ applied the 
method to XMM-Newton observations of this system, estimating the spin to be very 
near the maximal limit: a/M — 0.9891q qq2- ^^is case, the largest source of error 
comes from the unknown dependence of the line emissivity on the disk radius i?, 
which is usually modeled as a power law ^ . An analysis of systematic errors involved 
in the Fe Ka measurements was carried out by Reynolds and Fabian [769] ■ Their 

* The most pessimistic viewpoint on electromagnetic spin measurements is perhaps that of Miller 
and Turner: according to them, "with the uncertainty and debate on the physical processes that 
dominate the X-ray spectrum in this energy range [...] it is not currently possible to constrain the 
BH spin in a model-independent way" I771| (see also |772l '). 

tt Liedahl and Torres [773) focus on the basic general relativistic equations, scattering processes 
and atomic transitions. Reynolds and Nowak 17741 discuss theoretical problems in the modeling of 
relativistic disk lines. Fabian and Miniutti |775) and Miller |776| review relativistic X-ray emission 
lines from the inner accretion disks around rotating BHs and observational prospects for constraining 
the spin history of SMBHs (see also }777) for a similar review focusing on stellar-mass BHs). 
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main finding is that systematic errors can be significant for modest values of the 
spin, but they decrease for large spins, so large-spin measurements (such as the one 
of Ref. [607j ) should be more reliable. The method has been applied to other AGN 
[7551 [7571 [7551 [755] and even to stellar mass BH candidates [777] . A very accurate 
measurement of a/M = 0.935 ± 0.01(statistical)±0.01(systematic) has recently been 
claimed by Miller et al. [62l] and Reis et al. [iSO] for the stellar-mass BH GX 339-4 
in outburst. Ref. |612j discusses BH spin estimates for several different systems and 
their correlations with other physical properties of X-ray binaries. 

(c) Quasi-periodic oscillations (QPOs) are likely to offer the most reliable spin 
measurements for accreting systems that harbor stellar-mass BH candidates, once the 
correct model is known. Unfortunately, at present there are several competing models 
to explain QPOs (see |464l[32ir791j for a discussion and further references), and models 
based on different physical assumptions typically yield very different spin estimates. 
For example, Figure 12 in Psaltis's review [32] shows that the diskoseismology model 
of Ref. [792] and the parametric resonance model of Ref. [793] give very different BH 
spin estimates in the case of GRO J1655-40. At the very least, an indication that 
QPOs are related to dynamical frequencies near the ISCO comes from the fact that 
the QPO frequencies roughly scale with the inverse of the BH masses, as they should 
in general relativity. A relatively model-independent lower limit on the spin can then 
be obtained for GRO J1655-40: for this system one gets a lower limit of a/M > 0.25 
by requiring only that the observed 450 Hz oscillation frequency must be limited by 
the azimuthal frequency at the ISCO, a rather solid assumption |622j . 

(d) Statistical methods based on the radiative efficiency of AGN. A rather general 
argument constrains, in principle, the average properties of the spin distribution in 
AGN. The general idea is suggested by a glance at the efficiency rj = L^^c/Mc^ (i.e., 
the energy radiated iacc per unit accreted mass) in the right panel of Figure [T5| For 
a non-rotating BH 77 = 1 — (8/9)^^^ ~ 0.057, while for a maximally rotating Kerr BH 
77 ~ 0.42 is much larger. In a typical accretion system one can measure Lace (provided 
the distance is known), but not M, so there is no way of knowing rj accurately enough 
to estimate a/M . However, by observing high-redshift AGN one can estimate the 
mean energy radiated by SMBHs per unit volume. By considering SMBHs in nearby 
galaxies we can estimate the mean mass in SMBHs per unit volume of the current 
universe. If we assume that SMBHs grow mostly by accretion, by dividing these two 
quantities we can get the average radiative efficiency of AGN, hence their average spin 
(this is a variant of the famous "Soltan argument" [609] ). 

Elvis et al. [626j and Yu and Tremaine [627] use observational data to infer 
that the mean efficiency 77 ~ 0.10 — 0.15 on average, and ~ 0.2 (or larger) for 
the most massive systems. This is possible only if SMBHs have significant rotation. 
Their arguments have been supported by several later studies, some of which claim 
an average efficiency ~ 0.3 — 0.35 for quasars at moderate redshift [628|. I629[ 1630] . 
However there is no general consensus on the interpretation of the data. Shankar 
et al. [631] . Merloni and Heinz [632] and Daly [633] present arguments in favor of 
low average radiative efficiencies, hence low spins. Wang et al. [606] argue that the 
radiative efficiency has a strong cosmological evolution, decreasing from rj k, 0.3 at 
z ~ 2 to T] 0.03 at z « 0. The uncertainties in disentangling radiative efficiencies 
from quasar lifetimes were pointed out in Ref. [794] . Daly [610] provides spin estimates 
for 19 powerful FRII radio sources and for 29 central dominant galaxies (CDGs). For 
the first class of sources the spins seem to decrease from near-extremal values at z = 2 
to ~ 0.7 at z = 0, while for the (lower power) CDGs the estimated spins are in the 
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range 0.1 — 0.8. In conclusion, the jury is still out on the mean radiative efficiency 
(and mean spin) of quasars. The outcome of this debate is fundamental to constrain 
SMBH evolution models [482 1160111795] . 

(e) Other methods. Besides the methods listed above, there are other avenues for 
measuring BH spin that hold promise for the future. One of these is polarimetry. 
The idea is to exploit the fact that, in general relativity, the plane of polarization 
of BH disk radiation changes with energy. This is a purely relativistic effect, absent 
in Newtonian gravity, and the magnitude of the change of the plane of polarization 
can give a direct measure oi a/M |7961 1464j . A second idea is based on Very Long 
Baseline Interferometry (VLBI) imaging of the silhouette of the SMBH in M87. This 
system is at a distance of 16 Mpc, much farther away than the SMBH at the galactic 
center; however, because of its larger mass of ~ 3.4 x lO^M©, the apparent diameter 
of M87's horizon is 22 fj,as, about half as large as Sgr A*. Unhke Sgr A*, M87 exhibits 
a powerful radio jet, hence it holds promise for exploring the relation between the BH 
spin and the jet generation mechanism |797j . A third, indirect way to constrain BH 
spins is based on energetic considerations. The AGN outburst in the MS0735. 6+7421 
cluster's central galaxy implies that its putative SMBH grew by about 1/3 of its 
mass in the past 100 Myr, accreting matter at ~ 3 — ^Mq/yi, inconsistent with the 
Bondi mechanism. The energetics of the system could be explained instead by angular 
momentum released from a rapidly spinning SMBH with M > lO^^M© |611| . 
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